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ABSTRACT 


A  numerical  investigation  of  two-dimensional  natural  convection  flow 
and  heat  transfer  from  a  substrate-mounted  flush  heat  source  immersed  in 
a  liquid-filled  square  enclosure  was  conducted.  The  study  is  relevant  to 
direct  liquid-immersion  cooling  of  electronic  components.  A  control 
volume  based  finite-difference  model  that  accounts  for  conduction  heat 
transfer  within  the  substrate  and  heat  source  and  the  coupled  natural 
convection  in  the  fluid  was  utilized.  Numerical  predictions  were  obtained 
for  a  wide  range  of  Rayleigh  and  Prandtl  numbers,  substrate  to  fluid  and 
heat  source  to  fluid  thermal  conductivity  ratios  and  other  geometrical 
parameters  that  may  be  encountered  in  practice.  An  increase  in  Rayleigh 
number  lead  to  more  vigorous  flow  and  promoted  cooling.  No  noticeable 
effect  on  the  nondimensional  temperatures  was  observed  when  changing 
the  Prandtl  number  from  7  to  100.  Little  reduction  in  maximum 
temperatures  was  observed  when  substrate  and  component  to  fluid  thermal 
conductivity  ratios  were  increased  beyond  10  and  25,  respectively. 
Component  to  substrate  width  ratio  change  from  .25  to  .999  resulted  in 
approximately  linear  decrease  in  the  maximum  temperature.  A 
companion  experimental  study  of  three-dimensional  natural  convection 
transport  from  a  flush  mounted  array  of  heat  sources  in  water  was  also 
conducted.  Computed  temperatures  compared  favorably  to  appropriate 
experimental  data. 
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NOMENCLATURE 


Cp  specific  heat  at  constant  pressure  [J/kg-K] 

d  dimension  in  Fig.  1  [m] 

dt  substrate  thickness  [m] 

g  gravitational  acceleration  [m/s2] 

h  component  height  [m] 

H  enclosure  height  [m] 

k  thermal  conductivity  [W/m-K] 

L*  dimension  in  Fig.  1  [m] 

p  pressure  [N/m2] 

P  non-dimensional  pressure,  p/pU02 

Pr  Prandtl  number,  pcp/kf 

q  non-dimensional  heat  flux  [eq.7] 

Q  Heat  generation  rate  per  unit  length  [W/m] 

P  Power  to  the  heater 

Rc  protrusion  to  fluid  thermal  conductivity  ratio,  k</kf 
R„  substrate  to  fluid  thermal  conductivity  ratio,  k./kf 

Ra  Rayleigh  number,  gPQh3/ak^v 

51  non-dimensional  counter  clockwise  contour  distance  along  the 

component-fluid/component-substrate  interface  in  Table  2 

52  non-dimensional  counter  clockwise  contour  distance  along  the 

solid-fluid  interface  in  Table  1 
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Sb,  St  dimensions  in  Fig.  1  [m] 

Rl  Lead  resistance  [  Q  ] 

RP  Precision  resistor  resistance  [  ft  ] 

tfc  enclosure  wall  temperature  [K) 

T  non-dimensional  temperature,  (t-tc)/(Q/kf) 

u  vertical  velocity  component  [m/s] 

U  non-dimensional  vertical  velocity  component,  u/LJ0 

UD  reference  velocity,  (gPQh/kf)172  [m/s] 

v  horizontal  velocity  component  [m/s] 

V  non-dimensional  horizontal  velocity  component,  v/U0 

VL  Voltage  measured  over  the  leads  and  heater  [Volts] 

V-p  Voltage  over  the  precision  resistor,  leads  and  heater  [Volts] 
w  component  width  [m] 

x  vertical  coordinate  [m] 

X  non-dimensional  vertical  coordinate,  x/h 

y  horizontal  coordinate  [m] 

Y  non-dimensional  horizontal  coordinate,  y/h 

Greek  Symbols 

a  fluid  thermal  diffusivity  [m2/s] 

P  coefficient  of  thermal  expansion  [1/K] 

8X  non-dimensional  vertical  distance  between  grid  point  and  the 
control  volume  face 

6Y  non-dimensional  horizontal  distance  between  grid  point  and  the 
control  volume  face 
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p 

p 


fluid  density  [kg/m3] 
dynamic  viscosity  [kg/m-s] 
v  kinematic  viscosity  [m2/s] 

k  nondimensional  conductivity  [R/(RaPr)1/2] 

Subscripts 

c  component  (chip) 

f  fluid/liquid 

max  maximum 

s  substrate 

«ur  substrate  or  component  surface 
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L  INTRODUCTION 


A.  SEAT 


DU',  I W 


NT  OF  PROBLEM 


Heat  removal  from  electronic  equipment  is  an  area  of  extensive  research 
and  will  be  more  important  in  the  future  due  to  the  ever  increasing 
component  volumetric  heat  generation  rates.  The  advantages  of  liquid 
cooling  have  been  demonstrated  in  the  thermal  control  of  macroelectronic 
components  such  as  power  transistors  and  power  supplies  for  almost  40 
years.  However  liquid  cooling  of  microelectronic  equipment  has  only 
recently  gained  wide  attention  due  to  the  need  for  increased  power 
dissipation,  the  availability  of  inert  dielectric  liquids,  and  the  introduction 
of  immersion  cooled  supercomputers.  Bergles  and  Bar-Cohen  [Ref.  1]  and 
Nakayama  [Ref.  2]  have  described  various  cooling  techniques  currently 
under  investigation  to  meet  modem  component  cooling  requirements. 
Incropera  [Ref.  3]  discusses  convective  cooling. 

In  natural  convection  the  flow  is  induced  by  buoyancy  forces.  The 
buoyancy  force  is  generated  when  a  body  force  acts  on  a  fluid  in  which  a 
density  gradient  is  present.  Natural  convection  cooling  has  the  advantages 
of  simplicity  of  design,  low  operating  cost  due  to  no  outside  power 
requirements  and  minimum  maintenance,  absense  of  noise  and  high 
reliability.  Also,  in  the  event  of  a  mechanical  failure,  natural  convection 
may  be  the  only  method  of  cooling  available  [Ref.  4-5].  Natural  convection  in 
liquids  results  in  much  higher  cooling  rates  compared  to  those  of  air  and 
hence  numerous  studies  have  recently  been  performed  on  natural 
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convection  from  discrete  heat  sources  inside  liquid  filled  enclosures  [Ref.  6- 

103. 

Gaiser  [Ref.  6]  experimentally  investigated  natural  convection  liquid 
cooling  of  a  three  column  array  of  uniformly  heat  generating  simulated 
electronic  components  on  a  flat  plate  and  in  a  vertical  channel  by  using  foil 
heaters.  The  channel  was  made  by  placing  a  smooth  movable  shroud 
parallel  to  the  high  density  heater  surface.  Temperature  measurements 
and  flow  visualizations  at  various  power  levels  and  channel  spacings  were 
made.  Experimental  correlations  relating  local  heat  transfer  with  shroud 
wall  spacing  and  foil  heater  power  levels  were  developed.  Flow  was 
laminar  for  all  power  levels  studied.  For  the  unshrouded  case,  the  flow  of 
coolant  over  the  center  column  resulted  in  a  weak  fluid  entrainment  from 
the  side  columns  and  in  turn  enhanced  the  cooling  of  the  side  columns. 
For  the  shrouded  case,  the  fluid  entrainment  by  the  center  column  was 
much  greater  and  had  a  greater  cooling  effect  on  the  side  columns.  For  a 
channel  of  3  mm  or  less,  the  increased  frictional  resistance  caused  a 
significant  increase  in  the  component  temperature  and  the  correlations 
were  no  longer  followed. 

Joshi,  Willson  and  Hazard  [Ref.  7]  conducted  an  experimental 
investigation  of  steady  state  and  transient  natural  convection  from  a 
column  of  rectangular  heated  protrusions  in  a  vertical  channel  in  water. 
Flow  visualizations  and  component  surface  temperature  measurements 
were  carrried  out  for  several  power  dissipation  levels  and  channel 
spacings.  For  the  smallest  channel  spacing,  the  component  surface 
temperatures  increased  significantly  due  to  a  reduction  in  fluid  velocity. 
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During  transient  periods,  an  initial  diffusive  transport  was  observed 
followed  by  the  evolution  of  convective  effects.  No  overshoots  in  component 
temperatures  were  observed.  Steady  transport  responses  over  several 
component  power  levels  were  also  examined. 

Park  and  Bergles  [Ref.  8]  simulated  microelectronic  circuits  with  flush 
mounted  and  protruding  thin  foil  heaters.  Heat  transfer  coefficients  were 
obtained  for  various  heater  heights  and  widths  in  water  and  R-113.  The 
heat  transfer  coefficient  was  found  to  increase  with  decreasing  width, 
which  was  attributed  to  three  dimensional  boundary  layer  effects. 
Protruding  heaters  were  found  to  have  a  coefficient  about  15  percent  above 
flush  mounted  heaters.  Coefficients  for  upper  flush  heaters  were  found  to 
be  less  than  lower  heaters.  For  protruding  heaters,  the  upper  heaters  had 
higher  coefficients  than  the  lower  heaters. 

Bar«Cohen  and  Schweitzer  [Ref.  9]  provided  an  analytical  basis  for  the 
design  and  optimization  of  convective  immersion  cooling  systems. 
Analytical  development  and  experimental  verification  of  the  relations  for 
the  natural  convection  heat  transfer  coefficients  prevailing  along  the 
surfaces  of  immersed,  uniformly  heated  plates  in  both  symmetric  and 
asymmetric  configurations  were  conducted.  Milanez  and  Bergles  [Ref.  10] 
experimentally  studied  the  wake  temperature  distribution  produced  by  a 
horizontal  line  source  on  a  vertical  adiabatic  surface  and  compared  it  to  the 
similarity  solution.  They  also  examined  the  interaction  of  two  heaters. 
Velocities  induced  by  the  lower  heater  enhanced  the  heat  transfer 
coefficient  at  the  upper  heater,  however  this  effect  was  offset  by  the 
preheating  of  the  fluid  by  the  lower  heater.  Jaluria  [Ref.  11]  conducted  a 
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numerical  study  of  heat  sources  mounted  on  an  adiabatic  wall  under  the 
influence  of  natural  convection.  Sathe  and  Joshi  [Ref.  5]  conducted  a 
numerical  investigation  of  natural  convection  flow  and  heat  transfer 
arising  from  a  substrate  mounted  protruding  heat  source.  Their  study 
accounts  for  conduction  heat  transfer  within  the  protrusion  and  substrate 
and  the  coupled  natural  convection  in  the  fluid. 

Areas  of  current  concern  with  the  new  packaging  techniques  are  stress 
on  the  electronic  components  due  to  thermal  expansion  coefficient 
mismatch.  Temperatures  throughout  the  package  are  needed  so  thermal 
stress  analysis  can  be  conducted  and  high  stress  areas  can  be  located.  The 
complex  geometries  of  the  boards  and  the  coupled  conduction/convection 
heat  transfer  problem  lead  to  intricate  heat  transfer  paths.  The  difficulty  is 
magnified  when  the  circuit  board  spacing  is  small  and  the  flow  pattern  is 
populated  by  numerous  non-uniform  heat  sources.  An  examination  of  the 
heat  transfer  and  fluid  flow  characteristics  of  such  a  system  would  aid 
designers  in  the  development  of  better  cooling  schemes. 

B.  OBJECTIVES 

This  study  was  performed  to  determine  the  effects  of  changing  the 
various  parameters  involved  in  the  construction  of  electronic  packages  on 
the  heat  transfer  and  to  observe  the  effects  of  a  nonuniform  high  density 
grouping  of  discrete  heat  sources  in  a  vertical  channel  for  direct  liquid 
natural  convection  cooling.  Gaiser  [Ref.  6]  describes  the  high  density 
grouping  of  heat  sources. 
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The  study  consisted  of  two  parts.  First  a  computer  model  was  utilized  to 
conduct  a  two  dimensional  study  to  find  the  effects  of  various  governing 
parameters  on  the  natural  convective  heat  transfer  and  fluid  flow  behaviors 
for  a  substrate  mounted  flush  component  immersed  in  a  liquid.  Then  the 
three  dimensional  flow  over  a  simulated  printed  wire  circuit  board 
consisting  of  multiple  arrays  of  flush-mounted  heaters  was  observed 
experimentally.  Results  of  the  computer  model  were  compared  to  results 
from  past  and  current  experimental  studies. 

Specific  objectives  of  the  numerical  parametric  study  were: 

-  To  determine  temperatures  and  flows  throughout  and  around  a 
simulated  direct  liquid  natural  convection  cooled  flush  mounted  electronic 
component. 

-  To  study  the  effect  of  various  combinations  of  parameters  for  an 
electronic  component  flush  mounted  on  a  substrate. 

-  To  investigate  the  effect  of  two  component  interaction  on  the  heat 
transfer  and  flow. 

The  experimental  study  had  the  following  objectives: 

-  To  determine  the  effect  of  channel  spacing  on  fluid  flow  and  heat 
transfer. 
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-  To  observe  the  effect  of  nonuniform  powering  on  heat  transfer  and  fluid 
flow. 


-  To  observe  the  surface  temperature  patterns  during  uniform  and 
nonuniform  heating  using  liquid  crystal. 
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II.  NUMERICAL  MODEL 


A.  MATHEMATICAL  FORMULATION 

The  schematic  diagram  of  the  configuration  examined  is  shown  in  Fig. 
1.  A  flush  heat  source  mounted  in  a  vertical  substrate  is  immersed  in  a 
two-dimensional  fluid-filled  square  enclosure  of  height  H.  The  dark  shaded 
region  represents  the  heat  source  and  the  light  shaded  area  the  substrate. 
Uniform  volumetric  heat  generation  takes  place  in  the  flush  mounted 
heater.  The  enclosure  boundaries  are  maintained  at  a  constant 
temperature,  tc.  The  heater,  the  substrate  and  the  fluid  have  constant  but 

different  thermophysical  properties.  Assuming  a  steady-state,  laminar 
flow  with  no  viscous  dissipation  and  the  Boussinesq  approximation  to  be 
true,  the  dimensionless  governing  equations  can  be  written  as: 


Fluid  Region 


au  av 
ax  dY 


=o 


(i) 


d(U2)  d(UV) 

ax  dY 


.(WW)w|«+&4T  * 

ax2  3Y2  oX 


(2) 


aqiv) .  d(v2) 
ax  dY 


d(UT)  d(VT) 
ax  3Y 


(3) 

(4) 
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Figure  1.  Schematic  diagram  of  configuration  for  numerical  study 


:U4LL  MEL-  ? 


rnm&grw* 


Solid  Region  (heater) 

a^r  la2Tr  h 

ax2  3Y2  “  wRc  ®) 


Solid  Region  (substrate) 


ax2  +  3Y2 


=  0 


(6) 


where  X=x/h,  Y=y/h,  U=u/U0,  V=v/U0,  T=(t-tc)/(Q/kf),  P=p/pU02; 
U0=<gPQh/kf)1/2,  Pr=pCp/kf,  Ra=gPQh3/akfv,  R^Vkf  and  R,=k^kf. 

The  boundary  conditions  constitute  isothermal  enclosure  walls  at 
temperature  tc  and  the  no-slip  and  impermeable  wall  conditions  for  the 

velocity  components;  i.e  U=V=T=0  at  the  enclosure  walls.  Heat  fluxes  are 
appropriately  matched  at  the  interfaces  of  dissimilar  materials.  The 
following  governing  parameters  emerge  as  a  result  of  the  non- 
dimensionalized  governing  equations,  boundary  and  matching  conditions: 
Ra,  Pr,  Rc,  R,,  St/h,  d/h,  Lj/h,  Sj/h  and  w/d*.  The  definitions  of  the  symbols 

can  be  found  in  the  Nomenclature. 


B.  METHOD  OF  SOLUTION 

The  governing  equations  are  discretized  using  a  finite  difference  scheme 
wherein  the  control  volumes  for  the  temperature  and  pressure  are 
staggered  from  those  for  the  velocities.  Power  law  profiles  are  used  for  the 
spatial  variation  of  the  dependent  variables  to  ensure  realistic  results  for  a 


wide  range  of  the  grid  Peclet  numbers.  Interface  diffusivities  are 
calculated  using  a  harmonic-mean  formulation  in  order  to  handle  abrupt 
changes  in  the  material  properties.  The  details  of  the  discretization  process 
can  be  found  in  Patankar  [Ref.  12]. 

The  discretized  equations  were  solved  iteratively  using  the  line-by-line 
TDMA  (Tri-Diagonal  Matrix  Algorithm)  and  the  SIMPLER  procedure  as 
outlined  in  Ref.  12.  It  is  noted  that  even  though  separate  equations  are 
written  for  the  fluid  and  solid  regions,  the  solid  region  is  numerically 
simulated  by  letting  its  viscosity  become  very  large.  Thus  the  same 
momentum  equation  is  solved  throughout  the  computational  domain. 
Similarly,  for  the  energy  equation,  heat  sources  and  property  values  are 
used  implicitly  in  the  control  volume  formulation  wherein  mass, 
momentum  and  energy  balances  are  incorporated  for  individual  control 
volumes  [Ref.  12]. 

Test  computations  were  performed  on  a  series  of  grids  ranging  from 
30x30  to  73x73  control  volumes  to  determine  the  grid  size  effects.  The  values 
and  locations  for  the  maximum  temperature  and  peak  velocities  did  not 
change  appreciably  when  the  grid  was  refined  beyond  40x40  control 
volumes.  Most  calculations  were  performed  on  a  grid  size  of  50x50. 
However  for  Ra  >  106  the  solution  would  not  converge  on  a  50x50  grid  and  it 
was  found  that  by  increasing  the  control  volumes  in  the  area  of  the 
suspected  high  thermal  gradient,  convergence  could  be  achieved. 
Therefore  for  Ra  =  106  the  number  of  control  volumes  was  increased  to 
70x70.  Convergence  was  based  on  a  balance  of  the  rate  of  energy  generated 


10 


in  the  component  and  the  rate  of  energy  leaving  the  enclosure  walls. 
Convergence  was  based  on  the  energy  balance  being  better  than  1  percent. 

The  program  was  run  on  the  Mechanical  Engineering  VAX  station  2000 
Cluster  and  the  IBM-370/3033  computer.  The  greater  the  number  of  control 
volumes  the  the  greater  the  time  to  run  on  each  system.  Also  if  Ra  was  of 
the  order  of  106,  it  took  up  to  120,000  iterations  for  covergence  to  take  place. 
The  relaxation  factor  for  the  velocities,  temperatures,  pressure  and  buoyant 
force  were  adjusted  until  a  proper  combination  was  found  that  would  make 
the  solution  converge.  After  coming  up  with  a  proper  set  of  relaxation 
factors,  they  could  be  used  for  other  high  Ra  solutions  and  the  number  of 
iterations  required  was  reduced.  The  run  times  could  also  be  reduced 
when  previously  computed  solutions  were  used  as  an  initial  guess.  The 
VAX  system  took  approximately  10  times  longer  to  run  a  program  with  the 
same  number  of  iterations  and  control  volumes  than  the  IBM  system. 
However,  more  VAX  stations  could  be  used  on  a  continuous  basis  which 
made  using  the  VAX  system  slightly  preferable. 


HL  PARAMETRIC  STUDY 

Numerical  computations  were  conducted  for  various  Ra,  Pr,  R„  and  Rc 
values  that  would  be  found  in  actual  cooling  applications  [Ref.  5].  Ra 
numbers  from  10s  to  106  were  investigated  with  the  limit  of  106  due  to  the 
very  slow  convergence  of  the  computations  for  Ra  greater  than  106.  The 
values  of  Pr,  R,  and  Rc  were  varied  from  7  to  100,  0.1  to  100  and  1  to  100 

respectively.  The  geometric  parameters  used  in  the  computations  were 
St/h^,  d/h=2,  Lt/h=2,  St/h=3,  and  w/dt  varied  from  0.25  to  1.  A  baseline 
case  with  Ra=105,  Pr=25,  R,=10,  Rc=25  and  w/dt=0.5  was  conducted  and 
when  a  particular  parameter  was  being  studied  all  other  parameters  were 
kept  the  same  as  in  the  baseline  case.  Table  1  shows  the  breakdown  of  the 
cases  studied.  Appendix  A  shows  the  user  portion  of  the  finite  difference 
program  used  to  obtain  the  numerical  data.  Appendix  B  contains  the 
contour  program  used  in  making  the  following  figures. 

A.  EFFECT  OF  Ra 

The  temperature  contours  for  Ra=103  to  106are  shown  in  Figs.  2(a)-(d). 
These  are  nearly  symmetric  with  respect  to  the  enclosure  vertical 
centerline  for  all  Ra.  The  high  R,  (10)  value  tends  to  spread  the  heat  out  to 

each  side  even  though  the  heater  is  placed  on  the  right  side  of  the  substrate. 
As  Ra  increases,  the  extent  of  thermal  stratification  increases  in  the  lower 
portion  of  the  enclosure.  All  boundary  layers  get  thinner  and  a  larger 
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Table  1.  Baseline  parameters  (all  parameters  except  the  one  being  studied 
revert  to  these)  and  the  range  of  parameters 


Parameter 

Baseline 

value 

Range  studied 

Ra 

105 

103, 104, 105, 10e 

Pr 

25 

7, 25, 100 

R. 

10 

.1, 1, 10, 50, 100 

Rc 

25 

1, 25, 50, 100 

w/dt 

.5 

.25,  .5,  .75,  .999 

thermal  gradient  occurs  in  the  fluid  next  to  the  substrate  and  the  enclosure 
top  and  sides. 

Figures  2(a)-(d)  also  show  changes  in  the  streamline  patterns  with  more 
vigorous  flow  with  an  increase  of  Ra.  The  flow  travels  up  on  both  sides  of 
the  package.  The  flow  forms  two  cells  with  the  right  cell  being  somewhat 
larger  and  stronger,  due  to  the  direct  heater-fluid  interaction.  Flows  are 
nearly  symmetric  on  the  right  and  left  and  as  Ra  is  increased  the  eddy 
centers  move  up  and  toward  one  another.  Also  as  Ra  is  increased,  the  fluid 
in  the  lower  right  portion  of  the  enclosure  becomes  increasingly  stagnant. 

The  heat  loss  though  each  face  of  the  component  is  shown  in  Fig.  3.  As 
Ra  is  increased  the  percentage  of  the  heat  lost  directly  to  the  fluid  i.e.  from 
the  front  face,  is  increased  due  to  the  increased  convection.  As  Ra 
increases  the  percentage  of  heat  lost  from  the  top  and  bottom  of  the 
enclosure  decreases  as  shown  in  Fig.  4,  and  the  heat  lost  to  the  sides  of  the 
enclosure  increases.  The  heat  loss  through  the  enclosure  bottom  is  below  5 
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Figure  3.  Heat  loss  through  component  faces  for  various  Ka 


■  LEFT 
0  TOP 
E  RIGHT 
0  BOTTOM 


Figure  4.  Heat  loss  through  enclosure  walls  for  various  Ra 


percent  for  all  Ra  shown  but  it  decreases  with  an  increase  in  Ra  due  to 
increased  thermal  stratification  near  the  bottom.  There  is  a  greater 
percentage  of  heat  loss  through  the  left  wall  than  the  right  even  though  the 
right  wall  faces  the  component  and  the  flow  is  greater  on  the  right.  This  is 
because  the  left  wall  is  closer  to  the  substrate  and  heater  than  the  right 
wall.  It  is  interesting  to  note  that  the  top  wall  dissipates  almost  half  the 
heat  generated  by  the  component. 

Figure  5  shows  the  substrate  nondimensional  surface  heat  flux  (q) 
which  is  defined  as: 


Table  2  describes  the  values  of  S2  corresponding  to  its  location  on  the 
substrate.  Heat  fluxes  are  not  defined  at  the  comers.  S2=0  at  the  lower 
right  comer  of  the  substrate  and  goes  counter-clockwise  around  the  solid- 
fluid  boundary  including  the  component-fluid  boundary.  The  same  general 
trend  for  the  variation  of  heat  flux  with  S2  is  seen  for  all  Ra.  The  heat  flux 
goes  through  a  minumum  at  S2-0.5  after  which  it  rapidly  reaches  its 
maximum  value  at  S2=2,  at  the  base  of  the  component.  The  heat  flux  then 
decreases  going  up  the  component  and  when  the  substrate  above  the 
component  is  reached  at  S2=3  it  decreases  even  more  steeply,  going  through 
a  minimum  at  S2-4.5  prior  to  reaching  the  top  right  comer.  At  the  top 
right  corner  the  heat  flux  decreases  sharply  again  going  over  the  top 
substrate  face  and  comes  to  the  minimum  heat  flux  on  the  substrate  at  the 


Table  2.  Ranges  of  S2  corresponding  to  different  solid-fluid  interfaces 


Surface  no. 
in  Fig.  1 

X 

Y 

S2 

1 

0to2 

0 

0to2 

2 

2to3 

0 

2to3 

3 

3to5 

0 

3to5 

4 

5 

Oto-1 

5to6 

5 

5to0 

-1 

6  to  11 

6 

0 

-ItoO 

11  to  12 

midpoint  of  the  top  face  at  S2=5.5.  Reviewing  the  temperature  contours  in 
Figs.  2(aMd),  this  is  where  the  fluid  temperature  is  the  hottest  next  to  the 
substrate.  The  heat  flux  on  the  top  face  is  nearly  symmetric.  Proceeding 
down  the  substrate  left  face  the  heat  flux  goes  through  a  small  dip  but 
quickly  increases  as  one  approaches  the  back  of  the  component.  It  reaches 
a  local  maximum  heat  flux  at  the  back  of  the  base  of  the  component  at  S2  »  9 
and  then  decreases  to  a  local  minimum  at  S2-10.25.  It  then  increases  and 
comes  to  a  local  maximum  at  the  bottom  left  comer  at  S2=ll.  The  heat  flux 
on  the  bottom  face  shows  the  same  symmetry  as  the  top,  however  it  is  much 
higher  due  to  colder  fluid  in  the  vicinity  of  the  bottom  face. 

Nondimensional  substrate  surface  temperatures  (Tsur)  for  various  Ra 

are  shown  in  Fig.  6.  The  substrate  surface  temperature  is  computed  using 
the  harmonic  mean  formulation  as: 


T  = 

A  sur 


(*♦*) 


(8) 
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The  lowest  temperatures  are  at  the  bottom  of  the  substrate  and  the  highest 
are  at  the  top  of  the  component.  Local  maxima  are  evident  at  the  back  of  the 
substrate  due  to  direct  conduction  from  the  heater.  The  nondimensional 
temperatures  are  higher  for  the  lower  Ra.  The  higher  the  Ra,  the  greater 
the  buoyancy  induced  flow  and  hence  higher  heat  transfer  from  the 
component,  leading  to  lower  non-dimensional  temperatures.  Note  that  if  Q 
and  kf  are  held  fixed,  a  higher  T  implies  a  higher  actual  (dimensional) 

temperature.  In  such  a  case,  higher  Ra  can  be  obtained  by  using  a  liquid 
with  higher  p,  leading  to  higher  buoyant  forces,  enhanced  convection  and 
reduced  actual  temperatures.  Thus  the  dimensional  results  are  consistent 
with  the  non-dimensional  trends. 

Figure  7  shows  the  nondimensional  component  surface  heat  flux  for 
various  Ra  versus  SI.  Table  3  shows  the  corresponding  position  of  the 
component  to  SI  with  Sl=0  at  the  lower  right  comer  of  the  component  and 
increasing  going  counter-clockwise  around  the  component.  Fluxes  for  all 
Ra  show  similar  trends.  The  component  front  face,  in  direct  contact  with 
the  fluid,  has  the  lowest  heat  flux.  The  top  face  has  the  second  highest  heat 
flux  and  the  back  face  has  a  slightly  lower  heat  flux.  The  bottom  of  the 
component  has  generally  the  greatest  heat  flux.  The  largest  heat  flux  is 
near  the  bottom  right  comer. 

The  trends  of  Fig.  8  show  the  component  nondimensional  surface 
temperatures  to  be  similar  for  all  Ra.  Surface  temperatures  are  fairly 
uniform  around  the  component  with  the  highest  temperature  near  the  top 
on  the  front  face.  The  high  R,  and  R<.  values  allow  the  heat  to  be  uniformly 

spread  so  that  the  component  temperature  remains  fairly  uniform. 
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SI 

Interface’" 
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2  to3 

1 

3 

o  -.5 
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3 
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1  to  1.5 

C-8 

Back 

3to2 

5 

toO 

1.5  to  2.5 

C-8 

Bottom 

2 

-.51 

2.5  to  3 

C-8 

c-f  indicates  component-fluid  interface  and  c-s  indicates 
component-substrate  interface 
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Maximum  temperature  versus  Ra  is  shown  in  Fig.  9(a).  A  decrease  in  the 
maximum  temperature  from  0.122  to  0.042  is  observed  for  an  increase  in  Ra 
from  10s*  106.  Fig.  9(b)  shows  that  the  change  in  maximum  temperature 
versus  Ra  is  similar  for  all  w/dt.  Although  nondimensional  temperature 

decreases,  the  dimensional  temperature  may  go  up  because  the  non- 
dimensional  temperature  is  inversely  proportional  to  the  heat  transfer 
coefficient.  An  increase  in  Q  would  bring  about  an  increase  in  the 
dimensional  temperature  but  a  decrease  in  the  nondimensional 
temperature.  As  Q  increases,  the  natural  convection  flow  velocities 
increase  and  this  would  increase  the  heat  transfer  coefficient. 

R  EFFECT  OF  Pr 

Fluorinert  dielectric  liquids  used  in  immersion  cooling  have  Pr  values 
much  greater  than  1.  Computations  for  Pr=7,  Pr=25  and  Pr=100  were 
therefore  conducted.  These  revealed  that  with  other  parameters  remaining 
constant,  Pr  had  a  negligible  effect  on  temperature  contours,  streamlines, 
heat  flux  or  temperatures  as  is  seen  in  Figs.  10-17.  This  would  indicate  that 
the  results  of  experimental  studies  in  water  with  Pr=7  would  apply  to  the 
dielectric  liquids  as  well  for  the  same  value  of  Ra  and  other  parameters. 

C.  EFFECT  OF  SUBSTRATE  CONDUCnVITY  RATIO  Rg 

The  streamlines  and  temperature  contours  for  various  R,  values  when 
Ra=105  are  shown  in  Figs.  18(aMe).  For  .l^R^l  all  the  main  flow  is 

generated  by  the  heater  side  with  a  very  weak  cell  on  the  left  side  of  the 
enclosure.  When  R,  reaches  10  both  cells  are  of  about  the  same  strength 

and  the  streamlines  form  two  cells  that  start  on  opposite  sides  of  the 
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Figure  10(a).  Temperature  contours  and  streamlines  for  Pr*7 
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Figure  10(b).  Temperature  contours  and  streamlines  for  Pn=25 
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Figure  11.  Heat  loss  through  component  faces  for  various  Pr 
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Figure  12.  Heat  loss  through  enclosure  walls  for  various  Pr 


30 


« 


Figure  18(b).  Temperature  contours  and  streamlines  for  R,*l 


Figure  18(d).  Temperature  contours  and  streamlines  for  R^>50 
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substrate.  Temperature  contours  and  streamlines  for  10  100 

essentially  look  the  same  and  it  appears  that  there  is  little  change  in  the 
induced  flow  beyond  Rs=10.  Temperature  contours  show  a  very  high 
gradient  just  outside  the  component  for  R,=0.1  and  the  hot  spot  of  the 
component  is  in  the  top  left  corner.  For  R,=l  the  hot  spot  moves  towards  the 
center  of  the  component  and  as  R,  increases  the  component  hot  spot  moves 
towards  the  fluid. 

Figure  19  shows  the  percent  heat  loss  through  the  component  faces  for 
various  R,  with  Ra=105.  There  is  a  tremendous  decrease  in  the  percent 
heat  loss  to  the  fluid  from  the  front  face  when  Ra  is  increased  from  1  to  10 
after  which  it  shows  only  a  slight  decrease.  As  the  conductivity  of  the 
substrate  increases  it  is  easier  for  the  heat  to  flow  through  faces  other  than 
the  front.  The  substrate  acts  like  a  fin,  increasing  the  surface  area  where 
heat  loss  can  occur.  The  back  face  percent  heat  loss  increases  throughout 
the  range  of  R,  increase.  The  bottom  and  top  face  percent  heat  loss  increase 

when  Rs  is  increased  from  0.1  to  10  and  then  remain  fairly  constant. 

The  percent  heat  loss  through  the  enclosure  walls  for  various  Rg  is 

shown  in  Fig.  20.  The  bottom  wall  initially  shows  a  slight  decrease  when 
Rs  increases  from  0.1  to  1  due  the  left  cell  of  the  flow  in  the  fluid  bringing 

higher  temperatures  closer  to  the  bottom  wall  and  then  a  slight  increase 
over  the  remaining  range  of  R(  due  to  the  increased  heat  conduction 

through  the  bottom  portion  of  the  substrate  into  the  vicinity  of  the  bottom 
wall.  The  other  percent  heat  loss  changes  are  a  combination  of  the 
increased  thermal  conductivity  and  flow  changes. 
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PERCENT  HEAT  LOSS 


Figure  21  shows  the  substrate  surface  heat  flux  for  various  R,.  For 
Rg£  10  the  heat  flux  for  the  substrate  is  fairly  uniform.  For  Rs<10  most  of 
the  heat  flux  comes  from  the  front  face  of  the  component  for  S2  ranging 
from  2-3.  The  greatest  heat  flux  is  at  the  base  of  the  component  at  S2=2. 
The  same  pattern  follows  for  the  substrate  surface  temperatures  in  Fig.  22 
with  the  temperatures  being  uniform  for  R,  ^10  and  for  Rs<10  the 

temperatures  being  higher  over  the  component  surface,  S2=2-3,  with  the 
highest  temperature  at  the  top  of  the  component,  at  S2=3. 

Component  surface  heat  flux  is  shown  in  Fig.  23.  For  Rs<10  the  heat 

flux  through  the  component-fluid  interface  is  greatest,  for  Sl=0  to  1.  For 
R,^10  the  component-substrate  interface  shows  the  higher  heat  flux  with 

the  bottom  face  having  the  maximum.  Very  little  difference  occurs  when 
R,  is  increased  beyond  50.  The  component  surface  temperature  is  fairly 
uniform  for  all  Rs  due  to  the  high  conductivity  of  the  component  as  shown  in 
Fig.  24.  The  surface  temperature  decreases  greatly  for  R,=0.1  to  10  and 

then  decreases  slightly  over  the  remaining  range  with  very  little  change 
between  Rg=50  to  100.  As  the  substrate  conductivity  increases  there  is  less 

resistance  to  the  flow  of  heat  and  hence  the  surface  temperatures  decrease. 

Maximum  temperature  versus  Ra  is  shown  in  Fig.  25  and  shows  there  is 
little  benefit  to  increasing  R,  beyond  10. 

D.  EFFECT  OF  COMPONENT  CONDUCTIVITY  RATIO  R^. 

The  effects  of  Rc  on  the  temperature  contours  and  streamlines  are 
shown  in  Fig.  26  (a)-(d).  For  Rc=l  there  is  a  steep  thermal  gradient  in  the 
component  with  the  maximum  temperature  near  the  center  of  the 
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Figure  22.  Substrate  surface  temperature  for  various  R*,  Ra»10® 
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Figure  26(d).  Temperature  contours  and  streamlines  for  1^*100 
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component.  For  Rc=25  the  hot  spot  is  towards  the  fluid  face  and  the 
component  is  almost  at  a  uniform  temperature.  For  Rc=50  and  100  the 

temperature  contours  look  identical.  There  is  a  negligible  change  in  the 
streamlines  for  the  complete  range  of  Rc  studied. 

Figure  27  shows  the  percent  heat  loss  through  the  component  faces.  For 
Rc  increasing  from  1  to  25  there  is  a  large  decrease  in  the  percent  heat  loss 

from  the  front  face  and  a  slight  decrease  from  the  back.  A  large  increase  in 
the  percent  heat  loss  from  the  bottom  and  a  slight  increase  in  the  top  face 
make  up  for  the  decrease  at  the  front  and  back.  With  the  increase  in 
component  conductivity,  heat  flows  easier  through  the  component  and  a 
greater  percent  can  then  go  out  the  top  and  bottom.  The  substrate  has  a 
relatively  high  conductivity  and  the  top  and  bottom  subst  ite  areas  act  like 
fins.  For  R,.  increasing  from  25  to  100  there  is  a  slight  increase  in  the 

bottom  percent  heat  loss  and  a  very  slight  decrease  from  the  front  and  the 
back.  The  bottom  face  uses  the  substrate  below  as  a  fin  to  move  the  heat  to 
the  lower  cooler  region  of  the  enclosure.  The  heat  loss  through  the 
enclosure  walls,  shown  in  Fig.  28,  changes  insignificantly  which  is 
expected  due  to  no  change  in  the  fluid  flow  or  the  temperature  distribution 
in  the  fluid  over  the  entire  range  of  Rc  studied. 

The  component  surface  heat  flux,  Fig.  29,  shows  for  Rcsl  a  similar 
pattern  for  each  face.  Each  face  has  a  local  maximum  flux  at  the  center 
with  the  largest  maximum  flux  at  the  center  of  the  back  face  at  Sl=2.  The 
low  conductivity  of  the  component  limits  the  movement  of  heat,  so  the 
greater  flux  is  associated  with  the  surfaces  on  the  component  that  have  the 
most  volume  of  component  around  them.  These  locations  on  the  component 
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Figure  27.  Heat  loss  through  component  faces  for  various  Ra=105 
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Figure  28.  Heat  loss  through  enclosure  walls  for  various  Re,  Ra*l(P 
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surfaces  correspond  to  the  center  of  each  face.  For  25£RC£100,  the 
component  surface  heat  fluxes  are  nearly  identical  with  the  greatest  flux 
going  through  the  bottom  face  near  the  fluid-component  boundary.  The 
flux  is  greatest  here  since  the  substrate  extending  below  the  component 
bottom  surface  acts  like  a  fin.  In  addition  there  is  induced  flow  over  the 
substrate  below  the  bottom  of  the  component  augments  heat  transfer  from 
the  front  substrate  surface  below  the  heater. 

The  component  surface  temperatures  for  Rc=l  vary  widely  with  the 
maximum  occuring  at  SI  -0.7  as  seen  in  Fig.  30.  For  25£RC£100  the 

temperature  over  the  surface  is  fairly  uniform.  The  maximum 
temperature  versus  Rc  in  Fig.  31  shows  little  advantage  in  increasing  R* 

higher  than  25. 

E.  EFFECT  OF  COMPONENT  TO  SUBSTRATE  WIDTH  RATIO  w/dt 

Streamlines  and  the  temperature  contours  in  the  liquid  region  in  Figs. 
32(aMd)  show  only  minor  changes  when  w/dt  is  changed  from  0.25  to  0.999. 
This  occurs  because  the  high  Rs  and  Rc  result  in  low  solid  thermal 
resistance,  with  the  major  temperature  drop  occuring  in  the  fluid.  In  the 
component  region  the  hot  spot  moves  from  the  fluid  surface  to  the 
component  center  as  w/dt  increases  from  0.25  to  0.999. 

The  percent  heat  loss  through  the  component  faces  is  shown  in  Fig.  33. 
For  w/dt  increasing  from  0.25  to  0.999  there  is  approximately  a  10% 

increase  of  the  heat  loss  from  the  top  face  and  a  20%  increase  from  the 
bottom  face.  Both  the  top  and  the  bottom  percent  heat  losses  rise  because  of 
the  increased  heat  transfer  area  as  w/dt  increases  and  a  greater  volume  of 
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TEMPERATURE 


Temperature  Contours  Streamlines 


Figure  32(a).  Temperature  contours  and  streamlines  for  w/d^.25 
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Figure  32(b).  Temperature  contours  and  streamline*  for  w/d,aO^ 
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* 

Figure  32(d).  Temperature  contours  and  streamlines  for 
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the  component  next  to  them.  There  is  approximately  a  7%  decrease  in  heat 
loss  from  the  front  face  because  more  heat  is  being  taken  away  from  the  top 
and  bottom  and  the  front  is  actually  becoming  less  of  the  total  surface  area 
as  w/dt  is  increased.  There  is  approximately  a  23%  decrease  on  the  back 

face.  The  back  face  now  is  also  less  of  the  total  percent  area  of  the 
component  surface  and  the  back  has  less  of  the  substrate  acting  as  a  fin  as 
w/dt  increases  until  there  is  essentially  no  substrate  at  w/dt  =  0.999.  There 

is  no  difference  in  the  percent  heat  loss  from  the  enclosure  boundaries,  Fig. 
34,  as  was  suggested  from  the  temperature  contours  and  streamlines. 
Figure  35  shows  the  maximum  temperature  versus  w/dt.  There  is 

almost  a  linear  decrease  in  maximum  temperature  with  an  increase  in 
w/dt.  The  same  rate  of  heat  is  being  dissipated  by  a  larger  volume  which 

reduces  the  maximum  temperature. 
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PERCEENT  HEAT  LOSS 
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Figure  33.  Heat  loss  through  component  faces  for  various  w/d^  Ra=105 
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Figure  34.  Heat  loss  through  enclosure  walls  for  various  w/d,,  Ra^lO8 
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IV.  EXPERIMENTS 


A.  EXP] 


wriMiw 


NTAL  APPARATUS 


1.  Additions 

The  details  of  the  experimental  apparatus  can  be  found  in  Gaiser 
[Ref.  6].  It  consists  of  a  test  surface  with  3  columns  of  flush  heaters  with  15 
heaters  per  column  as  shown  in  Fig.  46.  The  test  surface  is  immersed  in  a 
large  plate  glass  tank  filled  with  deionized  water.  Each  heater  has  a 
thermocouple  attached  to  the  midpoint  of  the  back.  Voltages  of  the 
thermocouple  and  heater  leads  are  taken  by  a  data  aquisition  system  and  a 
HPCHewlett  Packard)  computer  determines  the  temperatures  and  power  of 
the  heaters.  Further  additions  and  changes  done  in  the  experimental  setup 
are  described  below. 

A  3-dimensional  traverse  assembly  with  controller  for  the  traverse 
motors  and  a  deaei  ting  system  for  the  water  were  installed.  The  traverse 
controller  interfaces  with  the  HP  computer  in  the  two  horizontal  directions 
and  a  program  was  developed  to  move  the  traverse  with  the  HP  computer 
integrating  the  temperature  acquisitions  with  a  probe.  The  traverse  must 
be  manually  moved  in  the  vertical  direction.  See  Appendix  C  for  the 
traverse  computer  program.  Deaeration  was  achieved  by  spraying  water 
through  a  nozzle  into  a  partially  evacuated  bottle.  The  water  was  then 
allowed  to  accumulate  under  the  vacuum  and  pumped  back  into  the  tank  at 
a  rate  similar  to  which  it  was  being  drawn  out.  Teflon  sheets  were  placed 
over  the  tank's  free  surface  to  reduce  the  aeration  of  the  water  after  it  had 
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gone  through  this  process.  When  the  deaeration  system  was  activated  and 
a  slight  (3"  or  greater  Hg)  vacuum  was  drawn,  the  bottle  would  collapse 
leaving  the  system  unoperable.  A  nonflexible  container  would  be  desirable 
for  a  better  operation  of  the  deaeration  system. 

2.  Changes 

A  change  was  made  in  the  way  the  power  input  to  the  heaters  was 
calculated.  Upon  taking  over  the  apparatus  it  was  noticed  that  the 
resistances  of  the  various  heaters  varied  by  as  much  as  8%  from  the  lowest 
resistance  reading  to  the  highest.  The  high  density  heater  board  was  taken 
apart  and  the  heaters  with  the  resistances  farthest  from  the  norm  were 
replaced.  The  new  heaters  were  of  a  different  batch  than  the  original 
heaters  and  they  all  had  a  resistance  of  10.4  ohms  originally  but  when  they 
were  measured  through  the  leads  after  they  were  installed  the  resistances 
again  varied  by  as  much  as  6%,  with  the  lowest  value  being  the  same  as  the 
heater  before  installation  and  the  highest  resistance  being  6%  above  the 
original  heater  resistance.  Since  all  the  heaters  were  the  same  before  being 
installed  it  was  assumed  that  the  excess  resistance  came  from  the  leads. 
The  excess  resistance  for  all  leads  was  then  determined  by  subtracting  the 
original  heater  resistance  from  the  resistance  of  the  heater  and  the  leads 
and  a  computer  program  that  took  this  into  account  was  implemented.  1% 
and  5%  precision  resistors  in  parrallel  with  the  original  1%  precision 
resistor  were  used  to  increase  the  current  of  the  heaters  that  had  less  power 
due  to  the  excess  lead  resistance.  The  following  relation  was  used  to 
determine  the  correct  power: 
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p.fV‘C.<V-.v‘4v-»J 

v  RP  J  (9) 

where  the  symbols  are  identified  in  the  Nomenclature.  Each  heater  and 
lead  resistance  had  to  be  treated  individually  and  in  some  cases  extra  line 
resistance  had  to  be  added  to  the  leads  to  make  the  heaters  have  all  the 
same  power.  A  uniform  and  nonuniform  power  program  were  developed 
for  the  HP  computer.  The  uniform  power  program  was  used  when  each 
heater  was  given  the  same  power.  The  nonuniform  power  program  was 
used  when  the  power  to  the  side  columns  was  different  than  the  center 
column  of  heaters.  Also  nonuniform  and  uniform  data  reduction  progams 
were  developed  to  obtain  the  temperatures  and  to  reduce  and  store  the  data 
simultaneously  for  each  case.  See  Appendixes  D-H  for  computer  programs 
utilized  in  the  acquisition  of  experimental  data. 

R  LIQUID  CRYSTALS 

Hallcrest  thermochromic  liquid  crystals  whose  optical  properties  of 
the  planar  texture  give  rise  to  the  selective  reflection  of  bright  colors  that 
change  with  temperature  [Ref.  13]  were  used  in  the  experiment.  The 
liquid  crystals  had  to  be  kept  refrigerated  prior  to  usage.  The  first  step  in 
liquid  crystal  application  was  to  obtain  a  dean  dry  surface.  A  black  backing 
was  then  applied  with  an  air  brush.  After  the  black  backing  had  dried,  the 
liquid  crystals  were  mixed  with  a  binder  with  3  parts  binder  to  1  part  liquid 
crystal  slurry  and  the  mixture  was  applied  over  the  backing  using  the  air 
brush.  Care  was  taken  to  keep  the  layers  very  thin  to  reduce  any 
interference  with  the  heat  transfer  process.  The  water  in  the  tank  had  to  be 
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kept  very  clean  and  the  water  filtering  system  was  run  for  two  days 
constantly  before  placing  the  test  surface  with  the  liquid  crystals  on  it  in  the 
water.  If  any  contamination  was  present,  only  a  dark  blue  color  would 
show  up  which  was  difficult  to  see  on  the  black  backing.  Even  with  the 
water  being  very  pure,  the  liquid  crystals  deteriorated  and  could  only  be 
expected  to  last  less  than  two  weeks.  It  is  believed  that  the  crystals 
somehow  react  with  the  polar  water  molecule  and  become  inactive.  In 
dielectric  fluid  this  had  not  been  a  problem,  which  was  evident  in  past 
thesis  research  [Ref.  14]. 

The  crystals  were  calibrated  against  temperatures  using  a  platinum 
resistance  standard  (Rosemont  Engineering  Model  920A  Commuting 
Bridge)  and  a  Model  913A  Calibration  Bath.  Two  ranges  of  liquid  crystal 
sensitivity  to  temperature  were  used.  After  standardization  for  the  colors 
which  appeared  below  30  degrees  C,  the  crystals  compared  within  0.5 
degrees  C  with  the  thermocouple  temperatures  of  the  heaters,  but  at  higher 
temperatures  (greater  than  40)  it  was  difficult  to  determine  the  change 
from  blue  to  black  and  an  agreement  between  the  crystals  and  the  heater 
thermocouples  was  about  5  degrees  C  with  the  thermocouple  temperature 
always  higher  than  the  color  indicated. 

A  problem  with  working  with  the  crystals  was  to  photograph  the 
resulting  surface  temperature  color  patterns.  It  is  very  difficult  to  shoot 
into  a  glass  container  filled  with  water.  All  light  sources  outside  the  tank 
must  be  turned  off  or  all  that  will  be  seen  in  the  picture  is  the  refection  of 
the  photographer.  The  correct  light  source  was  also  an  important 
consideration.  Florescent  lights  do  not  have  a  complete  spectrum  of  light 
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being  emitted  and  the  pictures  taken  under  these  conditions  did  not  show 
the  colors  that  were  actually  present.  Taking  pictures  with  a  spot  lamp  on 
the  face  with  a  high  speed  film  (ASA  1600)  did  not  work  because  the  area 
with  the  spot  ended  up  being  burned  out  in  the  picture  and  the  surrounding 
area  of  the  picture  appeared  black.  The  best  method  found  was  to  obtain  a 
light  source  that  contained  the  wavelength  band  necessary.  It  was 
obtainable  at  a  local  photograph  store  and  it  was  used  with  a  film  speed  of 
400  ASA  or  higher.  Care  had  to  be  taken  to  keep  the  light  source  above  the 
water  and  shining  inside  the  tank.  Shielding  was  placed  so  the  bulb  would 
only  light  up  inside  the  tank  staying  away  from  the  refecting  glass  and 
lighting  only  the  test  surface.  The  light  was  only  switched  on  when  taking 
pictures  to  avoid  setting  up  fluid  variations  by  warming  up  the  tank  water. 

Figure  36  shows  the  view  that  is  shown  in  the  pictures  to  follow.  It  also 
shows  the  shroud  and  the  test  surface.  The  shroud  is  a  dear  plastic  that 
causes  a  constriction  to  the  flow  and  simulates  another  circuit  board  or  a 
wall  in  an  actual  situation  in  electronic  equipment.  At  no  power,  the  test 
surface  appears  black  and  corresponds  to  ambient  temperature  (18  degrees 
C).  The  first  picture,  Fig.  37  is  a  picture  with  no  shroud  and  a  power  level 
of  .5  watts  per  heater.  The  actual  colors  that  were  present  are  brown  which 
appears  at  24.0  degrees  C,  green  which  appears  at  25.3  degrees  C,  and  blue 
which  appears  at  28.3  degrees  C  and  disappears  at  44.0  degrees  C  where 
the  surface  returns  to  black.  The  black  and  white  pictures  can  only  be  used 
to  show  higher  temperature  regions  with  the  lighter  regions  on  the 
photographs  indicating  where  the  hot  spots  are.  The  color  patterns  become 
more  evident  near  the  middle  and  top  compared  to  the  bottom  of  the  heater 
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Figure  36.  Schematic  of  the  view  that  is  seen  in  the  following  pictures 
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columns.  The  temperature  patterns  seem  symmetric  for  all  three  columns 
and  there  appears  to  be  no  thermal  interaction  between  the  three  columns. 
Heat  being  brought  up  by  natural  convection  is  brought  up  the  face  and 
some  is  transfered  back  into  the  plexiglass  substrate.  As  the  power  is 
increased  with  no  shroud  the  color  patterns  become  more  evident  at  the 
bottom  of  the  heater  column.  At  three  watts  per  heater  power  input,  Fig.  38 
shows  again  that  there  is  no  interaction  between  the  columns.  The  color 
pattern  above  the  top  heater  indicate  higher  temperatures  in  that  region 
compared  to  the  0.5  watts  per  heater  case. 

A  shroud  is  placed  6mm  in  front  of  the  test  surface  in  Fig.  39  and 
entrainment  of  the  flow  is  evident  by  the  shape  of  the  temperature  patterns. 
The  shroud  prevents  entrainment  flow  from  the  front  and  water  must  come 
from  the  sides  towards  the  center  column.  This  actually  enhances  cooling 
of  the  outside  top  heaters  for  a  shroud  spacing  of  6  mm  as  seen  in  Figs.  39  & 
40,  compared  to  when  there  is  no  shroud.  The  bottom  heaters  appear  not  to 
be  effected  by  the  shroud  placement.  In  Fig.  41  the  shroud  spacing  is  3  mm 
and  the  power  input  per  heater  is  0.2  watts.  The  color  patterns  now  become 
visible  at  a  much  lower  power  level.  The  flow  at  the  top  is  greatly 
constrained  with  the  high  temperature  patterns  actually  merging  over  the 
center  column.  Temperatures  are  increased  drastically  at  the  3mm  shroud 
spacing  and  Fig.  42  shows  that  at  a  0.5  watt  power  input  per  heater  the 
temperatures  are  as  high  as  with  a  1  watt,  6  mm  shroud  spacing  case  or  a  2 
watt,  no  shroud  case. 

Nonuniform  power  levels  for  the  side  columns  and  center  column  at 
different  power  input  were  looked  at  for  no  shroud  and  3mm  shroud  cases, 
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Figure  38.  No  shroud  case,  3  watts  per  heater  power  input 


Figure  39.  6  mm  shroud  spacing,  JS  watt  per  heater  power  input 
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Figure  40.  6  mm  shroud  spacing,  1  watt  per  heater  power  input 


Figure  41.3  mm  shroud  spacing,  J2  watt  per  heater  power  input 
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Figs.  43-45  ,  and  showed  the  same  patterns  as  the  uniform  power  cases. 
For  the  no  shroud  case  there  did  not  appear  to  be  any  interaction  between 
the  three  columns.  For  the  small  shroud  spacing  cases  there  was 
considerable  side  entrainment .  It  was  difficult  to  go  to  higher  power  levels 
in  the  low  shroud  spacing  cases  because  air  bubbles  formed  over  the  top 
center  region  and  disrupted  the  flow  and  made  the  temperatures  increase 
greatly  even  for  power  levels  of  the  order  of  0.5  watts  per  heater. 


I 
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Figure  42.  3  mm  shroud  spacing,  iS  watt  per  heater  power  input 
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Figure  44.  No  shroud  case,  0.5  &  1  watt  per  heater  power  input 

(center/sides) 
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Figure  45.  6  mm  shroud  spacing,  J2&.5  watt  per  heater  input  (center/side) 
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V.  NUMERICAL  COMPARISON  WITH  EXPERIMENT 


A  WHAT  TO  COMPARE  ? 

The  liquid  crystal  temperature  patterns  show  that  there  are  pronounced 
3-dimensional  effects  in  the  flow  and  temperature  patterns  with  low  shroud 
spacings,  especially  near  the  top  of  the  3  columns  of  heaters.  With  no 
shroud  and  wider  shroud  spacings,  the  heat  transfer  and  fluid  flow 
problem  is  primarily  2-dimensiona?,  especially  at  the  lower  center  portion 
of  the  3  columns  and  at  the  center  of  the  heaters.  The  central  plane  of  the 
bottom  heater  of  the  center  column  (section  AA  in  Fig.  46)  was  chosen  for  a 
2-dimensional  numerical  comparison. 

R.  0.2  WATT  NUMERICAL  COMPARISON 

The  first  attempt  to  model  the  temperature  measured  with  a 
thermocouple  at  the  middle  of  the  back  of  the  heater  for  the  0.2  watt  case 
agreed  only  by  approximately  25%.  Figure  47  shows  the  heater  surface 
temperatures  versus  SI.  SI,  as  in  the  parametric  study,  is  the  counter¬ 
clockwise  contour  distance  around  the  heater  starting  at  the  lower  right 
corner  of  the  heater.  The  experimental  value  of  the  thermocouple  is  located 
at  SI  slightly  greater  than  1.5.  To  get  closer  to  the  experimental 
temperature  the  numerical  model  was  examined  again.  A  two  heater 
model  was  developed  to  account  for  the  effect  of  upper  heaters.  Figure  48 
shows  the  corresponding  temperature  contours  and  streamlines.  With  the 
two  heater  model  the  agreement  was  within  approximately  13%,  as  seen  in 
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Figure  46.  Schematic  diagram  of  the  experimental  heater  studied 
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Fig.  49.  The  heaters  were  placed  within  the  grooves  in  the  plexiglass  test- 
surface  using  a  low  thermal  conductivity  epoxy.  During  the  liquid  crystal 
work  it  was  observed  that  the  epoxy  thermal  conductivity  was  lower  than 
that  of  plexiglass,  so,  a  thermal  conductivity  of  1/2  that  of  plexiglass  was 
used  in  the  numerical  model  to  see  what  effect  it  would  have  on  the 
comparison.  It  proved  to  be  a  slight  improvement  over  the  plexiglass 
thermal  conductivity  case  but  not  as  much  of  an  improvement  as  the  two 
heater  model,  as  seen  in  Fig.  50.  The  dimensions  of  the  heater  were  then 
looked  at  closely.  The  heater  is  made  up  of  a  Kapton  sheet  and  an  inconel 
sheet  together.  The  dimensions  of  the  Kapton  were  the  dimensions  of  the 
heater  that  were  being  used  but  the  inconel  was  actually  the  resistor  that 
was  causing  the  increase  in  temperature  and  its  dimensions  were  smaller 
than  the  Kapton.  The  height  of  the  region  which  contained  the  inconel  was 
0.0062  m  and  the  width  was  0.022  m.  When  these  new  smaller  dimensions 
were  used  in  the  numerical  model  the  agreement  was  much  better  as  seen 
in  Fig.  51.  Temperatures  for  both  thermal  conductivities  (0.14  W/m-K  and 
0.07  W/m-K)  are  plotted  in  Fig.  51  and  the  agreement  with  the  experiment  is 
within  approximately  9%  for  the  one  heater  model. 

The  heater  was  a  layer  of  Kapton  over  a  sheet  of  inconel  as  stated  above. 
When  the  heater  was  modeled  as  two  layers,  the  numerical  prediction  was 
within  7  %  of  the  experimental  measurement  as  shown  in  Fig.  52.  When 
the  two  heater  model  was  used  with  the  smaller  dimensions  for  the  heater 
and  the  heater  divided  into  two  layers,  as  mentioned  above,  Fig  53  shows 
that  the  temperature  variation  for  the  two  different  substrate  conductivities 
bracketed  the  experimental  result,  as  was  expected. 
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Figure  49.  First  attempt  at  2  heater  model 


Figure  52.  One  heater  model,  heater  fptit  into  a  Kapton  and  Inconel  layer 


C.  HIGHER  POWER  NUMERICAL  COMPARISONS 

The  measured  component  temperatures  at  higher  power  levels 
experimental  values  did  not  agree  as  closely  as  the  0.2  watt  power  case. 
Figure  54  shows  that  as  the  power  input  per  heater  is  increased  the 
disagreement  between  the  numerical  and  experimental  values  increases. 
At  0.5  watt  power  input  the  percent  difference  is  approximately  8  %.  One 
possible  reason  for  this  is  that  as  the  power  is  increased,  the  downstream 
heaters  have  a  greater  effect  on  the  upstream  heater  temperatures  due  to 
greater  conduction  through  the  substrate.  Thus  the  experimentally 
measured  temperatures  are  higher  by  a  greater  extent  compared  to  the 
numerical  model.  Thus  for  comparison  with  higher  power  cases, 
additional  heaters  must  be  added  both  on  the  side  and  on  the  top  in  the 
model.  Another  possible  reason  for  the  percent  difference  increase  with 
power  is  that  the  properties  for  the  water  used  in  the  model  were  taken  at 
the  average  temperature  of  the  heater  and  the  ambient  water.  Also,  the 
inconel  region  of  the  heater  was  taken  as  a  uniform  heat  source  in  the 
model  but  in  reality  is  a  ribbon  that  snakes  up  and  down  on  the  Kapton 
sheet.  At  higher  power  levels  this  could  cause  hotter  areas  near  the  ribbon 
and  this  could  cause  the  numerical  value  to  be  slightly  lower.  Still  another 
possibility  is  that  there  are  3-dimensional  effects  in  reality  that  the  2- 
dimensional  model  cannot  detect  and  these  effects  get  larger  with  increased 
power.  Entrainment  of  the  warmer  fluid  from  the  sides  columns  into  the 
center  region  of  the  heater  could  make  the  experimental  value  greater  than 
the  2-dimensional  model. 
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EXPERIMENTAL  VALUE 
SIBSTRATE  K..14  W/m-K 
SUBSTRATE  K..07  W/m-K 


D.  UPPER/LOWER  HEATER  NUMERICAL  EVALUATION 

The  numerical  model  for  two  heaters  provided  the  upper  and  lower 
heater  temperatures.  A  study  of  the  effects  of  the  lower  heater  on  the  upper 
heater  was  conducted.  Figures  55-59  show  the  progression  of  the  top  and 
bottom  heater  temperatures  as  the  power  is  increased  from  0.2  to  3  watts 
equally  in  both  heaters.  Figures  60-61  collect  the  temperature  patterns  at 
the  various  power  levels  studied  for  the  two  heaters.  Similar  patterns  are 
observed  as  power  level  increases  for  all  top  heaters  and  all  bottom  heaters 
when  compared  with  themselves.  When  comparing  the  lower  heater  to 
upper  heater  temperatures,  there  is  a  larger  increase  in  the  top  heater's  top 
face  temperatures.  The  back  face  temperatures  do  not  decrease  as  much  on 
the  top  heater  as  the  bottom  heater.  Also,  the  bottom  face  temperatures  of 
the  top  heater  do  not  decrease  below  the  temperature  of  the  bottom  of  the 
front  face  contrary  to  the  bottom  heater.  It  is  more  difficult  for  the  upper 
heater  to  lose  heat  because  of  the  natural  convection  bringing  up  the  heat 
from  the  lower  heater,  thus  the  upper  heater  temperatures  are  higher  than 
the  lower  heater  temperatures.  The  front  face  temperatures  of  both  heaters 
are  always  closer  than  the  other  faces.  As  power  level  increases,  the  top, 
back  and  bottom  faces  become  farther  from  their  counterparts  on  the  bottom 
heater. 

Figures  62-65  compare  the  heat  loss  through  the  heater  face  for 
increasing  power  levels.  The  only  face  that  shows  an  increase  for  the  top 
and  the  bottom  heater  is  the  front  face.  At  low  power  levels  the  percent  heat 
loss  through  the  front  face  with  increasing  power  is  its  largest.  Both 


82 


Figure  55.  0.2  watt  case  top  and  bottom  heater  surface  temperatures 


Figure  56. 0i5  watt  surface  temperatures  for  the  top  and  bottom  heater 
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Figure  57. 1  watt  per  beater  top  and  bottom  heater  surface  temperatures 


r 


BOTTOM  HEATER 
TOP  HEATER 


0 


1 


3 


4 


POWER  (WATTS) 


Figure  63. 0.2-3  watt  top  face  percent  heat  loss  for  top  and  bottom  heater 


Figure  64. 0.2  -3  watt  back  face  percent  heat  loss  for  top  and  bottom  heater 
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heaters  show  a  similar  pattern.  The  back  and  bottom  faces  percent  heat 
loss  both  decrease  for  increasing  power  levels 

The  effect  on  the  lower  heater  of  a  heater  above  it  is  shown  in  Figs.  66-67. 
For  0.2  watt  power  input  per  heater  shown  in  Fig.  66,  the  one  heater  model 
surface  temperatures  are  all  less  than  the  corresponding  2  heater  model 
surface  temperatures  for  the  bottom  heater  For  3  watts  power  per  heater, 
shown  in  Fig.  67  the  opposite  is  true.  The  surface  temperatures  for  the 
bottom  heater  of  the  2  heater  model  are  lower  except  at  the  top  of  the  heaters 
where  it  is  equal  to  the  one  heater  model.  The  increased  convection  for  the  2 
heater  model  when  the  heat  input  is  much  higher  causes  lower 
temperatures  at  the  bottom  of  the  lower  heater.  Thus  for  0.2  watts, 
conduction  from  the  upper  heater  to  the  lower  heater  causes  an  increase  in 
temperature  of  the  bottom  heater,  whereas  for  a  higher  power  input  of  3 
watts,  the  increased  convection  caused  a  cooling  effect  on  the  bottom  heater. 


Figure  65. 0.2-3  watt  bottom  face  percent  heat  loss  for  top  and  bottom  heater 


Figure  66. 0.2  watt  surface  temperatures  for  bottom  heater  for  1&2  heater 

models 
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1  HEATER  MODEL 

2  HEATER  MODEL 


VL  CONCLUSIONS 


During  the  course  of  the  numerical  and  experimental  study,  several 
important  conclusions  were  drawn  about  the  transfer  of  heat  from  a  flush 
heat  source.  A  constant  flux  of  temperature  boundary  condition  must  be 
used  carefully  as  the  numerical  study  showed  that  it  can  be  far  from 
uniform  when  using  the  conjugate  solution  to  solve  the  problem.  With  the 
advent  of  faster  computers  it  seems  that  all  heat  transfer  problems  should 
be  solved  with  a  combined  conduction/convection  procedure  to  improve 
accuracy  of  the  solution. 

Changes  in  the  Pr  between  7  and  100  did  not  change  the  nondimensional 
temperatures.  This  was  extremely  important  to  this  study  because  the 
experimental  work  has  been  done  in  water,  with  a  Pr  of  7  and  most  fluids 
that  are  currently  being  used  for  the  immersion  of  electronic  components 
are  fluorinerts  and  have  a  Pr  greater  than  10  and  less  than  100.  The 
experimental  results  should  therefore  be  applicable  to  these  fluids  as  well. 

Increasing  the  R,  provides  substantial  cooling  enhancement  when  Rs  is 

less  than  10  but  increasing  it  past  10  only  provides  slight  enhancement.  For 
the  designer  of  electronic  equipment,  this  should  help  in  the  decision  as  to 
how  high  a  conductivity  the  substrate  should  be.  When  a  choice  of  more 
than  one  material  is  available,  the  cheapest  material  with  an  R,  of  10  or 

greater  could  be  chosen. 

An  increase  in  Rc  provides  lower  and  more  uniform  component 
temperatures  primarily  for  Rc  less  than  25.  The  designer  can  once  again 
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look  for  the  cheapest  material  that  would  do  the  job  with  an  Rc  of  25  or 
greater  and  optimize  the  cooling  of  the  component  as  well. 

There  is  almost  a  linear  decrease  in  maximum  temperature  with 
increase  in  w/dt.  This  would  suggest  that  the  componet  should  be  made  as 

wide  as  the  substrate  if  the  design  permits. 

Finally,  liquid  crystals  proved  to  be  a  valuable  tool  in  determining  the 
areas  of  the  hot  flow  up  the  front  of  an  experimental  surface.  The  effect  of 
the  shroud  spacing  became  very  evident  as  the  spacing  decreased. 
Although  at  small  shroud  spacing  the  3-D  effects  of  the  fluid  flow  were 
evident,  at  large  spacing  and  no  shroud  conditions  and  near  the  bottom  of 
the  high  density  board  always  the  heat  transfer  process  could  be  modeled  as 
2-dimensional.  Very  good  agreement  of  numerical  and  experimental 
results  were  seen,  especially  at  the  low  power  levels. 
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VII.  RECOMMENDATIONS 


In  continuation  of  this  study,  it  is  suggested  that  the  same  configuration 
of  three  high  density  heater  columns  be  used,  however  it  is  recommend  that 
the  following  aspects  be  further  investigated: 

•Utilize  liquid  crystals  to  determine  the  spacing  where  shroud  effects  are 
negligible 

-Use  a  3  heater  model  to  determine  what  effects  this  would  have  on  the 
lower  heater  at  the  higher  power  levels 

-Analyse  the  effects  of  using  non-steady  power  inputs  on  the  heat 
transfer 

-Utilize  a  3-dimensional  numerical  model  to  see  if  3-D  effects  are 
significant,  especially  at  the  high  power  levels 

-A  nonflexible  container  should  be  acquired  for  the  dearation  system  to 
allow  the  use  of  the  test  surface  at  higher  power  levels 

-High  power  studies  to  determine  the  onset  of  turbulence  could  be  looked 
at  with  the  accompanying  effects  on  the  heat  transfer 
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APPENDIX  A 

USER  PORTION  OF  FINITE  DIFFERENCE  PROGRAMS 


C  ********** THESIS******* **************************** ************** 

c 

BLOCK  DATA 

LOGICAL  LSOLVE , LPRINT , LBLK , LSTOP 

COMMON  F  (  7  5 , 7  5 , 5 ) ,  P  (  75,75) , RHO (75,75), GAM (75,75) , CON (75,75) , 

1  AIP(75,75) ,AIM< 75,75 ) ,AJP( 75,75 ) ,AJM( 75,75 ) ,AP( 75,75 ) , 

2  X ( 7 5 ) , XU (75) , XDI F ( 7 5 ) , XCV (75) ,XCVS(75) , 

3  Y(75) ,YV(75) , YDIF(75) ,YCV(75) ,YCVS(75) , 

4  YCVR (75) , YCVRS (75) ,ARX(75) ,ARXJ{75) ,ARXJP(75) , 

5  R ( 7 5 ) , RMN( 7  5 ) , SX( 75 ) ,SXMN(75) ,XCVI(75) ,XCVIP(75) 

COMMON  DU (75, 75) ,  DV  (  7  5 , 7  5  ) , FV ( 7  5  ) ,FVP(75) , 

1  FX ( 7  5 ) , FXM (75) ,FY(75) ,FYM(75) , PT( 75 ) , QT( 75 ) 

COMMON/ I NDX/NF , NFMAX , NP , NRHO , NGAM , Ll , L2 , L3 , Ml , M2 , M3 , 

1IST, JST, ITER, LAST, TITLE ( 13) , RELAX (13) , TIME , DT, XL , YL , 

2IPREF , JPREF , LSOLVE ( 10 ) , LPRINT( 13 ) , LBLK( 11 ) , MODE , NTIMES ( 1 0 ) , 
3RHOCON , ZERO 
COMMON/CNTL/LSTOP 
COMMON/SORC/SMAX , SSUM 
COMMON/COEF/FLOW , DI FF , ACOF , ALFA 
DIMENSION  U(75,75) , V( 75 , 75 ) ,PC(75,75) 

EQUIVALENCE (F(1,1,1),U(1,1)) ,(F(1,1,2),V(1,1)),(F(1,1,3), 
1PC(1,1) ) 

DIMENSION  TH (75) , THU ( 7 5 ) , THDI F ( 7 5 ) ,THCU(75) ,THCVS(75) 
EQUIVALENCE ( X ( 1 ) , TH ( 1 ) ) , ( XU ( 1 ) ,THU(1 ) ) ,( XDI F ( 1 ), THDI F ( 1 ) ) , 

1 ( XCV ( 1 ) , THCV ) , ( XCVS ( 1 ) , THCVS ( 1 ) ) , ( XL , THL ) 

DATA  NFMAX,NP,NRHO,NGAM/5,6,7,8/ 

DATA  LSTOP, LSOLVE, LPRINT/1* . FALSE. , 10* . FALSE. , 13* . FALSE./ 

DATA  MODE, TIME, ITER/1, 0.  ,0/ 

DATA  RELAX, NTIMES/13*1 ., 10*1/ 

DATA  LBLK/1 1  * . TRUE . / 

DATA  DT, IPREF, JPREF, RHOCON/1 . E+ 1 0 , 1 , 1 , 1 . 0/ 
^******************************#*********************<1************ 

C  ENTER  DATA  INTO  LABELED  COMMONS  USING  DATA  STATEMENTS 

C****************************************************************** 

DATA  TITLE( 1 ) , TITLE ( 2 ) , TITLE ( 3 ) , TI TLE ( 4 ) , 

1  TITLE (6 )/'VELU' , 'VELV' , ' STFN ' , 'TEMP' , 

2 ' PRES'/ 

DATA  ( LSOLVE ( I ) , 1-1 , 4 ) , LSOLVE ( 6 ) , LPRINT ( 3 )/6* .TRUE ./ 

DATA  RELAX ( 1 ) , RELAX ( 2 ) , RELAX ( 4 ) , RELAX ( 6 ) , ALFA/. 04 ,.04,.4,.4,.l/ 
DATA  LAST/5000/ 

C  DATA  LAST/5/ 

END 

C************************************************************ 

SUBROUTINE  USER 

c************************************************************ 
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COMMON  F(75,75,5),P<75,75), RHO (75,75), GAM (75,75), CON (75,75), 

1  AIP(75,75) ,AIM(75,75) ,AJP(75,75) ,AJM(75,75) ,AP(75,75) , 

2  X( 75) ,XU( 75) ,XDIF(75) ,XCV( 75) ,XCVS( 75) , 

3  Y  (  7  5  ) ,YV(75) , YD I F ( 7 5 ) , YCV ( 7  5 ) ,YCVS(75)  , 

4  YCVR( 7  5 ) , YCVRS (75) ,ARX(75) ,ARXJ(75) ,ARXJP(75) , 

5  R( 75) ,RMN( 75) ,SX( 75) ,SXMN( 75) ,XCVI ( 7 5 ) , XCVI  P  (  7 5 ) 

COMMON  DU ( 7 5 , 7 5 ) , DV (75,75) , FV ( 7 5 ) ,FVP(75)  , 

1  FX ( 7  5 ) , FXM (75) , FY ( 7 5 ) , FYM ( 7 5 ) , PT ( 7 5 ) , QT ( 7 5 ) 

LOGICAL  LSOLVE , LPRINT, LBLK , LSTOP 

COMMON/I NDX/NF , NFMAX , NP , NRHO , NGAM ,Ll,L2,L3,Ml,M2,M3, 

1IST, JST, ITER, LAST, TITLE ( 13) , RELAX! 1 3 ) , TIME , DT , XL , YL , 
2IPREF,JPREF, LSOLVE (10) , LPR1NT( 1 3 ) , LBLK ( 1 1 ) ,  MODE , NTIMES ' 1 0 ) , 
3RHOCON , ZERO 
COMMON/CNTL/LSTOP 
COMMON/SORC/SMAX , SSUM 
COMMON/COEF/FLOW , DI FF , ACOF , ALFA 
DIMENSION  U(75,75) ,  V(  7 5 , 7 5  ) ,PC(75,75) 

EQUIVALENCE ( F ( 1 ,1,1),U(1,1)),(F(1,1,2),V(1,1)),(F(1,1,3), 
1PC(1,1) ) 

DIMENSION  TH (75) ,THU(75) ,THDIF(75) ,THCU(75) ,THCVS(75) 
EQUIVALENCE (X( 1 ) , TH ( 1 ) ) , ( XU ( 1 ) , THU ( 1 ) ) , (XDIF( 1 ) , THDI F ( 1 ) ) , 

1 ( XCV( 1 ) , THCV ) , ( XCVS ( 1 ) , THCVS ( 1 ) ) , ( XL , THL ) 

DIMENSION  T (  7  5 , 7  5 ) , STO ( 7 5 , 7 5  ) , I  BODY (  7 5 , 7 5 ) 

EQUIVALENCE  (T(1,1),F(1,1,4)) 

LOGICAL  LREAD, LWRITE 

DIMENSION  D(10),H(10),W(10),Q(10) 

DIMENSION  TTOP (75) , TSI DES ( 2 , 7 5 ) ,TBOT(75) ,QTOP(75) , QS I DES  (  2 , 7 5 )  , 
1QBOT (75) ,QHTR( 75) ,THBSID( 75 ) ,THTOP( 75 ) ,THBOT( 75) 

CHARACTER* 1  AREAD 

ENTRY  GRID 

c*********************************************** 

wr i te ( 36  ,  *  )  '  ks  «  10  ,ra=10e6  ' 

WRITE! 36  ,  *  ) 

C  SET  UP  THE  GRID 

NlX-11 
N2X-6 
N3X-9 
N4X-11 
N5X-17 
N6X-8 
N7X-10 
NlY-8 
N2Y-12 
N3Y-15 
N4Y-15 
N5Y-11 
N6Y-12 

C  HEATER  DIMENSIONS 
WIDTH- .02 
HI- . 0078 
DEPTH- . 5*H1 

C  SET  UP  THE  DIMENSIONS  OF  THE  ENCLOSURE 

ALl-3*Hl/Hl 
AL2-H1/H1 
AL3«( Hl-DEPTH )/Hl 
AL4-DEPTH/H1 
AL5- 2  *  Hi /H 1 
AL6-2*Hl/Hl 
AL7-H1/H1 


95 


BL1-2*H1/H1 

BL2-2*Hl/Hl 

BL3-H1/H1 

BL4-2*H1/H1 

BL5-2*Hl/Hl 

BL6-H1/H1 

DX1-AL1/N1X 

DX2-AL2/N2X 

DX3-AL3/N3X 

DX4-AL4/N4X 

DX5-AL5/N5X 

DX6-AL6/N6X 

DX7-AL7/N7X 

DYl-BLl/NlY 

DY2-BL2/N2Y 

DY3-BL3/N3Y 

DY4-BL4/N4Y 

DY5-BL5/N5Y 

DY6-BL6/N6Y 

C  MAKING  OF  THE  GRID 

XU ( 2 ) -0 
FIRST-3 
LLAST-3-1+N1X 
DO  11  I-FIRST, LLAST 
XU ( I ) -XU ( 1-1 )+DXl 

11  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N2X 

DO  12  I-FIRST, LLAST 
XU ( I ) -XU ( I -1 } +  DX2 

12  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N3X 

DO  13  I-FIRST. LLAST 
XU( I )-XU( 1-1 )+DX3 

13  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N4X 

DO  14  I-FIRST, LLAST 
XU ( I )-XU( 1-1 ) +DX4 

14  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N5X 

DO  15  I-FIRST, LLAST 
XU( I ) -XU { 1-1 )+DX5 

15  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N6X 

DO  16  I-FIRST, LLAST 
XU ( I )-XU( 1-1 )+DX6 

16  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N7X 

DO  17  I-FIRST, LLAST 
XU< I ) -XU { 1-1 ) +DX7 

17  CONTINUE 
Ll -LLAST 
XL-XU( LLAST) 

YV ( 2  )  -  0 
FIRST-3 
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n  n 


LLAST-FIRST-1+N1Y 
DO  19  I-FIRST, LLAST 
YV ( I )-YV( 1-1 )+DYl 

19  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N2Y 

DO  20  I-FIRST, LLAST 
YV ( I ) « YV ( I -1 ) +  DY  2 

20  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N3 Y 

DO  21  I-FIRST, LLAST 
YV ( I ) -YV ( I -1 ) +DY3 

21  CONTINUE 

FIRST-LLAST+1 
LLAST-LLAST+N4 Y 
DO  22  I-FIRST, LLAST 
YV ( I ) * YV ( I - 1 ) +  DY4 

22  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N5 Y 

DO  23  I-FIRST, LLAST 
YV ( I ) - YV ( 1-1 ) +DY5 

23  CONTINUE 
FIRST-LLAST+1 
LLAST-LLAST+N6 Y 

DO  24  I-FIRST, LLAST 
YV ( I ) - YV ( I -1 ) +DY6 

24  CONTINUE 
Ml-LLAST 
YL- YV ( LLAST) 

RETURN 

C 

ENTRY  START 

C  CIRCUIT  BOARD  CORNERS 
NA-N1X+N2X+2 
NB-N1X+N2X+N3X+N4X+1 
NC-N1Y+2 

ND-N1Y+N2Y+N3Y+N4 Y+l 
C  TEMPERATURE  PRINTOUT  LOCATION 
ITL-INT( N3Y/2 ) +1 
C  HEATER  LOCATIONS 
NQX-NB 
NQY1-NC+N2Y 

C  INITIAL  CONDITIONS 

DO  100  J« 1 , M 1 
DO  100  1-1 , LI 
U(I,J)«0. 

V(I,J)-0. 

T( I , J  )  -0 . 

STO ( I , J)«0. 

100  CONTINUE 

** ‘RELAXATION  FOR  BOUYANT  FORCE  (NOW  FOUND  AT  THE  TOP) 
ALFA- . 1 
ALFAM--1  .  -ALFA 
C  PROPERTY  ENTRY  POINT 

RA= 1 . E6 
PR-25 
CP-4165 
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on  no  noon  no 


RHOF-999 . 0 

BETA-162 . 020E-6 

G-9.807 

CONFLD- . 596 

CON P LX- 1 0  * CONFLD 

CONHTR-2 5* CONFLD 

DYNVI S- PR*  CONFLD/CP 

XKK-CONFLD 

ALPHA-CONFLD/( RHOF*CP ) 

Q1-RA*ALPHA* CONFLD* DYNVIS/( RHOF *G* BETA* Hi  * H 1 *H1 ) 

USTAR«SQRT( G*  BETA*Ql ‘Hl/CONFLD ) 

AMU- ( DYNV IS/(RH0F*H1 * US TAR ) ) 

DUMMY-CP*RHOF*Hl*USTAR 
C  DUMMY 1  USED  FOR  DIMENSIONAL  HEAT  FLUXES 
DUMMY 1-Ql * DUMMY/ ( HI  * CONFLD ) 

AMUP-CONFLD/DUMMY 
GAM 2- CONP LX/DUMMY 
GAM4-CONHTR/DUMMY 
XQ-01/(DEPTH*H1 ) 

EIN-Q1 

SOURCE- CONFLD/ (CP*DEPTH*  RHOF*USTAR ) 

PRINT  *,  '1  IF  WANT  TO  READ  FILE,  2  IF  STARTING  FROM  SCRATCH’ 

READ  *,  IFREAD 
I FREAD-2 

IF( IFREAD. EQ. 1 (THEN 
PRINT*,  'READING  DATA' 


*  *  *READ  DATA*** 

NTI ME-TOTAL  NUMBER  OF  TIME  STEPS 
I TIME- CURRENT  TIMESTEP 
rewind( 46  ) 

NTIME-1 

r ead ( 4  6 )  NTIME, Ml , LI 

PRINT  DATA  FOR  GRID 
read ( 46 )  ( YV( J ) , J-2 , Ml  ) 

read( 46 )  ( XU ( I ) , I -2 , LI ) 

read ( 46 )  ( Y( J ) , J-l , Ml ) 

read ( 46 )  ( X ( I ) , 1-1 , Ll ) 

READ  DATA  FOR  DEPENDENDT  VARIABLES 
I TIME- 1 

read ( 46  )  ITIME 
DO  3000  I - 1 , Ml 

3000  read ( 46 )  ( v( J , I ) , J-l , Ll  ) 

DO  3020  I - 1 , Ml 

3020  read( 46  )  (U( J , I ) , J-l , Ll ) 

DO  3030  I-l.Ml 

3030  read ( 46  )  ( T( J , I ) , J-l , Ll  ) 

DO  3040  I -1 , Ml 

3040  read ( 46  )  ( P( J, I ) , J-l ,L1  ) 

DO  3050  I - 1 , M 1 

3050  read ! 46  )  ( STO< J , I ) , J-l , Ll  ) 

DO  3060  1-1 , Ml 

3060  read ( 46  )  ( F ( J , I , 3 ) , J- 1 , Ll ) 

DO  3070  I-l.Ml 

3070  read ( 46  )  ( I  BODY ( J , I ) , J- 1 , Ll ) 
close ( 46  ) 

ELSE 
ENDI  F 
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RETURN 


ENTRY  DENSE 
RETURN 

ENTRY  BOUND 
DO  300  1-2, L2 
T ( I , Ml ) -0 
T( I , 1 ) -0 

300  CONTINUE 

DO  301  1-2 , M2 
T( LI , I ) -0 
T( 1 , I ) -0 

301  CONTINUE 
RETURN 

ENTRY  OUTPUT 

C  I F ( ITER . NE . 0 ) GO  TO  400 

C  PRINT  401 

C  401  FORMAT ( '  I TER ' , 6X , ' SMAX 1 , 8X , ' SSUM  7X V ( FRONT )' , 

C  1  6X, 'T( HEATER) ' ) 

C  400  PRINT  403, ITER, SMAX, SSUM, V(NQX+1, NQY1+2) , T( NQX , NQY1+2 ) 

C  403  FORMAT { I6,1P5E12.3) 

IF(ITER.NE.O) GO  TO  400 
wt ite (  36 , 401  ) 

401  FORMAT ( '  ITER ', 6X ,' SMAX ', 8X ,' SSUM  ', 7X ,’ V ( FRONT )’ , 

1  6X , 'T( THERMC) ' ) 

400  I F ( I NT ( ITER/25  ) * 2 5 . NE . ITER )GOTO  404 

write( 36 , 403 ) ITER, SMAX, SSUM, V(NQX+1 ,NQYl+ I TL ) , T( NQX , NQYl + I TL ) 

403  FORMAT ( I6.1P5E12.3) 

404  CONTINUE 

I F ( ITER. LT. LAST)  RETURN 
CALL  PRINT 
SSUM-0 . 

NCOUNT-O 

982  CONTINUE 

S2- ( YV ( 3 ) -YV ( 2 ) ) *Hl 

I F ( NCOUNT . EQ . 1 ) S2-( YV ( Ml )-YV(M2 ) ) *Hl 
LO  983  1-2 , Ll-1 

S1-(XU( 1+1 ) -XU { I ) ) *H1 

POWER-2*XKK*  SI  * ( T ( I , 2 ) *Ql/CONFLD ) /S2 

I F ( NCOUNT . EQ . 1 ) POWER-2*  XKK* SI  * ( T ( I , M2 ) +Q1/CONFLD ) /S2 
SSUM-POWER+SSUM 

983  CONTINUE 
NCOUNT-NCOUNT+ 1 
IF(NCOUNT.LT.1.5)  GO  TO  982 
NCOUNT-O 

984  CONTINUE 

S2« ( XU ( 3 ) -XU ( 2 ) ) *H1 

I F ( NCOUNT . EQ.l)S2«(XU(Ll )-XU( L2 ) )*Hl 
DO  985  I-2.M1-1 

Sl-( YV( 1+1 ) - YV ( I ) ) *H1 

POWER- 2* XKK*  SI  * ( T( 2 , I ) *Ql/CONFLD )/S2 

IF( NCOUNT. EO. 1 )  POWER- 2  * XKK* Si  * ( T ( L2  ,  I  ) *  01/CON FLD ) /S 2 
SSUM-POWER+SSUM 

985  CONTINUE 
NCOUNT-NCOUNT+1 

I F ( NCOUNT . LT .1.5)  GOTO  984 
wc i te ( 36  ,  *  ) 

wr  i te ( 36 , *  )  'ENERGY  IN  -  ',EIN 
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ENERGY  OUT 


SSUM 


wcite( 36 , * ) 
write( 36  ,  *  )  ' 
write( 36 , *  ) 
write (36,*)  'SURFACE  TEMPERATURES' 
write( 36  ,  *  ) 

DO  987  1-1,75 
TTOP ( I ) -0 
TBOT ( I )«0 
QTOP ( I ) -0 
QBOT( I )-0 
QHTR( I ) — 0 
THTOP( I )-0 
THBOT( I ) -0 
THBSIDt I )-0 
DO  987  J-1,75 

TSIDES ( 1 , 7  5 ) -0 
TSIDES ( 2 , 75 ) -0 
QS1DES (1,75)— 0 
QSIDES ( 2 , 7  5 ) -0 

987  CONTINUE 
C 

*  CALCULATION  OF  SUBSTRATE  SURFACE  TEMPERATURES 
DO  988  I-NA , NB 

TTOP ( I ) » ( T( I , ND ) * . 5* YCV (ND+1 )/CONFLD+T( I , ND+1 ) * . 5 * YCV ( ND ) /CONPLX ) / 
1(  . 5*  (  YCV ( ND+ 1 ) /CONFLD  +  YCV {  ND  ) /CONPLX ) ) 

TBOT ( I )«(T( I , NC ) * . 5* YCV { NC-1 )/CONFLD+T( I , NC-1 ) * . 5  * YCV ( NC )/CONPLX )/ 
1( . 5* (YCV( NC-1 )/CONFLD+YCV(NC) /CONPLX) ) 

988  CONTINUE 

DO  989  I-ND.NC.-l 

TSIDESd  ,  I  )-(T(NA,  I  )  *  .  5  *XCV{  NA-1  )/CONFLD+T( NA-1  ,  I  )  *  .  5*XCV(NA) 
1/CONPLX )  /  (  .  5  *  ( XCV ( NA- 1 ) /CONFLD+XCV ( NA ) /CONPLX ) ) 

TSIDES(2,I)-(T(NB,I )* .5*XCV(NB+l)/CONFLD+T(NB+l,I )*.5*XCV(NB) 
1/CONPLX )/( . 5* (XCVINB+l ) /CONFLD+XCV ( NB ) /CONPLX ) ) 

I F ( ( I . LT.NQY1 ) .OR. ( I . GT. NQY1+N3Y-1 ) )G0T0989 

TSIDES (2,I)-(T(NB,I ) * . 5*XCV(NB+1 )/CONFLD+T ( NB+1 , I ) * . 5* XCV ( NB ) 
1/CONHTR ) /( . 5* ( XCV ( NB+1 ) /CONFLD+XCV ( NB )/CONHTR ) ) 

989  CONTINUE 

wr i te ( 36 , *  )  'TOP' 

wr i te ( 36 , 4000 )  { TTOP ( I ) , I-NA , NB ) 

wri te { 36 , * )  'BACK  ','  FRONT' 

DO  990  I-ND , NC , — 1 

990  wr i te ( 36 , 4001 )  TSIDES ( 1 , I ) , TSIDES ( 2 , I ) 
wr ite ( 36 , * )  'BOTTOM' 

write( 36 , 4000 )  ( TBOT( I ) , I-NA ,NB ) 

4000  F0RMAT(1X,6E11.4  /  1X,6E11.4) 

4001  FORMAT ( 1X,E11.4,20X,E11.4) 

C 

C  SUBSTRATE  SURFACE  HEAT  FLUXES 
wr  i te( 36 ,  *  ) 

wr i te ( 36  ,  *  )  'SURFACE  HEAT  FLUXES-NONDI MENSIONAL ' 
write( 36  ,  *  ) 

DO  994  I-NA , NB 

QTOP ( I ) « ( T( I ,ND)-T( I , ND+1 )  )/(  . 5  * YCV( ND )/GAM2+ . 5* YCV ( ND+ 1 ) 

1/AMUP ) 

QBOT( I )-( T( I , NC ) -T( I , NC-1 ) )/(  . 5  * YCV ( NC ) /GAM2+ . 5  * YCV( NC- 1 ) 

1/AMUP) 

994  CONTINUE 

write(  36 , 4000  )  (QTOP( I ) , I-NA.NB) 

DO  995  I -ND , NC , -1 

QSIDES ( 1 , I ) - ( T ( NA , I )-T(NA-l , I ) )/( . 5* XCV ( NA ) /GAM 2+ . 5*XCV(NA-1 ) 
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n  n 


1/AMUP) 

QSIDES (  2 , 1 )-( T( NB , I )-T(NB*l , I) )/( . 5*XCV ( NB ) /GAM 2+ . 5*XCV(NB+1 ) 
1/AMUP) 

IF(  (  I .LT.NQY1 ) .OR. ( I . GT . NQY 1 +N3 Y- 1 ) 1G0T0995 

QSIDES(  2 , 1  )-(T(NB,  1  )-T(NB+l  ,  I)  )/(  .  5  *  XCV  (  NB  )  /GAM4  +  .  5  *XCV  (  NBt- 1  ) 
1/AMUP) 

995  CONTINUE 

DO  996  I -ND , NC , -1 

996  wr i te (  36 , 4001  )  QSIDES( 1 , 1 ) ,QSIDES( 2 , 1 ) 
write( 36,4000)  ( QBOTt I ) , I-NA , NB ) 

write( 36, * ) 

wr i te ( 36 , * )  'DIMENSIONAL  SUBSTRATE  HEAT  FLUXES’ 
wr i te ( 36 , * ) 

write( 36,4000)  (QTOP(I ) * DUMMY 1 , I-NA, NB ) 

DO  1050  I-ND,NC,-1 

1050  wr i te ( 36 , 4001 )  QSIDES ( 1 , I ) *DUMMYl , QSIDES ( 2 , I ) *DUMMY1 
write( 36, 4000 )  ( QBOT( I ) *DUMMYl , I-NA.NB) 


wr i te ( 36 , * ) 

write (36,*)  'LOCATION  OF  SUBSTRATE  TEMPERATURES  AND  FLUXES' 
wr i te ( 36  ,  *  ) 

write (36,*)  'TOP  ROW  ','1  Y','  &  X''S' 

write(36, ' (1X.E11.4  /  6E11.4  /  6E11.4)')  YV( ND+1 ) , ( X( I ) , I-NA, NB ) 
wr i te ( 36 , * ) 

write (36,*)  'BACK  AND  FRONT  X''S' 
write(36,'(lX,E11.4,2X,Ell.4)’)  XU(NA) ,XU(NB+1 ) 
wr i te  (  36  ,  * )  'Y  POSITIONS' 

DO  1020  I -ND , NC , - 1 

1020  write< 36, '  (5X.E11 .4 )  '  )  Y ( I ) 

wr i te ( 36 , * )  'BOTTOM  ROW  Y  AND  X''S' 

write(36, ' (IX, Ell. 4  /  6E11.4  /  6E11.4)')  YV ( NC ) , ( X ( I ) , I -NA , NB ) 
wri te( 36 , *  ) 

write (36,*)  'HEATER  SURFACE  TEMPERATURES' 

DO  1021  I-NA+N3X , NB 

THTOP ( I ) - ( T( I , NQY1+N3Y ) * . 5*YCV(NQYl )/CONHTR+T( I ,NQYl+N3Y-l ) * . 5* 

1 YCV ( NQYl +N3 Y ) /CONPLX )/( . 5* ( YCV ( NQY1 ) /CONHTR+ YCV ( NQY 1 +N 3 Y ) /CONPLX ) ) 
THBOT( I )-( T( I , NQYl ) * . 5*YCV(NQYl-l )/CONPLX+T( I , NQY1-1 ) * . 5* 

1 YCV (NQYl ) /CONHTR ) / ( . 5 * ( YCV( NQYl - 1 ) /CONPLX+ YCV ( NQYl )/CONHTR) ) 

1021  CONTINUE 

DO  1022  I-ND , NC , -1 

THBSID( I ) - ( T( NA+N3X , I ) * . 5* XCV ( NB-N4X )/CONPLX+T ( NB-N4X , I ) * . 5 
1 *XCV( NA+N3X ) /CONHTR )/( . 5* ( XCV ( NB-N4X )/CONPLX+XCV ( NA+N3X ) /CONHTR ) ) 

1022  CONTINUE 

wr i te( 36 , *  )  'HEATER  TOP' 

write (36, 4000)  (THTOP (I ) , I -NA+N3X , NB ) 

wri te ( 36 , *  )  'BACK  ','  FRONT' 

DO  1025  I-NQY1+N3Y-1 , NQYl , -1 
1025  wr ite(  36, 4001  )  THBSID( I ) , TSIDES ( 2 , 1  ) 
wri te( 36 , * )  'BOTTOM' 

write (36, 4000)  ( THBOT ( I ) , I -NA+N3X , NB ) 

HEATER  SUFACE  HEAT  FLUXES 
wr i te ( 36 ,  *  ) 

wr i te ( 36  ,  *  )  'HEATER  SURFACE  HEAT  FLUXES  -  NONDI MENS  I  ONAL ' 
wr i te ( 36 , *  ) 

IHTRC-1 

DO  997  I-NQX-N4X+1 ,NQX 

QHTR ( I HTRC ) - ( T ( I ,NQYl+N3Y-l )-T( I , NQYl +N  3 Y ) )/(  . 5* YCV ( NQY 1 +N3Y- 1 ) 
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1  /GAM4+ . 5*YCV(NQYl+N3Y  J/GAM2 ) 

IHTRC-IHTRC+1 

997  CONTINUE 

DO  998  I-NQX-N4X+1 ,NQX 

QHTR< IHTRC)-(T( I ,NQYl )-T( I ,NQYl-l ) )/( . 5*YCV(NQYl ) 

1  /GAM4+.5*YCV(NQY1-1)/GAM2) 

IHTRC-IHTRC+1 

998  CONTINUE 

DO  999  I-NQY1+N3Y-1 ,NQYl ,-l 

QHTR( I HTRC ) - ( T ( NQX-N4X+ 1 , I ) -T( NQX-N4X , I ) )/( . 5 *XCV ( NQX-N4X+ 1 ) 

1  /GAM4+ . 5*XCV ( NQX-N4X )/GAM2 ) 

IHTRC-IHTRC+1 

999  CONTINUE 

write( 36, 4000)  (QHTR( I ) , 1-1 , N4X ) 

HFTOT-O 

DO  1001  1-1 , N3Y 

wr i te( 36 , 4002  )  QHTR ( 2  *N4X+ 1 ) ,QSID£S( 2,NQY1+N3Y-I ) 

HFTOT-HFTOT+QSI DES ( 2 , NQY1+N3 Y- I ) * YCV( NQY1+N3Y-I ) 

1001  CONTINUE 

write(  36, 4000  )  (QHTR(N4X  +  I ) , 1-1 ,N4X) 

4002  FORMAT(1X,E11.4,10X,E11 .4) 
wr i te ( 36 , *  ) 

wli te ( 36 , *  )  'HEAT  GOING  DIRECTLY  TO  FLUID  -  ' , HFTOT  +  DUMHY 1  *  Hi 
wr i te ( 36 , *  ) 
write( 36  ,  *  ) 

wr i te { 36 , *  )  'DIMENSIONAL  HEATER  HEAT  FLUXES' 
write( 36,*) 

write(  36, 4  000  )  ( QHTR ( I ) * DUMMY 1 , 1-1 ,N4X) 

DO  1051  I- 1 , N  3 Y 

1051  wri te( 36, 4002)  QHTR( 2*N4X+I ) * DUMMY 1 , QSI DES ( 2 , NQY1+N3Y- I )*DUMMYl 
write( 36, 4000)  {QHTR(N4X+I ) *DUMMYl , 1-1 ,N4X) 
wr i te( 36 , * ) 

write( 36  ,  *  )  'LOCATION  OF  HEATER  TEMPERATURES  AND  FLUXES' 
wr i te ( 36  ,  *  ) 

wr  i  te  (  36  ,  *  )  'TOP  ROW  ','1  Y','  &  X"S' 

write( 36, ' ( IX , Ell . 4/6E1 1 . 4/6E1 1 . 4 ) ’ ) YV ( NQY1+N3Y ) , 

1 ( X ( I ) , I-NA+N3X,NB) 
wr i te ( 36  ,  *  ) 

wr i te ( 36  ,  *  )  'BACK  AND  FRONT  X''S' 

write( 36, ' ( IX, Ell .4 , 2X, Ell . 4 ) ' )  XU(NA+N3X) ,XU(NB+1 ) 
wr  i te { 36 , *  )  'Y  POSITIONS' 

DO  1030  I -NQYl +N3 Y-l , NQY1 , -1 
1030  wr i te ( 36, ' (5X.E11.4) ' )  Y<I) 

write (36,*)  'BOTTOM  ROW  Y  AND  X''S' 

write ( 36 , M 1X,E11 . 4/6E11 . 4/6  Ell , 4 ) ' )  YV ( NQYl ) , ( X ( I ) , I-NA+N3X , NB ) 
write( 3i, *  ) 

HPTOP-O 

HPBOT-O 

DO  1002  I-1.N4X 

HPTOP-HPTOP+QHTR ( I ) *XCV ( NQX-N4 X+ I ) 

HPBOT-HPBOT+QHTR ( N4X+I ) *XCV( NQX-N4X+ I ) 

1002  CONTINUE 
HPBACK-0 

DO  1003  I-2*N4X+1 , 2*N4X+N3Y 

HPBACK-HPBACK+QHTR ( I )*YCV(NQYl ) 

1003  CONTINUE 

HTRTOT-HFTOT+HPTOP+HPBOT+HPBACK 
PHTRF-IOO+HFTOT/HTRTOT 
PHTRT- 1 00  *  HPTOP/HTRTOT 
PHTRL-1 00*HPBACK/HTRTOT 


102 


n  n 


PHTRB-100*HPBOT/HTRTOT 

PHTOT-PHTRF+PHTRT+PHTRL+PHTRB 


MAXIMUM  TEMPERATURE 


wr i te  (  36  ,  * ) 

write  (36,*)  'HEATER  TEMPERATURES' 
wr  i  te ( 36,*) 

TMAX-0 

IXX-0 

jyy-o 

DO  1005  J-N3  Y ,1,-1 

write( 36,4000)  ( T( NQX-N4X+I .NQY1-1+ J ) , 1-1 ,N4X) 
DO  1005  K-1.N4X 

I F ( TMAX . LT . T( NQX-N4X+K , NQV 1-1+ J ) )  THEN 
TMAX«T{NQX-N4X+K,NQYl-l+J) 

IXX-K 

jyy-j 

ELSE 
ENDIF 
1005  CONTINUE 


1010 

1012 

1014 

1016 


wr i te ( 36  ,  *  ) 

wr i te( 36 , ' ( 1X,A, 3X.E11 .4 ) ' )  'MAX  TEMPERATURE  IS', TMAX 
wr i te ( 36 , *  ) 

wr i te ( 36 , *  )  'POSITION  FROM  LOWER  LEFT  IS' 
wr  i te ( 36 , *  )  'RIGHT', IXX 
wr  i te ( 36 , *  )  'UP’.JYY 
write( 36 , * ) 


wr ite ( 36 , *  ) 
wr i te ( 36 , *  ) 
write( 36  ,  *  ) 
write( 36 , *  ) 
wr i te ( 36  ,  *  ) 
wri te ( 36  ,  *  ) 
wri te( 36 , *  ) 
wr i te  <  36 , *  ) 
write( 36  ,  *  ) 
write( 36 , * ) 


'PERCENT  HEATER  ENERGY  LOST  TO  FLUID  -  ' , PHTRF 

'PERCENT  HEATER  ENERGY  LOST  TO  TOP  -  ' , PHTRT 

'PERCENT  HEATER  ENERGY  LOST  TO  BACK  -  ' , PHTRL 

'PERCENT  HEATER  ENERGY  LOST  TO  BOTTOM  -  ' , PHTRB 


'TOTAL  PERCENTAGE  -  ' . PHTOT 


'ENERGY  BALANCE' 


QLEFT-0 

DO  1010  1-2, M2 

QLEFT-QLEFT+T( 2, I )*Q1*YCV( I )/XDIFl 2) 

QUPPER-0 

DO  1012  1-2, L2 

QUPPER-QUPPER+T( I , M2 ) *Q1 *XCV ( I )/YDIF(Ml ) 

QRIGHT-0 

DO  1014  1-2, M2 

QRIGHT-QRIGHT+T( L2 , I ) *Q1*YCV( I )/XDIF(  LI ) 

QLOWER-O 

DO  1016  1-2, L2 

QLOWER-QLOWER+T ( 1,2) *Q1*XCV(I )/YDIF( 2) 

QOUT-QLEFT+QUPPER+QRIGHT+QLOWER 

PLEFT-QLEFT* 100/QOUT 

PUPPER-QUPPER* 1 00/QOUT 

PRI GHT-QRI GHT* 1 00/QOUT 

PLOWER-QLOWER* 1 OO/OOUT 

PTOT-PLEFT+  PUPPER+  PRI GHT+P LOWER 


write( 36  ,  *  ) 
wr i te ( 36  ,  *  ) 
wr i te ( 36  ,  *  ) 
wr i te ( 36 , *  ) 
wr i te ( 36  ,  *  ) 


'ENERGY  IN  EQUALS', EIN 
'ENERGY  OUT  EQUALS ’ , QOUT 
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on  on  noon 


wr i te ( 36  ,  *  ) 

'PERCENTAGE  OF  HEAT  LOSS  THRU  ENCLOSURE  WALLS' 

wr i te ( 36 ,  *  ) 

wr i te ( 36  ,  *  )  ' LEFT  WALL' , PLEFT 
wr i te ( 36 , *  )  'TOP  WALL' .PUPPER 
write(36,*)  ' RIGHT  WALL ', PRIGHT 
write (36,*)  'BOTTOM  WALL', FLOWER 

write( 36,  *  )  ' - - 

write (36,*)  'TOTAL  EQUALS ', PTOT 
wri te ( 36  ,  *  ) 

write ( 36 , * )  'PARAMETERS' 
write( 36, *  ) 

write( 36, * )  'H/Hl  -  ' , BL1+BL2+BL3+BL4+BL5+BL6 
write( 36, * )  'D/HI  -  ' ,BL2 
wr i te ( 36 , * )  'L/Hl  -  ',BL4 
wr i te ( 36 , * )  'ST/HI  -  ' , BL5+BL6 
wr i te ( 36 , * )  'SB/HI  -  ',BLl 
wr  i te ( 36  ,  *  )  ' DT/H1  -  ' ,AL3+AL4 
wr  i  te  (  36 , * ) 

write( 36,  *  )  ' RA  -  ' ,ra 
wri te( 36 , * )  'PR  -  ' , pr 
wr i te ( 36 , * )  'RS  -  ' , CONPLX/CONFLD 
wri te ( 36 , * )  'RC  -  ' , CONHTR/CONFLD 
wr i te ( 36 , * )  'W/Hl  -  ', DEPTH/HI 
wr i te ( 36, *  ) 

*  *  *  PREPARE  PLOTTING  DATA*** 

NTIME-TOTAL  NUMBER  OF  TIME  STEPS 
ITIME-CURRENT  TIMESTEP 
NTIME-1 

write( 46 )  NTIME , Ml , LI 

PRINT  DATA  FOR  GRID 
wr i te ( 46  )  ( YV( J ) , J-2 , Ml ) 
wr i te ( 46  )  ( XU ( I ) , 1-2 , LI } 
wr i te ( 4 6  )  ( Y ( J ) , J-1 , Ml  ) 
write( 46  )  ( X ( I ) , i-l , LI ) 

PRINT  DATA  FOR  DEPENDENDT  VARIABLES 
ITIME-1 

write( 46 )  ITIME 
DO  3001  I-l, Ml 
3001  wr i te ( 46 )  ( V ( J , I ) , J-1 , Ll ) 

DO  3021  I-l, Ml 

3021  wr i te ( 46 )  (U( J , I ) , J-1 , Ll ) 

DO  3031  I-l, Ml 

3031  wr i te ( 46 )  ( T( J , I ) , J-1 , Ll ) 

DO  3041  I-l, Ml 

3041  write(46)  ( P( J , 1  ) , J-1 , Ll ) 

DO  3051  I-l, Ml 

3051  wr i te ( 46 )  ( STO (J,I),J-1,L1) 

DO  3061  I-l, Ml 

3061  wr i te ( 46  )  { F( J , 1 , 3 ) , J-1 , Ll ) 

DO  3071  I-l, Ml 

3071  wr i te ( 46 )  ( IBODY< J , I ) , J-1 , Ll ) 
close ( 46 ) 

C 

RETURN 

C 

ENTRY  GAMSOR 
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IF(NF.EQ.3)  return 

DO  500  J-l, Ml 

DO  500  I-l.Ll 

GAM ( I , J ) -AMU 

I  BODY ( I , J ) -1 

I F ( NF . NE . 4 )  GO  TO  500 

GAM ( I , J ) -AMUP 

500  CONTINUE 

DO  501  I-NA, NB 
DO  501  J-NC , ND 
GAM( I , J ) -1 . 0E5 
IBODY ( I , J ) -0 
I F ( NF . NE . 4  )  GO  TO  501 
GAM (  I  ,  J  ) -GAM2 

501  CONTINUE 

I F ( NF . NE . 2 )  GO  TO  503 

DO  502  J-2,M2 

DO  502  1-2, L2 

IF( J . EQ . 2 )GO  TO  502 

TM-  (FY(J)*T(I,J) +FYM ( J ) *T( I , J-l  )  ) 

CON ( I , J )-ALFA*TM+ALFAM*STO( I , J ) 

STO( I , J ) -TM 

502  CONTINUE 

503  CONTINUE 

I F ( NF . EQ . 4 )  THEN 
DO  505  J-l , N4X 
DO  505  1-0 , N3Y-1 

GAM(NQX-N4X+J,NQY1+I )-GAM4 
CON(NQX-N4X+J,NQYl+I ) -SOURCE 
505  CONTINUE 
ELSE 
ENDIF 
RETURN 
END 
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on  oo  o  non  nnnnnnnnnnnnoon 


APPENDIX  B 

CONTOUR  PROGRAM 


MX*************************************************************** 
*********  PLOTTING  ROUTINE  FOR  A  SERIES  OF  CONTOUR  PLOTS  ****** 
***************************************************************** 

K  ....  CURRENT  NUMBER  OF  TIME  STEP  BEING  READ 

NT  -  TOTAL  NUMBER  OF  TIME-STEPS  ON  DATA  FILE 

NPTS  ...  TOTAL  NUMBER  OF  POINTS  (INPUT) 

NTYP  . . .  CURVE-TYPE 

1  . . .  PRESSURE 

2  .  .  .  STREAMFUNCTION 

3  ...  U-VELOCITY 

4  . . .  V-VELOCITY 

5  . . .  PRESSURE 


COMMON  /INDATA/  X( 1 02 ) , Y( 1 02 ) , ZC1 02, 1 02 ) 

COMMON  /GRIDII/  XX(  200  ) , YY ( 200  ) 

COMMON  /PL0TDT/  XPAGE , YPAGE , XPL OT , YPL OT 

COMMON  /BODCTR/  XCT< 102) , YCTC 1 02 ) , XBL ( 102) , YBL ( 1 02 ) , NCT 

COMMON  ITRK  (20000) 

DIMENSION  XI ( 10404 ), Y1 ( 1 0404 ),Z1(1 0404 ),ZZ( 200, 200) 
DIMENSION  U(102,102),V(102,102),T(102,102),P(102,102), 

2  IBODY(  102, 102) , ST  1(102, 102  > , ST2C  102, 102) 

DIMENSION  XU( 101 ) , YV( 101 ) 

REAL *8  XARRA Y( 102) 

CHARACTER*!  IDEC1 

LOGICAL  ZERO, CONT.ZERI, CONTI 


DEFINE  SOME  FUNCTIONS 

ZERO( 1 1 , J J )  =  ABS(U(II,JJ)).LT.0.1E-10  .AND. 

2  ABS(V(II,JJ)) .LT.0.1E-10 

CONT  ( 1 1 ,  J  J  )  =  (ZEP.  0(11  +  1  ,JJ  +  1  )  ,AND.ZERO(II+l  ,  J  J  —  1  )  .AND. 

2  ZER0(II-1, JJ+1). AND. ZERO (II-l, JJ-1)) 

3  . NEQV .  (ZERO(II+1,JJ+1).OR.ZERO(II+1,JJ-1).OR. 

4  ZERO(II-l, JJ+1) .OR.ZERO(Il-l.JJ-l)) 
ZERKII.JJ)  =  IBODYdl,  JJ)  .EQ.O 

CONTI  ( 1 1 ,  J  J  )  =  (ZERKII  +  l,  JJ  +  1).  AND.  ZERI  ( 1 1  + 1 ,  J  J-l )  .  AND. 

2  ZERKII-l,  JJ+1)  .  AND.ZERKII-l  ,  JJ-1)) 

3  .NEQV.  (ZERKII  +  l,  JJ  +  1)  .  OR  .  ZERI  ( 1 1  + 1 ,  JJ-1)  .OR. 

4  ZERKII-l,  JJ  +  1)  .OR.  ZERI  (II- 1,  JJ-1)) 

***  DEFINE  PAGESIZE  AND  FLOTSIZE  **X 

XPAGE  =  9.5 
YPAGE  =11.0 
XPIOT  =  6.5 
YPLOT  =  6.5 

***  READ  DATA  FOR  GRID  XXX 
870  REWIND  (3) 

READ  (3)  NT, LI, Ml 

READ  (3)  (XARRAY(I),I=2,L1) 

DO  8002  1=2, LI 
002  X'J(I)  =  SNGL  (  XARR  AY  (  I ) ) 

READ  (3)  (XARRAY( J) , J=2,M1 ) 

DO  8004  J=2,M1 

8004  YV ( J  )  =  SNGL ( XARRAY( J  ) ) 

READ  (3)  (XARR AY (I), 1=1, LI) 

DO  8006  1=1 , LI 
8006  X(I)  =  SNGL ( XARRAY( I  )  ) 

READ  (3)  (XARRAY(J),J=1,M1) 

DO  8008  J  =  1 , M 1 
8008  Y( J )  =  SNGL ( XARRAY( J  )  ) 

PREPARE  THE  GRID-DATA  FOR  THE  CONTOUR  PLOTTING 
NPTS  =  L 1 *M1 
DO  1300  J=1 , Ml 
DO  1300  1=1 , LI 
J1=L1*( J-l )+I 
1300  XI  (  J 1 )  =  X(I) 

DO  1310  J  =  1 , Ml 
DO  1310  1  =  1, LI 
J1=L1*( J-l  )  +  I 
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1310  Y1CJ1)  =  Y(J) 

C  ***  NOMINATE  DEVICE  *** 

C 

WRITE  (*,930) 

930  FORMAT  ( 1H1 , 39 ( 1 H*  )/lX, • PL  EASE ,  SELECT  PLOTTING  DEVICE •/ IX , 34( 1 K* 


2  //5X, 'TEKTRONIX  618  SCREEN  ....  (1)'/ 

3  5X, 'PRINTABLE  SHERPA  FILE  ...  (2)’/ 

4  5 X , 'IBM  COLOR  MONITOR  .  (3)'/ 

6  5X, 'COMPRESSED  METAFILE  .  (4)'/ 

5  5X, 'EXIT  PLOTTING  .  (99)  '  ) 

READ  (*,*)  IPLOT 

IF  <( IPLOT . LT . 1 ) .OR . (IPLOT .GT .4) >  GOTO  180 
IF  (IPLOT. EQ.l)  CALL  TEK618 


IF  (IPLOT. EQ.  2)  CALL  SHERPA  ('CONTOUR  ',’T’,2) 

IF  (IPLOT. EQ. 3)  CALL  IBM79 
IF  (IPLOT. EQ. 4)  CALL  COMPRS 
C 

C  ***  GET  THE  NUMBER  OF  TIME-STEPS  TO  BE  PLOTTED  AND  THE  INCREMENT  *** 
WRITE  (*,910)  NT 
910  FORMAT  ( 1H1 , 39( 1H*)/, 

2  IX, 'YOUR  INPUT  FILE  INCLUDES  ',13,'  TIME-STEPS '/l X . 39 ( 1H* )// 

3  IX, 'CHOOSE  A  METHOD  FOR  CONSTRUCTING  A  REGULAR  MATRIX:’/ 


4  5X, 'UNIFORM  GRID/NO  REFINEMENT  (102*102) .  (0)'/ 

5  5X , *  DI SSPL  A ' 1 S  GETMAT  ALGORITHM  (  1  02*1  02) .  (  1  )’/ 

6  5X, 'CUBIC  SPLINE  INTERPOLATION  (102*102) .  (2)'/ 

7  5X, 'LINEAR  INTERPOLATION  (200*200) .  (3)') 


READ  (*.*)  IMAT 
IF  (IMAT.NE.O)  GOTO  202 
LL  =  LI 
MM  =  Ml 
GOTO  204 

202  WRITE  (*,911  ) 

911  FORMAT  (IX, 'ENTER  THE  DIMENSIONS  OF  THE  MATRIX  TO  BE  CONSTRUCTED:' 

2  ,'  (IXD1M, IYDIM) ' ) 

READ  (*,*)  L  L , MM 

C  GENERATE  AN  EQUALLY  SPACED  GRID 
204  DEL X  =  X(L1)-X(1) 

DELY  =  Y(M1)-Y(1) 

DO  206  1=1, LL 

206  XXd)  =  Xd)  +  DELX*FL0AT(I-1)/FL0AT(LI-1) 

DO  208  J  =  1 , MM 

208  YY(J)  =  Yd)  +  DELY*FLOAT(  J-l  )/FLOAT(MM-l  ) 

C _ 

C  ***  LOOP  OVER  A  SERIES  OF  PLOTS  ***  BEGIN  OF  "200-LOOP" 

C 

DO  200  N=1 , NT 
READ  (3)  K 
WRITE  (*.912)  K 

912  FORMAT  (//IX, 'DO  YOU  WANT  ANY  PLOTS  FOR  THE  M3, 

2  '-TH  TIME-STEP?  (Y/N)’) 

READ  (*,914)  IDEC1 
914  FORMAT  (Al) 

IF  ( IDEC1  EQ. 'Y' )  GOTO  212 
C  ADVANCE  BY  1  TIME-STEP 
DO  210  1=1,6 
210  READ  (3) 

GOTO  200 

C  READ  DEPENDENT  VARIABLES  FROM  FILE  3 

212  DO  213  1=1.11 

READ  (3)  (XARRAY(J), J=1,M1) 

DO  213  J=1 . Ml 

213  Ud.J)  =  SNGL  (  XARRAY(  J )  1 
DO  214  1  =  1, LI 

READ  (3)  (XARRAY(J), J=1.M1> 

DO  214  J=1 ,M1 

214  V(I.J)  =  SNGL ( XARRAY ( J 1 ) 

DO  215  1  =  1, LI 

READ  (3)  (XARRAY(J),J=1,M1) 

DO  215  J=1 , Ml 

215  T (  I ,  J  )  *  SNGL ( XARRAY( J ) ) 

DO  216  1=1,11 

READ  (3)  (XARRAY(J), J  =  1 ,M1  ) 
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noon  no  noon 


DO  216  J=1 ,M1 

216  P(I,J)  -  SNGL ( XARRAYC  J  )  ) 

DO  217  1=1, LI 

READ  (3)  (XARRAYC J ) , J  =  1 ,M1 ) 
DO  217  J  =  1 , Ml 

217  ST1 ( I , J  )  =  SNGL (XARRAYC  J  ) ) 
DO  218  1=1, LI 


READ  (3)  (XARRAYCJ), J=1,M1 > 

DO  218  J=1,M1 

218  ST2( I , J )  =  SNGL ( XARRAYC  J ) ) 

DO  219  1=1, LI 

219  READ  (3)  C I BODYC I , J ) , J=1 , Ml ) 

IF  (N.NE.l)  GOTO  890 
WRITE  (*,939) 

939  FORMAT  C/1X, 'SHOULD  THE  CONTOUR  OF  THE  BODY  BE  IDENTIFIED  VIA:'/ 

2  6X, 'THE  VELOCITIES  .  Cl)’/ 

3  6X,  ’  THE  INPUT  DATA  .  (2)') 

READ  (*,*)  ICONT 

IF  (ICONT. EQ. 2)  GOTO  1500 


***  FIND  THE  CONTOUR  OF  THE  BODY  VIA  THE  VELOCITIES  *** 

LOCATE  A  CONVEX  CORNER  OF  THE  BODY 
DO  1320  J= A , Ml -3 
DO  1320  I=m,L1-3 
IF  C 2ER0C I , J ) )  GOTO  1322 
1320  CONTINUE 
1322  ICT  =  1 

XCTCICT)  =  XCI) 

YCTCICT)  =  YCJ) 

IP  =  1 
JP  =  0 

DO  1340  ICT=2 , 1 02 
DO  1324  IDUM=1,MAX(L1,M1) 

IF  (ZEROCI+IP, J+JP). AND. CONTCI+IP, J+JP))  GOTO  1330 
IF  (IP.EQ.O)  THEN 
JP  =  0 

DO  1326  IP=-1 , 1 , 2 

1326  IF  (ZEROCI+IP, J+JP) . AND . CONTC I+IP , J+JP))  GOTO  1332 
ELSE 

IP  =  0 

DO  1328  JP=-1 , 1 , 2 

1328  IF  CZEROCI+IP, J+JP). AND. CONTCI+IP, J+JP))  GOTO  1532 
EN  DI  F 

ANOTHER  CONTOUR  POINT  HAS  BEEN  FOUND 
1330  I  =  I+IP 
J  =  J+JP 
1324  CONTINUE 

THE  PREVIOUSLY  FOUND  POINT  IS  A  CORNER 
1332  XCTCICT)  =  XCI) 

YCTCICT)  =  YCJ) 

IF  (ABS(XCT(ICT)-XCT(1)).LT.0.1E-6.AND. 

2  ABSCYCT(ICT)-YCT(l)).LT.0.1E-6)  GOTO  1350 
I  =  I+IP 
J  =  J+JP 
1 340  CONTINUE 
1350  NCT  =  ICT 
GOTO  1400 


***  FIND  THE  CONTOUR  OF  THE  BODY  VIA  INPUT  *** 

LOCATE  A  CONVEX  CORNER  OF  THE  BODY 
1500  DO  1520  J  =4 , Ml -3 
DO  1520  I =4, L 1-3 
IF  ( ZERI ( I , J )  )  GOTO  1522 
1520  CONTINUE 
1522  ICT  =  1 

XCTCICT)  =  XCI) 

YCTCICT)  =  YCJ) 

IP  =  1 
JP  =  0 

DO  1540  ICT  =  2,102 
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on  on  no  on 


GOTO  1530 


DO  1524  IDUM* 1 , MAX( L 1 ,M1 ) 

IF  (ZERKI  +  IP, J+JP) .AND. CONTI (I+1P, J+JP ) ) 

IF  (IP.EQ.O)  THEN 
JP  =  0 

DO  1526  IP=-1,1,2 
1526  IF  (ZERKI  +  IP,  J  +  JP)  .AND. CONTMI  +  IP,  J  +  JP))  GOTO  1532 
ELSE 

IP  =  0 

DO  1528  JP=-1 ,1.2 

1528  IF  (ZERKI  +  IP,  J  +  JP). AND. C0NTK1  +  IP, J  +  JP))  GOTO  1532 
ENDIF 

C  ANOTHER  CONTOUR  POINT  HAS  BEEN  FOUND 
1530  I  =  I+IP 
J  =  J+JP 
1524  CONTINUE 

C  THE  PREVIOUSLY  FOUND  POINT  IS  A  CORNER 
1532  XCT(ICT)  =  XC I ) 

YCT(ICT)  =  Y(J) 

IF  (ABS(XCT(ICT)-XCT(l))  .LT.0.1E-6  .AND. 

2  ABS(YCT(ICT)-YCT(l)).LT.0.1E-6)  GOTO  1550 
I  =  I+IP 
J  =  J+JP 
1540  CONTINUE 
1550  NCT  =  ICT 

COMPUTE  THE  COORDINATES  FOR  BLANKING  (IN  INCHES  FROM  PHYSOR) 
1400  DO  1410  ICT= 1 , NCT 

XBLCICT)  =  £XCT(ICT)-X<1))/(X(L1)-X(1))XXPL0T 
1410  YBLCICT)  =  (YCT(ICT)-Y(1))/(Y(M1)-Y(1))XYPL0T 

xxx  FIND  OUT  WHICH  CURVES  TO  PLOT  xxx 
890  WRITE  ( x , 900 ) 

900  FORMAT  (lHl.'WHAT  WOULD  YOU  LIKE  TO  PLOT?'/ 


2  5X, 'TEMPERATURE .  (1)'/ 

3  5X, ’STREAMLINES .  (2)'/ 

4  5X, 'U-VELOCITY  .  C  3 ) *^ 

5  5X, 'V-VELOCITY  .  (4)’/ 

6  5X, ’PRESSURE  .  (5>V 

7  5X, ’NEXT  TIMESTEP  .  (0)’/ 

8  5X, ’  RE-INITIAL ISE  .  (-99)’/ 

9  5X, ’EXIT  PLOTTING  .  (99)  ’  ) 

READ  (x,x)  NTYP 


IF  (NTYP. EQ. 99)  GOTO  180 
IF  (NTYP. EQ. -99)  GOTO  870 
IF  (NTYP.EQ.O)  GOTO  200 

ASSIGN  THE  VALUES  OF  THE  Z-FUNCTION  WITH  THE  CHOSEN  DATA 
GOTO  (1100,1102,1104,1106,1100)  NTYP 

TEMPERATURE 

1100  DO  1101  J* 1 , Ml 
DO  1101  1=1 , LI 

1101  Z(I.J)  =  T( I ,  J ) 

GOTO  1120 

C  STREAMFUNCT1 ON 

1102  DO  1103  J=1 ,M1 
DO  1103  1*1,11 

1103  Z(I,J)  =  ST2( I , J )  +  .0005222 
GOTO  1120 

C  U-VELOCITY 

1104  DO  1105  J  =  1 , Ml 
DO  1105, 1=1, LI 

1105  Z(I,J)  =  U(I,J) 

GOTO  1120 

C  V-VELOCITY 

1106  DO  1107  J*1,M1 
DO  1107, 1*1, LI 

1107  Z(I,J)  *  V(I,J) 

GOTO  1120 

C  PRESSURE 

1108  DO  1109  J  =  1 , Ml 
DO  1109  1*1, LI 
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no  on  non  no  non  o *  * x *  *  *  * *  oooo  o  no 


1109  Z(I,J)  =  P(I,J) 

COMPUTE  AND  DRAW  CONTOUR-LINES 
1120  CALL  CONTDD  ( NT YP , IHAT , XI , Y1 , Z1 , ZZ , L 1 , Ml , L 1 *M1 , L L , MM ) 
GOTO  890 

200  CONTINUE 

END  OF  "200-L  OOP" 


xxx  EXIT  PLOTTING  ROUTINE  **X 
130  CALL  DONEPL 
STOP 
END 

:xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx> 

SUBROUTINE  CONTDD  (NTYP , IMAT, XI , Y1 , Z1 , ZZ , L 1 , Ml , NTTS, LL , MM) 

THIS  SUBROUTINE  COMPUTES  A  REGULAR  MATRIX  (ON  AN  EQUALLY  SPACED 
GRID)  EITHER  BY  DISSPLA'S  GET MAT-ALGORITHM  OR  BY  A  CUBIC  SPLINE 
INTERPOLATION.  SUBSEQUENTLY  CONTOUR-LINES  ARE  DRAWN. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx« 

COMMON  /BODCTR/  XCT ( 1 02 ) , YCT( 1 02 ) , XBL ( 1 02 ) , YBL ( 1 02 ) , NCT 
COMMON  /INDATA/  XC 1 02 ) , Y( 1 02 ) , Z< 1 02 , 1 02 ) 

COMMON  /GRIDII/  XX( 200 ) , YY( 200 ) 

COMMON  /PLOTDT/  XPAGE, YPAGE, XPL OT, YPLOT 
COMMON  WRK  (20000) 

DIMENSION  X1(NPTS),Y1(NPTS),Z1(NPTS),ZZ(LL,MM) 

DIMENSION  21(200,200) ,C<4, 102) 

CHARACTERXl  LCONT , LGRI D 


xxx  INPUT  THE  INCREMENT  BETWEEN  CONTOUR  LEVELS  xxx 
WRITE  (*,900) 

900  FORMAT  (IX, 'SHOULD  DISSPLA  PICK  THE  INCREMENT  BETWEEN  CONTOUR  ', 
2  'LEVELS?  (Y/N)') 

.  READ  (*,910)  LCONT 

910  FORMAT  (Al) 

IF  (LCONT. EQ. 'Y')  GOTO  925 
WRITE( x,  920  ) 

920  FORMAT  (IX, 'INPUT  THE  INCREMENT  BETWEEN  CONTOUR  LEVELS:') 

READ  (*,*)  ZINCR 

925  WRITE  ( x , 930 ) 

930  FORMAT  (IX, 'DO  YOU  WANT  A  GRID  TO  BE  DRAWN?  (Y/N)’) 

READ  (*,910)  LCRID 

xxx  SET  UP  ALPHABET  xxx 
CALL  RESET  CALL') 

CALL  SW1SSM 
CALL  HWSHD 
CALL  CHRPAT ( 16  ) 

CALL  BASAL  F  ( 'L/CSTD' ) 

CALL  MIXALF  ( 'STAND' ) 

CALL  SETCLR  ( 'BLUE'  ) 

CALL  SHDCHR( 90 . ,1, .002,1) 

xxx  DEFINE  PLOT  AREA  xxx 

CALL  PAGE  (XPAGE, YPAGE) 

IF  (IPLOT.EQ.l)  CALL  HWSCAL  ('SCREEN') 

CALL  PHYSOR  (1.25,1.0) 

CALL  AREA2D  ( YPLOT , XPL OT ) 

xxx  TEMPORARY  FIX  FOR  CONTOURING  PROBLEMS 
1  =  10 
IF  1  =  1 

CALL  THKFRM  ( .01) 

xxx  LABEL  AXES  xxx 

CALL  HEIGHT  (0.18) 

CALL  SETCLR  ('GREEN') 
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ooo  nn  ooon  non  non  o  on 


IF  (NTYF.E8.1)  CALL  HEADIN  ( *  £ T ) EMPERAT URE  ( C. 10NT0URS* ’ , 
1  100,1.5,1) 

IF  (NTYP.EQ.2)  CALL  HEAD1I1  ( ' (S)TKEAMLINE3$ ' » 

1  100,1.5,1) 

IF  (NTYP.EQ.3)  CALL  HEADI N  ( • ( U- V ) EL OCI TY  ( C  )0NT0UR5S ' , 

1  100,1.5,1) 

IF  (NTYP.EQ.9)  CALL  HEADIN  (  •  ( V- V  )EL  OCIT  Y  (C)OIITOURSS  '  , 

1  100,1.5,1) 

IF  (NTYP.EQ.5)  CALL  HEADIN  ('(P)RESSURE  (C)ONTOURSS' 

1  ,100,1.5,1) 

CALL  RESET  ( 'HEIGHT' ) 

CALL  SETCLR  ('BLUE') 


XXX  SET  UP  AXES  *xx 

CALL  XNAME  < 'Y  $ '  ,  100) 

CALL  YAXANG  (0.0) 

CALL  YNAME  ( 'X  $  '  , 100 ) 

XINC  =  0.1*(X(L1)-X(1)) 

YINC  =  0 . 1X(Y(M1 )-Y( 1 ) ) 

CALL  GRAF  £ Y ( 1 ), YINC , Y( Ml ), X( 1 ), XINC , X( L 1 ) ) 

CALL  GRID  (0,0) 

CALL  BCOMON  (20000) 

DRAW  A  GRID  (IF  REQUESTED) 

IF  (LGRID.EQ. 'N' )  GOTO  90 
DO  30  1  =  1, LI 

30  CALL  RLVEC  ( Y( 1 ) , X( I ) , Y( Ml ) , X( I ) , 0000 ) 

DO  32  J  =  1 , Ml 

32  CALL  RLVEC  ( Y( J ) , X< 1 ) , Y( J ) , X( Ml ) , 0000 ) 

XXX  PREPARE  CONTOURS  PLOTS  xxx 
90  GOTO  (100,200,300)  IMAT 

xxx  INPUT  MATRIX  IS  REGULAR  /  NO  GRID  REFINEMENT  REQUIRED  xxx 

DO  10  J=1 , MM 
DO  10  1  =  1, LL 
10  ZZ(I,J)  =  Z(J,I) 

GOTO  900 

XXX  CONSTRUCT  THE  MATRIX  BY  DISSPLA’S  ALGORITHM  GETMAT  xxx 

TRANSFORM  THE  Z  2-D  FUNCTION  TO  1-D  FUNCTION 
100  DO  110  J  =  1 , Ml 
DO  110  1  =  1. LI 
J1  =L 1 *( J-l >+I 
110  Zl(Jl)  =  Z(I,J) 

CONSTRUCT  A  REGULAR  MATRIX 
CALL  BGNMAT ( MM, L L ) 

CALL  GETMAT(Y1,X1,Z1,NPTS,0) 

CALL  ENDMAT (ZZ, 0) 

GOTO  900 

xxx  CONSTRUCT  THE  MATRIX  BY  CUBIC  SPLINE  INTERPOLATION  xxx 

200  DO  210  J= 1 , Ml 
DO  212  1=1,11 
212  C C 1 , I )  =  Z(I,J) 

CALL  CUBSPL  (X,C,L1,0,0) 

DO  219  1  =  1, LL 

219  ZI(I,J)  =  PPVALU  (X,C,L1,9,XX(I),0) 

210  CONTINUE 

DO  220  1  =  1, LL 
DO  222  J  =  1 , Ml 
222  C ( 1 , J )  =  ZI(I,J) 

CALL  CUBSPL  (Y,C,M1,0,0) 

DO  229  J  =  1 , MM 

229  ZZ ( J , I )  =  PPVALU  ( Y , C , Ml , 9 , YY( J ) , 0 ) 

220  CONTINUE 
GOTO  900 
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c 

C  ***  CONSTRUCT  THE  MATRIX  BY  LINEAR  INTERPOLATION  *** 

C 

300  DO  310  J  =  1 .Ml 
DO  310  1=1 . U 

CALL  INTERV  ( X, L 1 , XXC I ) , L EFT , MPL AC ) 

310  ZI(I,J)  =  ZCLEFT,J)  +  < Z< L EFT+ 1 , J  ) -Z< L EFT , J ) )*( XXC I ) -XC L EFT ) )/ 

2  (XCLEFT»1)-X(LEFT)) 

C 

DO  320  1  =  1, LL 
DO  320  J=1 , MM 

CALL  INTERV  ( Y ,M1 , YYC J ) , LEFT , MFL AG ) 

320  ZZ(J,I)  =  ZI ( I , LEFT )  +  ( ZI C I , L EFT + 1 ) -ZI < I , L EFT ) )* C YYC J ) -Y< L EFT ) )/ 

2  (YCLEFT  +  D-YCLEFT)) 

C 

C  ***  DRAW  THE  BODY  CONTOUR  AND  BLANK  THE  INSIDE  *** 

400  CALL  LINEAR 

CALL  SETCLR  C'CYAN') 

CALL  THKCRV  CO. 025) 

CALL  CURVE  C YCT , XCT , NCT , 0 ) 

CALL  RESET  C 'THKCRV') 

IF  (NTYP.EQ.l  .OR.  LGRID.E9. »Y' )  GOTO  410 
C  CALL  SHADE  C YCT , XCT , NCT , 45 . ,  .  1 ,  1 , 0 , 0  ) 

CALL  BLPOLY  C YBL , XBL , NCT , 1 ) 

C 

C  ***  COMPUTE  THE  CONTOUR-LINES  *** 

410  CALL  HEIGHT  CO. 12) 

CALL  SETCLR  C 'RED' ) 

IF  (LCOHT.EQ.'Y')  CALL  CONMAK  t ZZ , MM, L L , 5HSCAL E ) 

IF  C  LCONT . EQ  .  ' N ' )  CALL  CONMAK  C ZZ , MM, L L , ZI NCR ) 

CALL  COHLINCO,5HSOLID,6HLABELS, 1 ,5) 

C 

C  ***  DRAW  CONTOUR  LINES  *** 

CALL  CONANGC 1 35 . ) 

CALL  CONMIMC 1 .5) 

C  CALL  RASPLNC  0.25) 

CALL  C0NTURC1,6HLABELS,4HDRAW) 

CALL  RESET  C 'THKCRV') 

C 

C  TERMINATE  THE  CURRENTLY  DISPLAYED  GRAFH 
CALL  ENDGRCO) 

CALL  ENDPLCO) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C  c 
SUBROUTINE  CUBSPL  C TAU, C , N, IBCBEG, IBCEND ) 

C  ( 
C  THIS  SUBROUTINE  COMPUTES  THE  CUBIC  SPLINE  INTERPOLANT  FOR  C 
C  DIFFERENT  BOUNDARY  CONDITIONS.  C 
C  SOURCE:  CARL  DE  BOOR,  A  PRACTICAL  GUIDE  TO  SPLINES  C 
C  C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c  TAU  .  ABSCISSAE  OF  THE  DATA  POINTS 

C  C(1,I)  ....  ORDINATES  OF  THE  DATA  POINTS 

C  N  .  NUMBER  OF  DATA  POINTS 

C  IBCBEG  ....  INDICATOR  OF  BOUNDARY  CONDITION  AT  TAU(l) 

C  I BCBEG=  0 :  NOT-A-KNOT-CONDITION 

C  IBCBEG=1:  PRESCRIBED  SLOPE 

C  IBCBEG=2 ■  PRESCRIBED  SECOND  DERIVATIVE 

C  IBCEND  ....  INDICATOR  OF  BOUNDARY  CONDITION  AT  TAU(N) 

C  CCJ.I)  ....  POLYNOMIAL  COEFFICIENTS  OF  SPLINE  INTERPOLANT 

C  PCX)  =  C(l,I)+H*(C(2,I)+HX(C(3,I)+H*C(4,I)/3.)/2.) 

C  WHERE  H  =  X  -  TAUCI) 

C  _ 

C 

DIMENSION  C(4,N),TAU(N) 

C  _ 

C 

L  *  N-l 

C  COMPUTE  FIRST  DIFFERENCES  AND  FIRST  DIVIDED  DIFFERENCES 
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non  no  no  no  non  o  non  no  no  on  o  o  ooo 


C  CC  SERVES  AS  A  TEMPORARY  STORAGE) 

DO  10  M=2,N 

CC  3 , M )  =  TAU(M)  -  TAU( M-l  ) 

10  CC4,M)  =  CCCl.M)  -  CC1,M-1))/CC3,M) 

XXX  CONSTRUCT  FIRST  EQUATION  FROM  FIRST  BOUNDARY  CONDITION  *** 

IF  (IBCBEG-1)  11,15,16 

11  IF  CN.GT.2)  GOTO  12 

NOT-A-KNOT  COND  AT  LEFT  END,  AND  N.EQ.2 
CCA, 1  )  =  1.0 
CC  3 , 1 )  =  1.0 
CC  2 , 1 )  =  2 . 0*CC  4, 2) 

GOTO  25 

NOT-A-KNOT  COND  AT  LEFT  END,  AND  N.GT.2 

12  CC  4 , 1 )  =  CC  3 , 3 ) 

CC  3 , 1 )  =  CC  3 , 2 )  +  CC  3y  3  ) 

CC  2 , 1 )  =  CCCC3,2)+2.0*C(3,1))*CC4,2>*CC3,3)+CC3,2)K*2*CC4,3)) 

2  /CC  3, 1 ) 

GOTO  19 

SLOPE  PRESCRIBED  AT  LEFT  END 

15  CC  4, 1 )  =  1.0 
CC  3, 1 )  =  0.0 
GOTO  18 

SECOND  DERIVATIVE  PRESCRIBED  AT  LEFT  END 

16  CC  4 , 1 )  =  2.0 
CC  3, 1 )  =  1.0 

CC 2 , 1 )  =  3 . 0*CC  4 , 2 )  -  CC3,2)/2.0*CC2,1) 

18  IF  CN.E0.2)  GOTO  25 

IN  CASE  OF  INTERIOR  KNOTS  SET  UP  EQU,  EXECUTE  FORWARD  PASS  OF  GAUSS 

19  DO  20  M=2, L 
G  =  -CC3,M+1)/CC4,M-1) 

CC  2, M)  =  GKCC2.M-1)  +  3 . 0*C CC  3 , M' VC  • , M+ 1 HCC 3, M+l ) *CC 4 , M ) ) 

20  CC  4 , M)  =  G*CC  3 , M-l )  +  2.0*(C(T,M)  +  C(3,M+1)) 

***  CONSTRUCT  LAST  EUQATI0N  FROM  SECOND  BOUNDARY  CONDITION 
IF  CIBCEND-1)  21,30,24 

FOR  PRESCRIBED  SLOPE  ENTER  DIRECTLY  THE  BACK  SUBSTITUTION 

21  IF  CCN.EQ.3).AND.CIBCEEC  EQ  0))  GOTO  22 

NOT-A-KNOT  COND  AT  RIGHT  END,  AND  N.GE.3, 

AND  CCN.GT.3)  .OR.  C  NOT-A-KNOT  COND  AT  LEFT  END)) 

G  =  CC  3, N-l  )  +  CC  3 , N  ) 

C  C  2 , N )  =  CCCC3,N)  +  2.0*G)*CC4,N)*CC3,N-1) 

2  +  CC3,N)*K2KCCCl,N-l)-CCl,N-2))/CC3,N-l))/G 

G  =  -G/CC4.N-1) 

C C  4 , N )  =  CC  3, N-l ) 

GOTO  29 

NOT-A-KNOT  COND  AT  RIGHT  END,  AND  C C N . EQ . 3) . AND . C NOT-A-KNOT  AT  LEFT) 
.OR.  CC N.EQ.2) . AN D.C NOT-A-KNOT  COND  AT  LEFT  END)) 

22  C  C  2 , N  )  =  2 . 0  * C  C  4 , N ) 

CC  4 , N  )  =  1.0 
GOTO  28 

SECOND  DERIVATIVE  PRESCRIBED  AT  RIGHT  END 

24  C C  2 , N )  =  3 . 0*C  C  4 , N )  +  C ( 3 , N )/2 . 0 *C C 2 , N ) 

C  C  4 , N )  =  2.0 
GOTO  28 

25  IF  CIBCEND-1)  26,30,24 

26  IF  CIBCBEG.GT.O)  GOTO  22 

NOT-A-KNOT  COND  AT  RIGHT  AND  LEFT  END,  AND  N.EQ.2 
C C  2 , N )  =  CC4.N) 

GOTO  30 

28  G  =  -1  .0/CC4.N-1) 

xxx  SOLVE  TRIDIAGONAL  SYSTEM  FOR  THE  SLOPES  CC2,J)  xxx 
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C  FORWARD  PASS  ^ELIMINATION  OF  LOWER  DIAGONAL) 

29  C(4,N)  =  GKC(3,N-1)  +  C(A,N) 

C(2,N)  =  (G*C(2,N-1)  +  C(2,N))/C(4,N) 

C  BACKWARD  PASS  ( =  BACK  SUBSTITUTION) 

30  DO  AO  J=L, 1,-1 

AO  C(2,J)  =  ( C( 2, J )  -  C(3,J)*C(2, J+1))/C(A, J) 

C  _ 

c 

C  ***  GENERATE  CUBIC  COEFFICIENTS  IN  EACH  INTERVAL 
C 

DO  50  1=2. N 
DTAU  =  C( 3, 1  ) 

DIVDF1  =  <C(1,I)  -  C( 1 , 1-1 ) )/DTAU 
DIVDF3  =  C(2,I-1)  +  CC  2, 1 )  -  2.0*DIVDF1 
C( 3, 1-1 )  =  2.0*(DIVDF1-C(2,I-1)-DIVDF3)/DTAU 
50  C(A,I-1)  =  (DIVDF3/DTAU)*(6 .O/DTAU) 

C 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

REAL  FUNCTION  PPVALU  ( BREAK , COEF, L , K, X, JDERIV ) 

C 

C  THIS  FUNCTION  CALCULATES  THE  VALUE  OF  THE  JDERIV-TH  DERIVATIVE 
C  OF  A  PIECEWISE  POLYNOMIAL  FUNCTION. 

C  SOURCE  j  CARL  DE  BOOR,  A  PRACTICAL  GUIDE  TO  SPLINES 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


c 

C  BREAK  .  ABSCISSAE  OF  BREAKPOINTS 

C  COEF  .  COEFFICIENTS  OF  PIECEWISE  POLYNOMIAL  INTERPOLANT 

C  L  .  NUMBER  OF  BREAKPOINTS 

C  K  .  ORDER  OF  POLYNOMIAL  FUNCTIONS 

C  X  .  ABSCISSA  AT  WHICH  TO  EVALUATE 

C  JDERIV  .  ORDER  OF  DERIVATIVE  (0  =  FUNCTION) 

C  _ 

C 


DIMENSION  BREAK(L),COEF(K,L) 

C  _ 

c 

PPVALU  =  0.0 
FMMJDR  =  K  -  JDERIV 

C  DERIVATIVES  OF  ORDER  K  OR  HIGHER  ARE  IDENTICALLY  ZERO 
IF  (FMMJDR. LE. 0.0)  GOTO  99 
C 

IF  CL.EQ.l)  THEN 

C  SINGLE  POLYNOMIAL  FUNCTION 
I  =  1 
ELSE 

C  FIND  INDEX  I  OF  LARGEST  BREAKPOINT  TO  THE  LEFT  OF  X 
CALL  INTERV  ( BREAK , L , X , I , NDUMMY ) 

ENDIF 

H  =  X  -  BREAK(I) 

C 

C  EVALUATE  JDERIV-TH  DERIVATIVE  OF  I-TH  POLYNOMIAL  PIECE  AT  X 
DO  10  M=K, JDERIV+1,-1 

PPVALU  =  ( PPVALU/ FMMJDR )*H  +  COEFCM.I) 

10  FMMJDR  =  FMMJDR  -  1.0 
C 

99  RETURN 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  INTERV  (XT, LXT,X, LEFT.MFLAG) 

c 

C  THIS  SUBROUTINE  DETERMINES  THE  INTERVAL  INDEX  WITH  RESPECT 
C  TO  THE  BREAKPOINT  SEQUENCE. 

C  SOURCE:  CARL  DE  BOOR,  A  PRACTICAL  GUIDE  TO  SPLINES 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c  XT  .  ABSCISSAE  OF  BREAKPOINTS 
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nnonoOOU  OOU  O  (JO  OO  OO  O  UUUU 


LXT  .  NUMBER  OF  BREAKPOINTS 

X  .  ABSCISSA  WHOSE  INTERVAL  INDEX  IS  TO  BE  DETERMINED 

LEFT  .  INDEX  OF  THE  LARGEST  BREAKPOINT  TO  THE  LEFT  OF  X 

MFLAG  .  RETURN  CODE:  MFLAG  =-l :  X.LT.XT(l) 

MFLAG  =  0:  XT(1) .LE.X.LT.XT(LXT) 
MFLAG  =  1 :  X.GE.XT(LXT) 


DIMENSION  XT C  LXT ) 


RECALL  THE  INDEX  OF  THE  PREVIOUS  CALL  TO  THIS  SUBROUTINE 
SAVE  ILO 
DATA  ILO/1/ 

I  HI  =  ILO  +  1 
IF  ( I  HI  .LT.LXT)  GOTO  20 
IF  C  X . GT . XT ( LXT ) )  GOTO  110 
IF  (LXT.LE.l)  GOTO  90 
ILO  =  LXT  -1 
IHI  =  LXT 

20  IF  (X.GE.XT(IHI))  GOTO  40 
IF  (X.GE.XT(ILO))  GOTO  100 

DECREASE  ILO  TO  CAPTURE 
ISTEP  =  1 
31  IHI  =  ILO 

ILO  =  IHI  -  ISTEP 
IF  tILO.LE.13  GOTO  35 
IF  (X.GE.XT(ILO))  GOTO  50 
ISTEP  =  ISTEP*2 
GOTO  31 
35  ILO  =  1 

IF  (X.LT.XTCl))  GOTO  90 
GOTO  50 

INCREASE  IHI  TO  CAPTURE  X 

40  ISTEP  =  1 

41  ILO  =  IHI 
IHI  =  ILO  +  ISTEP 
IF  (IHI .GE.LXT)  GOTO  45 
IF  (X.LT.XT(IHI)l  GOTO  50 
ISTEP  =  ISTEP*2 
GOTO  41 

45  IF  ( X . GE . XT ( LXT I )  GOTO  110 
IHI  =  LXT 

NARROW  THE  INTERVAL  ( XT< I LO ) , XTC IHI  )  ) 

50  MIDDLE  =  (ILO  +  I H I )/ 2 

IF  (MIDDLE. EQ. ILO)  GOTO  100 

IF  ( X . L  T . XT ( MI  DDL  E ) )  GOTO  53 
ILO  =  MIDDLE 
GOTO  50 

53  IHI  =  MIDDLE 
GOTO  50 

***  SET  OUTPUT  AND  ASSIGN  RETURN  CODES  *** 

X  LIES  TO  THE  LEFT  OF  THE  SMALLEST  BREAKPOINT 
90  MFLAG  =  -1 
LEFT  =  1 
RETURN 

C  REGULAR  RETURN  FROM  THIS  SUBROUTINE 
100  MFLAG  =  0 
LEFT  =  ILO 
RETURN 

C  X  LIES  TO  THE  RIGHT  OF  THE  LARGEST  BREAKPOINT 
110  MFLAG  -  1 

LEFT  =  L  XT- 1 

RETURN 

END 
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APPENDIX  C 
TRAVERSE  PROGRAM 


10  M  M  !!!!!  I  !!!!!!!!!!!!!!!!  I  !!!!!!  j  !  j  I  !!!!!!!!{)  I  !<!  I  !!!  M  !!!!  I  I  M  M  H 
20  !  !  PROGRAM  TRAV 

30  ! !  PROGRAM  USED  TO  MOVE  TRAVERSE 

31  !  ! 

32  !  ! 

40  I | I  I !)! I !!!!!!! M  !!!!!!!!!!!!!!!!!!  j  j  !{  j  !l  If!  j  j  !  Ml  f  !!  I  !!!  t  <  f  H  *  !!!  H 

42  Isc-9 

43  ASSIGN  @Datacoi*im  TO  Isc 

44  ASSIGN  STerminal  TO  1 

4G  CONTROL  Isc, 3i 1200 

47  CONTROL  I  sc ,4 ; 2+4+0+ 1 G 

49  CONTROL  Isc  ,5; 10 

50  COM  /Path/  @Datacomm 

52  OUTPUT  @Datacomm; "EZT" 

53  OUTPUT  SDataconm; "C  S 1 000  .SM2000 , 1 1 500  , IM1 500  ,R “ 

54  END 
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APPENDIX  D 


TEMPERATURE  ACQUISITION  PROGRAM 

!  i !  m  m  m  ! !  m m m m m  m m m M  i ii  m  ! i  m  i  m 1 1 1  m  m  ! 

!  !  TEMPERATURE  ACQUISITION  PR06RAM 

M  (STEADY  STATE) 

M  PROSRAM  PERFORMS  THE. REQUIRED  FUNCTIONS  NEEDED  TO  OBTAIN  DATA 
! I  FROM  THE  ACQUISITION  SYSTEM  AND  CONUERT  THE  THERMOCOUPLE  QUT- 
M  PUT  U0LTA6ES  INTO  TEMPERATURE  UALUES.  THE  TEMPERATURES  ARE 
I !  STORED  INTO  DATA  FILES  FOR  FURTHER  REDUCTION. 

I  UARIASLE  DEFINITIONS 
! !  UOLTS  :  THERMOCOUPLE  OUTPUT 

i  !  TEMP  :  HEATER  TEMPERATURES 

M  >  M  !  I  M  M  H  M  !  M  M  M  M  M  M  M  !  M  M  M  M  1  M  1  M  !  M  M  M  M  M  !  H  M  M  H  < 
REAL  Ucits(EO) 

REAL  Temp (50) 

ICREATE  BOAT  "THREDI5:  700,!"  45 
! ASSISN  @PatH3  TO  “THRED15:  700  ! " 

PP2NT  OATE$( TIMEDATE ) 

PRINT  “  " 

PRINT  “  1.0  UATT  INPUT  POWER “ 

PR Z NT  "  45  STRIPS  NO  SHROUD" 

PRINT  *  SDAT  FILE  THREDI5:,700  !“ 

PQTKjT  "  " 


£9  PRINT  "  ” 

93  ||  i  m  i  j  i  i  i  i  j  i  i i  m i i ii i  |  i  m  i  i  SURFACE  t !  -20 

95  OUTPUT  709;  "CONFMEAS  DCU  .  !00-M9,USE  0" 

QC  C Qp  T  m 1  TO  20 

99  ENTER  709;Uoits( I  ) 

59  Temp' I  >- .0005797+! 25925. ! 329*Uol ts< I ) >-( 607799. 2457* ( Uel t s(  I )*UcI t s 
2 !  952034 .3354* ( 'Jolts',  I  M3  )  >+(  9370910995.  !8?4»(Uciis<  I  )"4  )  >+  .039 

100  NEXT  I 

520  |  |  |  |  !  i  i  i  i  m  i  |  j  m  i  j  i  !  !  |  |  i  i  i  i  SURFACE  i2 1  -140 

550  OUTPUT  709; "CONFMEAS  DCU  ,200-2 1 9  .USE  0" 

550  FOR  1*2 1  TO  40 

570  ENTER  709 ; Ucl t  s ( I  > 

590  Temp'  I  )-. 000E797+! 25925 .  1328H'olts<  I  )  >-(607789. 246?*<Ucltsi  I  .•*ue:u 
21952034 .33S4*( Ucl ts( I  >A3> )+( 8370810995. 1 874«<Uolis< I )“4 ) >+.039 
590  NEXT  I 

600  I  M  !  I  l  M!  I  I  Ml  I  I  I  I  !  !  M  I  I  !  I  SURFAC  14 1  -445 
610  OUTPUT  709; "CONFMEAS  DCU ,300-304  .USE  0" 

520  FOR  1-41  TO  45 

856  ENTER  709;Uoits< I > 

857  Temp' I >-. 0006797+( 25825 . 1 329*Uol t s< I >  >-( 607799. 2467* ( Ucl tsv I )*Uoi-ts 
2! 952034.3364* (Celts'  I  )"3 >>+( 8370810996 . 1 874*< Uolts( I >“4  > >+.039 

859  NEXT  I 

979  III  PRINT  OUTPUT  TEMPERATURES 

879  PRINT  '  “ 

880  FOR  1-1  TO  45 

88!  PRINT  "HEATER  TEMPERATURE  IN  DES.  C  -";Tenp(I> 

692  NEXT  I 

893  'OUTPUT  0Path3;Temp( * > 

894  PRINT  "  " 

885  PRINT  "  " 
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APPENDIX  E 

UNIFORM  POWER  PROGRAM 


1  1  l  l  M | i l ll | !  |  |  |  l  |  i  l  |  !  i  M  !  l  I  !  i ! | ! I  M ! I ! ! ! ! I l M !  I  *  M  !  i  !  1  1  '  1  1  1  1  1  1  '  m 

2  i  !  !  POWER  INPUT  PROGRAM  '  '  ! 

4  i  i  M  i|  |  M  (  M  !  U  (  |  I  (  |  (  (  |  I  (  I  M  I  1  !  !  !  !  M  !  !  !  !  !  I  M  M  !  M  !  I  !  !  »  1  1  I  I  '  i  1  I  •  *  ' 

5  (OUTPUT  709s "LOCAL  CLR" 

7  OUTPUT  709 i “RST " 

10  REAL  Uolts<51  ) 

20  REAL  Pou(5I  ) 

30  REAL  Ptot  ,Pavo  , Perth  f 

920  OUTPUT  709 1 "CONFMEAS  DCU  ,308-31 5  ,USE  0” 

930  FOR  1-1  TO  8 

931  J-9-I 

940  ENTER  709;Volts(J  ) 

950  PRINT  "  “ 

952  NEXT  I 

953  OUTPUT  709s "RST" 

954  OUTPUT  709s "CONFMEAS  DCU,316,USE  0“ 

955  ENTER  709sUolts<  46  ) 

955  OUTPUT  709s "RST" 

971  OUTPUT  709s "CONFMEAS  OCU  .403  -4 1 7  .USE  0" 

972  FOR  I-S  TO  23 

973  ENTER  709sUolt3< I  ) 

977  NEXT  I 

976  OUTPUT  705s "RST" 

979  OUTPUT  709s "CONFMEAS  DCU  .400-40 1  .USE  0" 

980  FOR  1-47  TO  48 

961  ENTER  709sUolts( I  ) 

982  NEXT  I 

983  OUTPUT  709s "RST" 

985  OUTPUT  709s "CONFMEAS  DCU  ,503-51 7  .USE  0" 

985  FOR  1-24  TO  36 
987  ENTER  709sUolts(  I  ) 

989  NEXT  I 

990  OUTPUT  709s "RST" 

991  OUTPUT  709s "CONFMEAS  DCU ,500-50 1  .USE  0" 

992  FOR  1-49  TO  50 

993  ENTER  709sUolt3<  I  ) 

994  NEXT  I 

996  OUTPUT  709s "RST” 

997  OUTPUT  709s "CONFMEAS  DCU  .602-608  .USE  0" 

998  FOR  1-39  TO  45 

999  ENTER  709sUolts< I  ) 

1000  NEXT  I 

1001  OUTPUT  709s "RST" 

1002  OUTPUT  709s "CONFMEAS  DCU, 600  .USE  0* 

1003  ENTER  709iUolts<5! ) 

1004  OUTPUT  709s "RST" 

1005  PRINT  "  * 

1006  PRINT  *  " 

1050  PRINT  "INPUT  U0LTA6E  D.C.  UOLTS  PAN  1"sUolts(46> 

1051  PRINT  "INPUT  UOLTAGE  D.C.  UOLTS  PAN  2"iUoits<47) 

1052  PRINT  "INPUT  UOLTAGE  D.C.  UOLTS  PAN  3 " * Uo 1 1 s ( 48  ) 

1053  PRINT  "INPUT  UOLTAGE  D.C.  UOLTS  PAN  4"sUolts<49> 
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1054  PRINT  'INPUT  VOLTAGE  D.C.  VOLTS  PAN  S'«Volts<50> 

1055  PRINT  'INPUT  VOLTAGE  D.C.  VOLTS  PAN  6" iVolts(Sl  > 

10E0  PRINT  “  ' 

1070  PRINT  '  “ 

1080  Resist-2.0 
1109  R2*Resist*Resist 

1119  Pow(2  )-<  Volts(46  )-Volts<2  )  )*(  Volts(2  )«< Resist +  .4  )-Volts(46  )• .4  )/R2 
1 149  Pow<5  )-(Volts<46  )-Volts<5  >  )«<Volts(5  >*(Resist+.2 )-Voits(4S )*.2 )/R2 
1 179  Pow(9)-(Volts<47>-Volts<9}  )•< Vol t s< 9 >*( Resist+ .3 >-Vol ts( 47 >• .3 >/R2 
1189  Pou( 1 0  )-( Volts<47  )-Volta( 10  )  >«( Volta< 1 0 >*< Reaist+ .45 )-Volts(  47 >• .451/R2 
1 194  Pou< 18 >-( Volts(48  )-Volts< 1 8  )  >*<  Volta< 1 8  >*( Resist*. 55 )-Voits(  48 >• .55 )/R2 

1197  Pow( 23  )-( Volt a(48  )-Volts<23  )  )«(Volts(23  >*(Resist+ .55 >-Vol ts( 48 >*.55  )/R2 

1 198  Pou( 24 )-( Volt 3(49  )-Volts<24  >  >•( Volt s(24  >*(Resist+. 1  )-Vol ts< 48 >* . 1  >/R2 

1199  Pow(2S>-(Volt3(49>-Volts(25n*(Volts(25>*(Resist  +  .3>-Volts(48>».3>/R2 
1238  Pou( 29)-(Volts(49  )-Vol ts<  29>)*( Vo lts(29>*( Resist*. 45  )-Vol ts( 49)*.45)/R2 
1241  Pow<  32  )-< Volt 3(50  )-Volts( 32  )  )«(Volts(32  >*<Resist+.E )-Volts(50 )*.6  )/R2 

1243  Pow<  34  )»( Vol ts( 50  )-Vol t s( 34  >  >*(Volts(34  >•< Resist*. 3 )-Vol t s( 50 )* . 3 )/R2 

1244  Pou<  35  )-( Volt a(50  )-Volt s(35  )  >*(Volts(35 >«(Resist+.B >-Volis<50  >*.6  5/R2 

1245  Pow(3G  )-( Vol  t  s<  50  )-Vol ts( 3B )  )*(Volts(3B  )*(Resi3t+.B >-Volts( 50 >*  .6 )/R2 
124E  Poui<  37  )«( Vol  ts(  50  )-Vol ts<  37  ))*(  Vo  lts(  37  )•(  Resist^. B  )-Volts(  50  >•  .5  )/R2 

1247  Pou<  39  >■<  Vol ts<  5 1  )-Vol ts< 39  )  )*( Vol t s( 39  )«( Resist*. 3 >-Vol t s(  51  )*.3>/R2 

1248  Pou(40)-(Volts<51  )-Volts(40))«{Volts(40)«(Resi5t  +  .3)-Volt3(51  >*.3>/R2 

1250  Pou( 42  )-( Volt s(S1  )-Volts<42  )  >*(Volts<42 >*(Resist+.55 >-Volts( 51  >*.55  >/R2 

1251  Pow( 43  )-<  Vol 1 5 ( 51  )-Volts<43>>«(Volts(43)»<Resist+.4S)-Volt5(51  >*.45>/R2 

1252  Pou( 44 )-( Volts( 51  1-Volt 3(44  )  )•( Volts(44  )«( Resist* ,S5  >-Vol t s( 51  >*.55!/R2 
1281  Pou(2  )-< Volt s(46  )-Volts<  2  )  >*< Volts(2  >•( Resist +  .4 )-Volts(4B )• .4 )/R2 
1284  Poy( 5  )■( Vo  1 1  s ( 4G >-Vol ts( 5  )  >*< Vol t s( 5  >*< Resist* .2  )-Vol t s< 48 >*.2  )/R2 
1290  Pow<  4  )- ( Vol t  s( 4S  )-Voi t  s ( 4  ) >*  ( Vol ts(4 )*(Resist+. 1  )-Volts( 46  )• . 1  >/R2 

1294  Pow(9>«(Volta(47>-Volts(9>>*(Volts(9>»(Resist+.3>“Volts(47>*.3>/R2 

1295  Pou(25)-(Volts(49  )-Vol ts<  25  )>*<Volts(25)*<Resist+.3 )-Vol ts(  49  )•  .  3  )/R2 

1297  Pow(29)-<Volts(49  )-Vol t a ( 29))*(Volts(29)*(Resist+.2>-Volts(49)*.2  )/R2 

1296  Resist-1. 972 

1299  R2»Resist*Resist 

1300  Pou( 15  '-(Volt s(47  )— Vo 1 1  s ( 15  >  )*( Volts< 1 5 >* < Resis t+ . 2  )-Volts(47 >• .2  ),R2 

1301  Resist-1 .9184 

1302  R2-Resist»Resist 

1303  Pow( 14  >-( Vol t s( 47  )-Vol t s< 14>)*< Vo lts(14)*< Resist*. 1  )-Vol ts( 47  >• . 1  >/R2 

1304  Resist-1 .886 

1305  R2-Resist*Resi3t 

1306  Pow(4 >-(Volts(4B  )-Volts(4  ) )«( Volt s( 4  )*( Resist* . 1  )-Volts(46  >• .  1  >/R2 

1307  Pow< 6 >-( Volt s( 4B )-Vol ts( 6  )  >*( Volts( 6  )•< Resist*. 1  >-Volts( 4B >• . 1  )/R2 

1 308  Pow( 30  )-( Vol t  s( 49  )-Vol ts<  30>>*(Volts(30>*(Resist+.4  )-Vol t s( 49 )* . 4  )/R2 

1 309  Resist-! . 862 

1310  R2-Resist*Resist 

1311  Pow(27)-(Volts(49  )-Volts( 27  )  )•( Vol ts( 27  >•( Resist +.5  )-Volts( 49 )• .5  )/R2 
’312  Pow<  23 )-( Vol t s( 49  )-Vol t s<  28  )  )*( Volts< 28  >*( Resis t*0.  >-Volts< 49 >*0.  >/R2 

1313  Pow( 31  )-( Vol t s< 50  )-Vol ts( 31  >  )• ( Vo  1 1 s ( 3 1  )*(Resist+.S5  )-Volts(50  )*.B5  )/R2 

1314  Powf 33  )-( Volt s(50  >-Volts<  38  ) >*( Volts( 38  >•( Resi st* .3 )-Vol t s< 50 >• .3 )/R2 
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1315  Resist-1 .852 

1316  R2-Resist*Resiat 

1317  Pow<  24  )«<Volts(49)-Volts<24  ) >•< Vol t sv 24  >♦< Resi st  +  . 1  >-Volts< 49 >* .  l  )/R2 

1318  Pow<33)-<Volts(S0>-Volts<33)>«<Volts<33>*<Resist+.7>-Volts<50>*.7>/R2 

1319  Resist-1 .8333 

1320  R2-Resist#Resist 

1321  Pow< 1 1  >-< Volts(47)-Volts< 11  )  >•< Volts<  1  1  )•< Resist* . 5  >-Vol ts< 47 )* .5  )/R2 

1322  Pou( 12  >-( Vol ts<  47  >-Volts( 12  >>*< Volt s( 12  )•< Resist +. 1  )-Volts< 47 >• . i  '.'R; 

1323  Pou( l3>-(Volts(47)-Volts< 13))*(Volts( 1  3  >•< Resist* .2  )-Vol t s<  47  >«  .2  >/R2 

1325  Pouf  4 1  )-< Volts( 51  )-Volts<41  ))«(Volts(41  )«(Resist  +  .1 )-Volts(51  )•.)  )/R2 

1326  Pou( 45 >-(Volts(5l  )-Volts(45  >  >•<  VoltsUS  >*< Resist*. 6  )-Volts(S1 )».6  )/R2 

1327  Resist-1 .8182 

1328  R2-Resist*Resist 

1329  Pow(21  )-<  Volts<  48  >-Volts( 21  >  >*<  Volts<  21  >*( Resist* . 7  )-Vol t s< 48 >* . 7 >/R2 

1 330  Res i  st- 1 .8 

1331  R2-Resist*Resist 

1332  Pow<  8  >-( Vol t « ( 46  )-Volts(8  )  >* < Vol t s< 8  >♦ < Resis t+ . 2  >-Vol t s< 46 )«.2  )/R2 

1 333  Pow( 1 9  >■< Vol t s( 48  )-Vol t s< 19))«<Volta!19)«(Resist+,2 )-Vol ts<  48  >* . 2  >/R2 

1334  Poy(20)-(Volts(48)-Volts(20)  >*< Vol t s( 20 >*< Resi st* . 7 >-Vol t s( 43 >* . 7 >/R2 

1335  Poui( 22  >-< Volts<  48  >-Volts<22  )  )•<  Vol  ts(22  >*< Resist*.  2  >-Volts<  48  )».2  >/R2 

1336  Res l st- 1 . 6667 

1337  R2-Resist »Resi st 

1338  Pow< 1  >-< Vol ta( 46 >-Vol ts< 1  )  > * ( Vo 1 1 » < 1  >*( Resist +  .3 )-Vol ts(  48 )* .3 1/R2 

1 339  Pow<  3  )-< Volts< 46 >-Vol t a ( 3  )  )*( Volts <  3  )•( Res  1st +.7 >-Volts( 46 )• .7 )/R2 

1 340  Poui<7  )-<  Vol  t  s<  46  )-Vol ts(  7  ))»<Volts<7)*(Resist+.2  )-Volts( 46  >*.2  >/R2 

1341  Pou< 16  )-<  Voltsf  48  )-Volts( 16  )  )*( Volts< 16  )*< Resist* .4 )-Volts( 48 >».4 )/R2 

1 342  Pou( 17  )-<  Volt s(48  >-Volts<  17  )  )*( Volts< 17  )*( Resist*. 4  >-Volts(48  )« .4  )/R2 

1343  Pou<  26  >-<  Vol ts<  49  )-Volts(26  )  )*<  VoJts(26  )•( Resist +.2 )-Vol t s(  49 )• . 2  >/R2 

1344  FOR  I-l  TO  45 

1345  PRINT  'BLOCK  t  “»I,"  POWER  «  *»PowU),“  WATTS- 

1346  WAIT  .5 

1347  NEXT  I 

1348  Ptot-0 

1348  FOR  1-1  TO  45  "" 

1350  Ptot«Ptot+Pow< I  ) 

1351  NEXT  I 

1352  Pavg«Ptot/45 

1353  PRINT  -  " 

1354  PRINT  “AVERAGE  POWER  -  “iPavg,"  WATTS- 

1355  PRINT  ”  " 

1356  FOR  1-1  TO  45 

1 357  Perdi f-AB5( ( Pavg-Pow< I )  )/Pavg  )« 1 00 

1358  PRINT  “HEATER  t"jl,*  PERCENT  DEVIATION  “iPerdif 

1359  WAIT  .5 

1360  NEXT  I 

1361  END 
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APPENDIX  F 

UNIFORM  DATA  REDUCTION  PROGRAM 

10  I  II  !  I  I  I  !!  I  I  !  I  M  II  M  I  I  I  M  I  I  I  !  I  I  I  I  II  M  !  M  I  I  I  I  I  I  I  II  I  I  I  II  II  I  M  II  II  I  I  I  I  I  I  I  I  I 

20  H  TEMPERATURE  ACQUISITION  AND  REDUCTION  PROGRAM  n 

30  !'  <  STEADY  STATE )  n 

31  !'  PROGRAM  PERFORMS  THE  REQUIRED  FUNCTIONS  NEEDED  TO  OBTAIN  DATA  u 

32  ! !  FROM  THE  ACQUISITION  SYSTEM  AND  CONVERT  THE  THERMOCOUPLE  OUT-  i ! 

33  1 !  PUT  VOLTAGES  INTO  TEMPERATURE  VALUES.  THE  TEMPERATURES  ARE  n 

34  ! i  PRINTED  OUT  AND  THEN  REDUCED  TO  COMPARE  WITH  THEORETICAL  VALUES  ' ' 

40  I  l  i  i  l  l  i I  i  l  I  I  ll  l  l  i  |  i  l  i  j  i  |  l  i  i  i  i  M  I  l I i |  I  | | ( j I  i  i  i  i  i  ii  ii  i  i i  |  i  i i  i i I  |  i  i  l  |  i  i  i  i  i 

50  REAL  Volts<60) 

60  REAL  T( 60  ) 

65  REAL  Ndt(4S),Ndt1(45),Lnu<45)  ,Lgr( 45 ) 

66  REAL  Power ,Tamb  ,G  ,Kpq  ,Delx  , As  .Per ,Ta .Qconv ,Lbar 

67  REAL  Qf lux .Del t  ,Tf i In  ,Bet a  .Spvol  .Dynvi sc .Qcond 
66  REAL  K i nv l sc  ,K f  ,Pr  .The t a  ,Hx  ,Grm ,Hux 

6S  PRINTER  IS  CRT 

70  PRINT  "ENTER  OUTPUT  DATA  FILE  NrME " 

71  INPUT  Filename* 

72  PRINT 

73  PRINT  "ENTER  POWER  VALUE  IN  WATTS" 

74  INPUT  Power 

75  PRINT 

76  PRINT  "ENTER  SHROUD  SPACING  IN  mm" 

77  INPUT  SpaceS 

76  (PRINT 

79  'CREATE  BDAT  * THRED I  5 :  .700 , 1 " .45 

80  'ASSIGN  0Path3  TO  " THREO 1 5 :  .700 , 1" 

81  'PRINT  DATESf  TIMEOATE  ) 

e2  'PRINT  Power."  WATT  INPUT  POWER" 

84  'PRINT  "  45  STRIPS  NO  SHROUD" 

65  'PRINT  “  BDAT  FILE  THRED I  5  :  .700  .  t " 

66  'PRINT  “  “ 

89  I  I  i  i  |  i  I  i  i  i  i  i  I  i  ii  i  ii  i  i  i  m  i  i  i  i  SURFACE  f 1-20 

90  OUTPUT  709 ( "CONFMEAS  OCV  .  100- 11 9, USE  0" 

91  FOR  1-1  TO  20 

92  ENTER  709;Volts<  I  ) 

93  T (  I )• . 0006797  +  1 25825. 1 328* Vo  1 1 s < I > )- i 607789 . 2467* i Vol t s (  I >«Volts( I ) ) )-(  219 

52034. 3364.  (Vo  It s< I  )"3)  >+<8370810996. 1 874 • ( Vol t s (I >"4  >  >+ . 039 

94  NEXT  I 

95  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  (  i  ii  i  i  ii  i  i  SURFACE  *2 1 -140 

96  OUTPUT  709; "CONFMEAS  DCV  .200-2 1 9  .USE  0" 

97  FOR  1-21  TO  40 

98  ENTER  709;Volts<  I  ) 

99  T ( I  >-.0006797  +  <  25825. 1328* Vol ts< I  >  >-< 607789. 2467* ( Vol t s (  I >*Vclts(  I >  >  >-(219 
52034 . 3364 •(  Vol t s(  I  >*3  >  )  +  <  8370810996. 1 874»< Vol t s< I >* 4 ) )+ . 039 

100  NEXT  I 

101  i i i i i ii . .  i  i  i  i  i  i  i  SURFAC  141-445 

102  OUTPUT  709; "CONFMEAS  OCV  .300-304  .USE  0" 

103  FOR  1-41  TO  45 

104  ENTER  709 ; Voi t s<  I  > 

105  T (  I  >-.0006797  +  (  25825. 1 32S*Volts< I > >-( 607789. 2467* ( Vol ts< I >*Vol tsl I  )  > >-( 21 S 

52034. 3364.  (  Volt s(  I )"3) >+<8370810996. !874«(Volts< I >"4  >  >+.039 

106  NEXT  I 


107  M!  PRINT  OUTPUT  TEMPERATURES 

108  IPRINT  "  ‘ 

109  'FOR  1-1  TO  45 

110  ! PRINT  "HEATER  *"|I."  TEMPERATURE  IN  DE6.  C  - “ ; T< I  ) 

1 1 1  INEXT  I 

112  ! OUTPUT  SPath3 1 T( •  ) 

115  'PRINT  ’  “ 

IIS  IPRINT  ’  BATH  TEMPERATURES  (TOP  TO  BOTTOM)" 

117  IPRINT  "  ' 

118  OUTPUT  709 : "CONFMEAS  DCU .317-31 9 .USE  0" 

119  FOR  1-57  TO  59 

120  ENTER  709iUolts< I  ) 

121  T( I  )-.0006797+< 25825. 1328»Uolts< I >  >-(507789. 2467* (Uolts< I >*Uolts( I >  >  >-(219 
52034 . 3364* ( Volts! I  ) ‘3 >)+< 837081 0996 . 1 874*< Uol ts< I >”4)  >+.036 

122  (PRINT  "TEMP.  DES.  C  "|T(I> 

124  NEXT  I 

126  T <  22  )-25 

127  T ( 29  )-25 

133  'PRINT  "INPUT  AMBIENT  TEMPERATURE" 

134  Tamb«( T (57)+T(58)+T(59))/3 

135  PRINTER  IS  701 

136  PRINT  "  ” 

137  PRINT  DATES! TIMEDATE) 

138  PRINT  “DATA  REDUCTION  OUTPUT- 

139  PRINT  "  " 

140  PRINT  “POWER  (WATTS)  :"  .Power 

141  PRINT 

142  PRINT  " SHF.ul -  POSITION  (mm)  Spaces 

143  PRINT 

144  PRINT  " AMEIENT  TEMPERATURE  (C)  : ”  .Tanb 

145  PRINT 

146  PRINT 

147 

148  'ASSIGN  @Path  TO  FilenaMeS 

14S  'ENTER  @Pa t h i  T(  •  ) 

150  ' 

151  6-9.81 

1  52  Kpg- .1421 

153  Delx-. 006731 

154  As- . 0239* . 0078 

155  Per-2«. 0239+2*. 0078 

158  Lbar-As/Per 

157  AS-''Qcond" 

158  BS«“X/L" 

159  CS* “NOT  ” 

160  DS-"LOG(Nux >" 

161  ES- “LOG<  Grx*  ) " 

162  FS-" DELT" 

163  6$-"Hx“ 

164  HS-'Nux” 

165  IS-"Grx*" 

166  J$-”TeMp" 

168  PRINTER  IS  CRT 
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169  t  !  }  I  I  |  { | | j | | { !  |  ! I |  I  |  |  !  |  i  |  i  {  M  m  M  i  i  {  i  I  (  !  |  M  |  ( l f | 

170  !  ENTER  CONDUCTION  LOSSES  ! 

171  ( | | I i li  I  II  I  1  I  f  I  i u  M  i  M  I  |  I  !  |  |  !  j  !  |  II  |  )  j  |  |  I  |  M  I  !  i  | 

172  PRINT  “INPUT  CONDUCTION  LOSSES  AS  COMPUTED  FROM  ELLPACK” 

173  INPUT  Qcond 

174  I  I | | i | | | | li | l l l l | | |  i  li  M  l  l  |l  |  l  |  l  l  M  l  l  l  l  |  ll  I  I  l  M  I 

175  CREATE  BOAT  Filenames  ,45 

176  ASSIGN  ®Path3  TO  Filenames 

177  PRINTER  IS  701 

178  FOR  1-1  TO  45 

179  PRINT  “  " 

180  PRINT  “  ” 

181  PRINT  ’  SURFACE  #“  .1 


182  PRINT  “  " 

183  PRINT  USING  "2X .AAAAA ,3X . AAA ,3X ,AAA ,3X  .AAAAAAAA  ,2X  .AAAAAAAAA  ,3X  ,AAAA  ,4X  ,AA 
,4X  .AAA  ,3X  ,AAAA  ,2X  ,AAAA" i AS ,B$ ,C$ ,D$  ,E$ ,F$ ,G$  ,HS  ,IS , J* 

184  Qconv-Pouer-Qcond 

185  Of lux-Qconv/As 

186  Del t-T ( I )-Tamb 

187  Tf i lm-273. 1 S+( T< I  1+Tamb  >/2 

188  Beta»( -7 . 944858E-8*  <  Tf i lm “2  )  )+< 5 . 7356479E-5»Tf i lm )-. 00978 1 0563 

189  Spvol «( 4 . 69699E-9* <  Tf l lm'2  ) )+< -2 .S374SE-6»Tf i lm 1+ . 00 1 34 1 903 

1 90  Dynvisc-( 3.234851 1E-7*(Tfi lm“ 2 1 )  +  <-. 0002 1 474487*Tf l lm  1+ . 036 1 66792 

191  Kinvisc-Spvol*Dynvisc 

192  Kf«< 1 .41818181 8E~3*T film  )+. 1866 

193  Pr»( 4 , 65706E-3* (  Tf i lm*2  1  >+ < -2 . 922094*Tf i lm  1+463.3319 

194  Theta!  =>Power7(  .0220*Kf  ) 

195  Ndt 1 ( I  )»De it /Theta  1 

196  Nd-Ndt 1 < I  ) 

197  Temp-T(I) 

199  Theta-Qf  lux»Lbar-/K  f 

200  Ndt ( I  )-Oelt/Theta 

201  H.x-Qf  lux/Delt 

202  6rm«(G*Beta»0flux«(Lbar)‘4)/(Kf»Kinvisc“2  ) 

203  Nux«Hx »Lbar/K f 

204  Lgr< I  )-LGT< Grm ) 

205  Lnu( I  )-LGT( Nux  1 

206  Lq-Lqi-U) 

207  Lnn-Lnu( I  ) 

209  Xtol-( 3 1 . 75- (6.985+ ( 1-16  >*1 .270)1/31.75 

210  PRINT  USING  “IX ,0 . DODD ,2X ,D . ODD , 1 X ,D . DDD  ,2X  ,D . DDD  ,5X  ,D . DDD  , 4X  ,DD . ODD  ,2X  ,DD 
DD.D  ,1X  .DD.DDD  .DDDOD.D  ,1X  .DD.DD"  i Qcond  ,XtoI  ,Nd  ,Lnn  ,Lg  , Celt  ,H.v  ,Nux  , Grm  .Temp 

215  NEXT  I 

225  OUTPUT  «Path3sT< •  1  ,Ndt ! ( •  1  .Ndt ( *  )  ,Lnu< » >  .Lgr( ♦  1 

226  PRINTER  IS  CRT 
235  END 
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APPENDIX  G 

NONUNIFORM  POWER  PROGRAM 


2 

4 

i ! !  POWER 

1  1  t  1  II  II  1  1  1  1  M  M  1  1  1  1 

INPUT  PROGRAM 

i i i t  m ii i i i n i i i ii i i i i  n  i  ii  i 

5 

•OUTPUT  709; "LOCAL  CLR“ 

7 

OUTPUT  709; “RST" 

10 

REAL  Vol t  s( 5 1  ) 

20 

REAL  Pow( 51  ) 

30 

REAL  Ptot  ,PavQ  .Perdi f 

920 

OUTPUT  709; “CONFMEAS 

DCV  .308-315  , USE  0" 

930 

FOR  1*1  TO  8 

931 

J-9-I 

940 

ENTER  709;Volts( J  ) 

950 

PRINT  ”  “ 

SE2 

NEXT  I 

953 

OUTPUT  709; "RST" 

954 

OUTPUT  709; "CONFMEAS 

0CV.316.USE  0" 

955 

ENTER  709;Volts<  46  ) 

956 

OUTPUT  709; "RST" 

971 

OUTPUT  709; "CONFMEAS 

OCV  ,403-4t7  .USE  0" 

972 

FOR  1=9  TO  23 

973 

ENTER  709 ; Uol t s ( I  ) 

977 

NEXT  I 

973 

OUTPUT  709; “RST" 

979 

OUTPUT  70S; "CONFMEAS 

DCV  .400-401  ,USE  0" 

980 

FOR  1*47  TO  48 

951 

ENTER  709;Volts< I  ) 

962 

NEXT  I 

9B3 

OUTPUT  709; "RST" 

995 

OUTPUT  709; "CONFMEAS 

DCV  .503-517, USE  0" 

966 

FOR  1-24  TO  38 

967 

ENTER  709; Volts< I ) 

969 

NEXT  I 

990 

OUTPUT  709; "RST" 

991 

OUTPUT  709; "CONFMEAS 

DCV  ,500-501  .USE  0" 

992 

FOR  1-49  TO  50 

993 

ENTER  709 ;  Vol t  s(  I  ) 

994 

NEXT  I 

996 

OUTPUT  709; "RST" 

997 

OUTPUT  709; "CONFMEAS 

DCV  .602-608  .USE  0" 

998 

FOR  1-39  TO  45 

999 

ENTER  709;Volts<  I  ) 

1000 

NEXT  I 

1001 

OUTPUT  709; "RST" 

1002 

OUTPUT  709; "CONFMEAS 

DCV, 600  .USE  0" 

1003 

ENTER  709; Volts! 51 ) 

1004 

OUTPUT  709; "RST" 

1005 

PRINT  "  " 

1006 

PRINT  '  " 

1050 

PRINT  "INPUT  VOLTAGE 

D.C.  VOLTS  PAN  1 " ;Volts(46  > 

1051 

PRINT  "INPUT  VOLTAGE 

D.C.  VOLTS  PAN  2";VoltsU7> 

1052 

PRINT  "INPUT  VOLTAGE 

D.C.  VOLTS  PAN  3“ ; Vol t s< 48 > 

]  053 

PRINT  "INPUT  VOLTAGE 

D.C.  VOLTS  PAN  4" ;Volts(49  ) 
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1054  PRINT  "INPUT  U0LTA6E  D.C.  UOLTS  PAN  5"iUolts(S0> 

1055  PRINT  "INPUT  UOLTAGE  D.C.  UOLTS  PAN  6"iUolts(5l  > 

1060  PRINT  "  “ 

1070  PRINT  "  “ 

1080  Res ist«2 . 0 
1109  R2*Resist»Resist 

1119  Pout  2  )-(Uolts(46  )-Uol t  s( 2  )  )« <  Uol t st  2 )•< Resist* .4 )-Uol ts( 46 )• . 4  >/R2 
1149  Powt  5  )»( Uol t st  46 >-Uolts<5  >  )* ( Uol t s ( 5 )* ( Res 2 st  + . 2 >-Uol t s ( 46 )• .2  >/R2 
1 179  Pow<  9  )*( Uol t  s ( 47  )-Uolts( 9  )  >•< Uoltst 9 >« ( Res ist+ . 3  >-Uol t s(  47 >* . 3  >/R2 
1 189  P ou ( 1 0  )*( Uol t  st  47  )-Uol t  s ( 10)  )*<Uolts( 1 0  )• ( Res l st + . 45 >-Uol t s < 47 >• . 45 1/R2 
1 194  Pow ( 18  >-< Uol t s(48  )-Uol t s < 1 8  >  )*( Uol t s< 1 8  )• ( Resi s t+ . 55  )-Uolts(48  >*.55 1/R2 

1 197  Pout  23  )-< Uol t s(48 >-Uolts<  23  )  >*< Uol tst 23 )*( Resi st*  .55 )-Uolts<48 ) • .55 )/R2 

1 198  Pou<24 >«( Uoltst 49  )-Uol tst 24 >)»( Uol tst 24 >•< Resist*. 1  >-Uolts< 48 >• . 1  >/R2 

1199  Pow(25)*(Uolts(49  >-Uolts<  25  )  >* ( Uol t s ( 25 >•< Resi st* . 3 >-Uol tst 48  >* - 3  )/R2 
1238  Powt 29  )-< Uoltst 49  )-Uol tst 29)  >•< Uoltst 29 >• ( Res i s t + . 45 )-Uol ts( 49 >• . 45 >/R2 
1241  Foul  32  )«<  Uol t st  50  )-Uol t  st  32  )  )* ( Uol tst  32 >•< Resi st* .6  )-Uol tst 50 >* . 6  >/R2 

1243  Pow( 34 >»<  Uol t  s <  50  ) -Uol tst 34  )  >•< Uoltst 34  )*( Res i st+ . 3 )-Uol t s< 50 )* . 3 >/R2 

1244  Pout  35  >■( Uol tst  50  )-Uol ts( 35  )  )♦< Uoltst 35  >*( Resist* .6  >-Uolts(  50 )» .6 >/R2 

1245  Pout  36 >*( Uol t s(50 >-Uolts< 36  )  )*( Uol tst 36  >»(Resi5t  +  .6  >-Uolts( 50 >* .6  >/R2 

1246  Pout  37  )«< Uol tst 50  >-Uolts(37  >  )• ( Uol t st 37 >« < Res l s t + . 6 )-Uolts(  50 )• .6 >/R2 

1247  Powt  39  >*( Uo 1 1  st  5 1  )-Uol t s< 39  )  >• ( Uol t s ( 39  )• ( Res l s t  + . 3  )-Uol t s< 5 1  >*.3>/R2 

1249  Pout  40  )*( Uol t  s ( 5 1  )-Uolts(40)  >*< Uol tst 40 >*( Resist*. 3 )-Uel tst 51  )*.3)/R2 

1 250  Pout  42  >-<  Uoltat 51  l-Uoltst 42  )  )*<Uolts( 42  >•< Resist* .55  )-Uolts( 51  ) • . 55  )/R2 

1251  Pou(43  )-(Uolts(5l  )-Uolts(43)  >*< Uol tst 43 )*(  Resist* .45 >-Vol t st 51 )*.45)/R2 

1252  Pout 44  )■( Uol t s ( 51  )-Uol t s! 44  )  >*< Uo 1 1 st 44 )*( Res i s t  + .55 )-Uol t s ( 5 1  >*.55)/R2 
1281  Pout  2  )-<  Uoltst  46  )-Uolts<  2  )  )•( Uol t s ( 2 >* ( Res l st + . 4  >-Uol t s ( 46 )• .4  )/R2 

1 284  Pout  5  )-<  Uoltst  46  )-Uol tst  5  )  )•( Uoltst S >« ( Res ist + . 2 )-Uol t st  46 >• .2 )/R2 
1290  Pou(4)«(Uolts<46)-Uolts(4))»(Uoltst4)«(Resist+.l  )-Uol t st  46  )• . 1  ) / R  2 

1294  Pow(9)“(Uolts(47)-Uoltst9)  )*t Uoltst 9  )*( Resist+ .3  )-Uolts( 47  )* .3  )/R2 

1295  Pout  25  )* ( Uol t s ( 49  )-Uol t st  25  )  )*  <  Uol t st  25  )•( Resi s t  + . 3 )-Uol t s ( 49  >• . 3 )/R2 
1 297  Pout  29  )-t Uoltst  49  )-Uol tst  29  ))•( Uoltst  29  >•< Res i st+ . 2  )-Uo 1 1 s ( 49  )*.Z  >/R2 
I29e  Resist* 1 . 972 

1299  R2*Res l s t »Res i s t 

1300  Pout  15  >“( Uoltst  47  )  — Uoltst  15  )  )•( Uoltst 1 5 >* t Resi s t + . 2  )-Uol t s t 47  )• . 2  )/R2 

1 30 1  Res i st *  1 . 9 1 84 

1302  R2*Res l s t *Res i s t 

1503  Pout  1  4  >•*<  Uoltst  47  >  —  Uoltst  1  4  )  >♦  (  Uoltst  14)*iResist+.1  >— Uol  tst  47  >•  .  1  )/RI 
1 204  Resi s  t » 1 .886 

1305  R2*Res l s t *Res l s t 

1 306  Pout  4  )»( Uol t si  46  )-Uclt st  4  )  >•<  Uoltst  4  )•( Resist+ . 1  )-Uol tst  46  )• . 1  )/R2 
1207  Pout 6)*(Uolts<46)-Uolts<6))«(Uolts<6  )*<Resist  +  . 1  )-Uol t s ( 46  )« . 1  >/R2 

1303  Pout30)"(Uolts(49)-Uolt5(30;)»(Uolt5(30)«(Resist+.4)-Uolts(49)#.4)/R2 

1309  Res l s t  *  1 . 862 

1310  R2*Res i s t »Res i s t 

1311  Pout  27  )-( Uol tst  49  >-Uc 1 1 s ( 27  )  )•< Uoltst  27  )*( Resi st+ .5  >-Uol t st  49  >« . 5 )/R2 

1312  Pout  28 )-(  Uoltst 49  )-Uol tst 28  )  )• < Uol t st 28 >• ( Res l s t +0 . >-Uoltsi 49  1*0.  )/R2 

1313  Pou<31  >*( Uol tst  50  >— Uoltst  31  )  > « <  Uo 1 1  a  <  31  )*( Resi st+ . 65  )-Uol t st 50  )• . 65  )/R2 

1314  Pout  38 )*<  Uol tst  50  )-Uol t s ( 38))«(Uolts(38 >• ( Res l s t  +  . 5  >-Uol t s ( 50 >• .3  )/R2 

1315  Resist  *  1 . 852 

1316  R2*Res l s t »Res i s t 

1317  Pout  24  >-< Uoltst 49  )-Uo 1 1 s < 24  )  )•( Uoltst  24  )•( Resi st+ . 1  )-Uol tst  49  >• . 1  >/R2 

1318  Pout  33  )•( Uol tst  50  )-Vol t  s< 33)>*t Uoltst 33  )*(Resist  +  . 7)-Uo ltst 50  )*.7)/R2 

1319  Resist-1 .8333 

t  jZ  0  R2*Res i s t «Res  i  s t 

1 321  Pout  II  >-<  Uol tst  47  )-Uol tst  1 1  )  )•( Uol tst 1!  )*( Resist +  .5  )-Uol tst  47  )• .5  )/R2 

1322  Pout  12  )«< Uol tst  47 >-Uol t st 12  )>•<  Uoltst 12 )* (Resist* . 1  )-Uol tst  47  >•  .  1  )/R2 

1 323  Pout  1 3 )»<  Uoltst  47 )-Uol tst  1 3  )  >*( Uol tst 1 3 )*( Resist* .2 )-Uolts(  47 !• .2 )/R2 

1325  Powt  4 1  )■( Uol t st  5 1  )-Uo 1 1 s ( 4 1  )  ) • t  Uol t s ( 4 1  )•( Resist*. I  )-Uo 1 1 s <  5 1  > • . 1  )/R2 

1  326  Pout  45 )-( Uo 1 tst 51  )-Uolts( 45  >  >•<  Uol tst  45 )*< Resist*. 6 >-Uolts<  51  )».6  >/R2 
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1327  Resist-! .8182 

1328  R2-Resist*Resist 

1329  Pouf  2 1  )-( Volts(48>-Volts<21  >)*<Vo2ts(2!  >* ( Resist +. 7 >-Vc2 ts<48 >* .7 >/R2 

1330  Resist-1 .8 

1331  R2-Resist«Resist 

1332  Pow(8 >-(Volts< 4g )-Volts<  8  >  )*< Volts<  8 )*( Resist +  .2  )-VoIts( 46  )*.2  >/R2 

1333  Pou( 1 9 )-< Vo2ts( 48 >-VoIts< 1 9  >  >«< Volts( 1 9  >•< Resi st+ .2 )-Vo2ts<48  )*  .2  )/R2 

1334  Pow( 20  >-< Volts( 48  )-Vo 1 1 s < 20  )  )*< Vol ts( 20  )*< Resist* .7  )-Volts< 48  )*.7  )/R2 

1335  Pow< 22  >«<  Voltst 48 )-VoIts<  22 )  )*<  Vol ts( 22  >«<  Resist+ .2  )-Volts< 48  )* .2  )/R2 

1336  Resist-1 .6867 

1337  R2-Resist*Resist 

1338  Pow( 1  )-( Vo2ts< 46  >-Volts<  1 ) >•< Volts( 1 )• ( Resist* .3  >-Vol ts< 46 )• .3  >/R2 

1339  Pow( 3  )-< Volts<  46  >-Volts<  3 )>*<Volts<3>*< Resist^. 7 >-Volt s( 46 )* .7  )/R2 

1 340  Pow<  7  )-<  Volts(46  )-Volts<7 ) >*< Volts( 7  )♦( Resist+ . 2 >-Volts( 46 >* .2  >/R2 

1341  Pow< 1 6  )-( Vol t  s< 48  )-Volts< 16))*<Volts(16  >• (Resist*. 4  >-Volts<  48 )* .4 )/R2 

1 342  Pow< 17)-(Volts<48  )-Volts< 1 7  )  )♦< Volts! 17  >*( Resist +  .4  >-Vol t a C 48 >* . 4  >/R2 

1343  Pow( 26  )•< Vol ts( 49)-Volts( 26  )  )*< Volts<26 >*<Resist+.2  >-Volts<  49  )*  .2  >/R2 

1344  FOR  1-1  TO  45 

1346  PRINT  “BLOCK  t  "si,"  POWER  -  *»Pou<I),*  WATTS” 

1347  WAIT  .5 

1348  NEXT  I 

1350  Piot-0 

1351  FOR  1-1  TO  15 

1352  Ptot-Ptot+Pou( I > 

1353  NEXT  I 

1354  FOR  1-31  TO  45 

1355  Ptot»Ptot+Pow( I  ) 

1356  NEXT  I 

1357  Pavgsides-Ptot/30 
1353  PRINT  *  " 

1359  PRINT  "AVERAGE  POWER  SIOES-  " iPavgsides  , "  WATTS" 

1360  PRINT  "  " 

1361  Ptot-0 

1362  FOR  1-16  TO  30 

1363  Ptot»Ftot+Pow( I  ) 

1364  NEXT  I 

1365  Pa vgcenter«Pt ot / 1 5 

1357  PRINT  "AVERAGE  POWER  CENTER-  Pevgcenter  , "  WATT5" 

1366  PRINT  "  ” 

1370  PRINT  “  " 

1371  FOR  1-1  TO  15 

1372  Perdi f- ABS( < Pavgsides-Pou( I )  >/Pavgsi des )» 1 00 

1373  PRINT  "HEATER  •* il,“  PERCENT  DEVIATION  "iPerdif 

1374  WAIT  .5 

1375  NEXT  I 

1376  FOR  1-16  TO  30 

1377  Perdi f«ABS< ( Pavgcenter-Pou< I ) ) /Pevgcenter  >*100 

1378  PRINT  "HEATER  #”«I,“  PERCENT  DEVIATION  "iPerdif 

1379  WAIT  .5 

1380  NEXT  I 

1381  FOR  1-31  TO  45 

1382  Perdi f»ABS< ( Pavgsides-Pou< I >>/Pavgsides  >*100 

1383  PRINT  "HEATER  fsl,”  PERCENT  DEVIATION  "iPerdif 

1384  WAIT  .5 

1385  NEXT  I 

1386  ENO 
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APPENDIX  H 

NONUNIFORM  DATA  REDUCTION  PROGRAM 


10  Ml  I  !  H  !  H  H  !!!!  H  !!!!!!!!!!  H  !!!!  1  !!  !  H  !!!!!!!  i!  !  !!  !!  H  I  i  |  j  t  |  MM!!  M  i 

20  ! !  TEMPERATURE  ACQUISITION  AND  REDUCTION  PR06RAM  ! ! 

30  M  (STEADY  STATE)  !! 

31  M  PROGRAM  PERFORMS  THE  REQUIREO  FUNCTIONS  NEEDED  TO  OBTAIN  DATA  I! 

32  ! !  FROM  THE  ACQUISITION  SYSTEM  AND  CONVERT  THE  THERMOCOUPLE  OUT-  ! ! 

33  ! !  PUT  VOLTAGES  INTO  TEMPERATURE  VALUES.  THE  TEMPERATURES  ARE  !! 

34  ! !  PRINTED  OUT  AND  THEN  REDUCED  TO  COMPARE  WITH  THEORETICAL  VALUES  ! I 

40  I  j  t  I  |  I  I  I  j  I  I  I  |  I  I i  | i | j  j | I j i I (|  i i i I j ! ( (  I  l  I  |  !  j  |  j  M  j  (  I  !  M  j  !  I  M  u  (  I  !  I  I  i  I  l  1  i  M 

50  REAL  Volts<60> 

G0  REAL  T<  60  ) 

G5  REAL  Ndt ( 45  )  ,Ndt 1 ( 45  )  ,Lnu ( 45  > ,Lqp(  45 ) 

GE  REAL  Power  1  ,Power2  .Tarnb  ,G  ,Kpo  ,0eix , As  .Per  ,Ta .Gconv .Lbar 
B7  REAL  Qf lux ,Delt  .Tfilm  .Beta  .Spvol .Dynvisc .Qcond 
G8  REAL  Kinvisc  ,Kf  ,Pr  .Theta ,Hx ,6rm ,Hux 
E9  PRINTER  IS  CRT 

70  PRINT  "ENTER  OUTPUT  DATA  FILE  NAME- 

71  INPUT  Filenames 

72  PRINT 

73  PRINT  "ENTER  TWO  POWER  VALUE  IN  WATTS" 

74  INPUT  Power!  .PowerZ 

75  PRINT 

7E  PRINT  "ENTER  SHROUD  SPACING  IN  mm" 

77  INPUT  Space* 

78  IPRINT 

79  ICREATE  BOAT  "THRED1 5 :  .700 . 1” ,45 

80  i ASSIGN  8Path3  TO  “ THRED 1 5 :  ,700 . 1 " 

8!  IPRINT  DATESUIMEOATE  ) 

83  IPRINT  Power."  WATT  INPUT  POWER" 

84  IPRINT  *  45  STRIP5  NO  SHROUD" 

85  'PRINT  "  BOAT  FILE  THRED 1 5 :  ,700 , 1 " 

86  IPRINT  '  " 

89  I  I ! I  I  I !! I  I ! II  I  II  i  M  I  I  I  i  I  I  l  i  i  SURFACE  *  1 -20 

90  OUTPUT  709; "CONFMEAS  DCV  .  100-1 19. USE  0” 

91  FOR  I-l  TO  20 

92  ENTER  709;Volts< I > 

93  T ( I >-.000G797+<  25825 . 1 3Z8*Vol t s ( I  ) )-( E07789 . 2467*( Voi ts( I >*Volts< I ) ) )-(  2 1 S 
52034. 33G4*( Volt s( I ) '3 > )  +  < 83708 1 0995 . 1 874* ( Volt sC I >"4 ) >+.039 

94  NEXT  I 

95  HI  i  M  I  !  I  I  I  M  I  I  I  !  I  !  1  I  i  I  I  I  I  !  SURFACE  12 1  -140 

96  OUTPUT  709; "CONFMEAS  DCV  .200-2 1 9 .USE  0” 

97  FOR  1-21  TO  40 

98  ENTER  709 ; Vol t a( I  ) 

99  T( I >-.0006797+< 25825. 1328»Volts< I >>-<607789. 24G7*< Vol ts(  I )«Volts<  I >>>-(219 
52034. 3364* (Volts'  I  )"3  )  >  +  (8370810996 . 1 874*< Volt5( I  )"4  > >+.039 

100  NEXT  I 

101  I ! M  M  »  « I i I i M  I  M  I  M  i  I  I  I  !  1  SURFAC  #4 1 -*4S 


102  OUTPUT  709; "CONFMEAS  DCV  .300-304 .USE  0" 

103  FOR  1-41  TO  45 

104  ENTER  709;Volts< I > 

105  T ( I )-.0006797+( 25825. 132 8* Volt a( I >  >-( 607789 . 2467* ( Vol ts(  I >*Voits< I ) >>-(219 
52034. 3364* ( Volt s( I >*3  >  >+< 837081 0996 . 1 874* ( Vol ts< I >'4 ) >+.039 

106  NEXT  I 
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107  !!!  PRINT  OUTPUT  TEMPERATURES 

108  'PRINT  '  * 

109  'FOR  1-1  TO  45 

110  'PRINT  ’HEATER  *"»I.*  TEMPERATURE  IN  DEG.  C  -“;T(I) 

111  (NEXT  I 

1 12  'OUTPUT  @Path3»T( *  ) 

115  ! PRINT  ■  " 

116  'PRINT  “  BATH  TEMPERATURES  (TOP  TO  BOTTOM >" 

117  .'PRINT  *  “ 

118  OUTPUT  709 1 “CONFMEAS  DCU  ,31 7-31 9  .USE  0“ 

119  FOR  1-57  TO  59 

120  ENTER  709;Uolts< I ) 

121  T( I >-. 0006797+ ( 25825. 1328*Uolts< I >  >-< 607789. 2467*(Uolts( I  )*Uolts( I )  >  )-<  21 9 
52034. 3364M Uolts( I  K3 >)+< 837081 0996, 1874«< Uolts< I >‘4  )  )+.036 

122  (PRINT  "TEMP.  DEG.  C  “  * T < I > 

124  NEXT  I 

125  T ( 22  )-25 

126  T ( 29  )-25 

133  (PRINT  "INPUT  AMBIENT  TEMPERATURE" 

134  Tamb»( T( 57  )+T(58)+T(59))/3 

135  PRINTER  IS  701 

1 36  PRINT  "  " 

137  PRINT  DATES' TIMEDATE  ) 

138  PRINT  "DATA  REDUCTION  OUTPUT- 

139  PRINT  "  " 

140  PRINT  'POWER  LEUELS  (WATTS)  .Power  1  ,Pouier2 

141  PRINT 

142  PRINT  “SHROUD  POSITION  (m)  Space* 

143  PRINT 

144  PRINT  "AMBIENT  TEMPERATURE  (C)  :",Tamb 

145  PRINT 

146  PRINT 

147  ' 

148  'ASSIGN  @Path  TO  Filenames 

149  ! ENTER  @Path  *  T( *  > • 

150  ' 

151  G-9.81 

152  Kpg- .1421 

153  Del*;-. 006731 

154  As- . 0239* . 0078 

155  Per=2* .0239+2* .0078 

156  Lbar— As/Per 

157  A$="Qcond“ 

159  B$-"X/L" 

159  CS--NDT" 

160  DS-"LOG<  Nux  )  “ 

161  ES«"LOG(G rx*)‘ 

162  FS--DELT" 

163  6S-"Hx  “ 

164  H$» “Nux ” 

165  IS-’Grx*" 

166  JS-"Temp" 

168  PRINTER  IS  CRT 

1 5g  i  |  i  i  m  j  |  i  i  m  i  i  i  M  i  !  m  i  I  |  M  i  (  (  M  |  ii  m  i  |  (  M  i  i  |  (  i  i  i 

170  '  ENTER  CONDUCTION  LOSSES  ( 

171  II  i i i l I  I  II  I  I  II  I  II  i  |l  I  l|  l  |  II  i  i  i  i  i  l  I  I  II  i  i  i  li  i  |  |  |  l  i 

172  PRINT  “INPUT  SIDE  COLUMN  CONDUCTION  LOSSES  AS  COMPUTED  FROM  ELLPACK" 

173  INPUT  Qcond 

174  | i j! | i | i i i | | i i  i  i  i  I  I  i  |  |  ii  I  I  M  1  M  i  I  i  I  M  |  i  I  i  l  |  I  |  !  I  I 

175  CREATE  BOAT  Fi ienameS .45 

176  ASSIGN  0Path3  TO  Filenames 

177  PRINTER  IS  701 
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SURFACE  t\I 


178  FOR  1-1  TO  15 

179  PRINT  *  • 

180  PRINT  “  * 

181  PRINT  ‘ 

182  PRINT  '  ” 

183  PRINT  USINS  "2X  .AAlffl^WC  ,AAA ,3X ,AAA ,3X  .AAAAAAAA ,2X .AAAAAAAAA ,3X ,AAAA ,4X ,AA 
.4X  ,AAA  ,3X  ,AAAA  ,2X ,AAAA" »A$ ,B$ ,C* ,D* ,E$ ,F$ ,SS ,HS  ,1$ ,JS 

184  Qconv-Powerl -Qcond 

185  Qf lux-Qconv/As 
188  Deit-T< I )-Tanb 

187  Tfi in-273. 1 5+<  T< I )+Tamb )/2 

188  Beta-( -7. 944858E-8*<  T f iln"2  ) )+( 5 . 73S64?9E-5*Tf i In >-.0097810583 

189  Spvol-<4.69699E-9«(Tfilm*2  )  >+< -2 . 53745E-6*Tf i In  >+.001341  903 

190  Dynvisc-< 3.234851 1E-7*(Tf iln“2  >>+<-. 0002 1 474487*Tf i In >+. 036 166792 

191  Kinvisc-Spvoi«Oynvisc 

192  Kf-( 1 -41818181 8E-3*Tf i In )+. 1868 

193  Pr-< 4.65706E-3* < Tf i ln*2  )  )+< -2. 922094«Tf i In >+463.33! 9 

194  Thetal -Power  1 /< .0220*Kf ) 

195  Ndt 1 ( I )-De It /Thetal 

196  Nd-Ndt 1(1) 

197  Temp=T ( I ) 

198  Theta-Qf lux*Lbar/Kf 

199  Ndt ( I  )»Del t/Theta 

200  Hx-Of lux/Dei t 

201  Grm«(  S»Beta*Qf  lux»  (  Lbar  )  “4  )/<  Kf  *h'invi  sc “2  ) 

202  Nux*Hx«Lbar/Kf 

203  Lgr< I  )-LGT(  Grn ) 

204  Lnui I  )*LGT< Nux  ) 

205  Lnn»Lnu( I ) 

206  Lg«Lgri I ) 

208  Xtol«<31 . 75-<  6 . 98S+< 1-1 6 )* 1 . 270 ) )/31 . 75 

209  PRINT  USING  “IX  ,D . OODD ,2X ,0 .000 , 1 X ,D. OOO  ,2X ,D.DDO ,5X ,D . ODD ,4X ,DD . DDD  ,  2X ,DD 
DD.D  ,  1 X  .DD.DDD ,DDDDD . 0 , 1 X  ,00. DO" ; Qcond  ,Xtol  ,Nd  ,Lnn  ,Lg  .Del t ,Hx ,Nux ,Grn ,Tenp 

210  NEXT  I 

211  PRINTER  IS  CRT 

217  I  |  !  )  n  M  I  M  i  i  i  i  m  i  I  |  |  m  |  |  |  m  |  !  |  M  M  (  t  (  (  I  I  |  (  i  it  j  ( 

218  i  ENTER  CONDUCTION  LOSSES  ! 

219  !  !  |  !  I  !  I  |  M  I  m  i  t{  {  {  (  {  |  i  )  |  j  {  j  i  m  i  m  I  i  }  |  j  |  i  j  m  i  it  i  i 

220  PRINT  "INPUT  CENTER  COLUMN  CONDUCTION  LOSSES  AS  COMPUTED  FROM  ELLPACK" 

221  INPUT  Qcond 

222  (  I  I  i  I  j  I  i  I  I  |  I  j  i  j  i  i  i  j  i  i  |  i  i  i  i  I  i  i  i  i  i  i  i  i  i  i  i  |  |  |  M  !  !  M  ! 

225  PRINTER  IS  701 

226  FOR  1-16  TO  30 

227  PRINT  "  “ 

228  PRINT  "  “ 

229  PRINT  “  SURFACE  #" ,1 

230  PRINT  "  “ 

231  PRINT  USING  "2X  .AAAAA  ,3X ,AAA  ,3X ,AAA ,3X  .AAAAAAAA  ,2X ,  AAAnAAAAA ,3X ,AAAA  ,4X  ,An 
,4X ,AAA  ,3X  ,AAAA  ,2X  ,AAAA" |A$  ,8$  ,CS  ,D$  ,E$  ,F$ ,GS  ,H*  .IS  , J* 

232  Qconv-Power2-0cond 

233  Of lux-Oconv/As 

234  Oelt-T( I )— Tamb 

235  Tfilm-273. 15+(T(I  HTamb >/2 

236  Beta-< -7. 944858E-8*<  Tf i lm“2 > >+( 5 .73S6479E-5«Tfi In >- .009781 0563 

237  Spvol-<4.69699E-9*(Tfiin"2  )  )+< -2 . 5374SE-6«Tf i In >+ .001 34 1 903 

238  Dynvisc-( 3.234851 1E-7*( Tf iln“2  ) >  +  <- .0002 1 474487»Tfi In  )+. 036 166792 

239  Kinvisc-Spvol*Dynvisc 

240  K f-( 1 .418!81818E-3*Tfiln  >+.1866 

241  Pr-< 4 . 65706E-3* ( T  f i ln“2 ) )  +  < -2 . 922094* Tf i In  1+463 . 33 1  9 

242  Theta1»Power2/< .0220»Kf ) 
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243  Ndt 1 < I  )*De It /Theta) 

244  Nd-NdtKI) 

245  Tenp-T  < I > 

246  Theta-Qf lux*Lbar/Kf 

247  Ndt< I  )«Delt/Theta 

248  Hx-Qf lux/Dei t 

249  Srn*<  G*Beta«Qflux«(Lbar)'>4)/(Kf»Kinvjsc“2  ) 

250  Nux-Hx*Lbar/Kf 

252  Lgr  ( I )-LGT  < Grm  ) 

253  Lnu< I  )-LGT(Nux  ) 

254  Lnn»Lnu< I > 

255  Lg-Lgr(I) 

256  Xtol-<  31  .75-<  6 . 985+<  1-16  >••!  .270  )  >/3l  .75 

257  PRINT  USING  " 1 X ,0 . DODD .2X ,0 .DOQ . 1 X ,D . ODD ,2X .0 .ODD ,SX ,D. ODD ,4X .DO. DDD  ,2X  ,DD 
00.0 ,1 X  .00.000  .DDODD.O , 1 X .OD.OQ" sOcond ,Xtol ,Nd  .Lnn ,Lg .Del t ,Hx ,Nux .Grm.Temp 

258  NEXT  I 

259  PRINTER  IS  CRT 

260  (Ml!!!!!  |  !!!!!!!(  |  MM!!!!!!  m  him  |  i  j  i  i  m  I;  t  m 

261  !  ENTER  CONDUCTION  LOSSES  ! 

262  l  I  l  I  I  i  l  l  !  I  I  I  ;  j  |  |  ;  |  I  !  |  |  !  I  I  !  {  (  i  i  |  I  I  (  f  j  i  j  {  !  j  I !  Ij  i  i  | 

263  PRINT  “INPUT  SIDE  COLUMN  CONDUCTION  LOSSES  AS  COMPUTED  FROM  ELLPACK* 

264  INPUT  Qcond 

265  (MM!!  I II  H  M  H  !*!!!!!  !*  H  M  H  !!!!  i  !!!!!!!!  !!  !* 

266  PRINTER  IS  701 

267  FOR  1-31  TO  45 

268  PRINT  “  * 

269  PRINT  *  “ 

270  PRINT  *  SURFACE  #“ .1 

271  PRINT  “  “ 

272  PRINT  USING  “2X .AAAAA ,3X .AAA ,3X , AAA ,3X .AAAAAAAA ,2X .AAAAAAAAA ,3X  . AAAA  ,4X  ,AA 
,4X .AAA.3X  .AAAA ,2X .AAAA" tA$ ,BS ,CS ,D$ ,ES ,F$ .6$ ,HS .1$ ,J$ 

273  Qconv-Pouerl-Qcond 

274  Qf  lux-Qconv/As 

275  Del t»T( I )-Tanb 

276  Tf jln-273. 15+<T< I )+Tanb >/2 

277  Beta-( -7 . 944858E-8* (  Tf  i lrn'  2  )  >+<  5 .7356479E-5*T f  i In  )-. 00978 1  0563 

27S  Spvol«(4.69699E-9*(Tfiln‘2  )  >+( -2 . 53745E-6*Tf i In )+ .00 1 34 1 903 

279  Dynvisc«< 3.234851 1 E-7* ( T f i ln'2 ) )  +  < - .0002 1 474487*Tf i In >+.036166792 

280  Kinvisc=5pvol*0ynvisc  < 

261  Kf*( 1 .41 818181 8E-3*Tfi In )+. 1866 

262  Pr-( 4. 65706E-3*(Tfiln“2  >  >+< -2 . 922094»Tf i In >+463 . 33 1 9 

263  Thet a  1 -Power  1 / ( . 0223»K f > 

264  Ndt 1 i I  )-Del t/Thetal 

285  Nd-Ndt 1(1) 

286  Tenp-T < I ) 

287  Theta-Qf lux*Lbar/Kf 

288  Ndt ( I >-Del t/Theta 

289  Hx-Qf lux/Dei t 

290  Grn-(  G*Beta»Qflux*(Lbar  ) "  4  )/(Kf*Kinvisc''2  > 

281  Nua-Hx  *Lbar/K  f 

292  Lgr( I >-LGT ( Grn ) 

293  Lnu< I >-LGT < Nux  ) 

294  Lnn-Lnu( I  ) 

295  Lg-Lgr < I > 

297  Xtol-<31.75-(6. 985+ ( 1-1 6  )* 1 . 270 > )/31 .75 

298  PRINT  USING  “ 1 X .D.DDOO ,2X ,0. DOD , IX ,D. DDD ,2X .D.DDD ,5X .D.DOD ,4X .DD.DDD ,2X ,DD 
DO.D  ,  I X  ,00. DDD, DDODD.O  ,  1 X ,0D .DD“ tOcond ,X  tol ,Nd  ,Lnn ,Lg .Del t  ,Hx ,Nux ,Grn .Tenp 

2S9  NEXT  I 

300  OUTPUT  ®Path3 i T<  *  > ,Ndt 1 < • >  .Ndt  <  *  >  ,Lnu( *  >  ,Lgr( •  ) 

301  PRINTER  IS  CRT 

302  ENO 


\ 

I 


l 
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